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ABOUT  THESE  PROCEEDINGS 


Due  to  the  large  number  of  papers  present"'.'!  at  the  Cbnference,  I  have 
divided  the  proceedings  into  two  manageable  volumes.  At  the  beginning  of 
each  volume  is  a  list  of  all  papers  contained  in  both  volumes,  in  actual 
order  of  presentation  at  the  Conference.  This  is  also  the  sequence  of  the 
published  papers  within  these  proceedings. 

For  the  sake  of  completeness,  both  volumes  contain  the  Attendees  List, 
Cbnference  Agenda,  Lists  of  Papers,  Conference  History,  Acknowledgments 
and  this  explanation. 

Every  paper  is  preceded  by  an  abstract  in  a  standard  format.  Some  papers 
may  also  have  the  original  abstract  included.  Further,  you  may  recall  the 
Q&A  session  we  had  at  the  end  of  each  presentation.  In  cases  where  a 
technical  Interchange  did  take  place,  the  questions  and  answers  are 
documented  at  the  end  of  each  pertinent  paper.  Every  paper  did  not  have 
a  Q&A  session,  and  1  have  included  all  Q&A  sheets  that  were  handed  to  me 
at  the  end  of  each  presentation.  Except  for  a  few  minor  editorial  changes, 
the  information  on  these  sheets  has  not  been  significantly  altered. 
Obviously,  these  sheets  are  as  "good”  as  the  inputs  you  provided. 

In  summary,  I  hope  the  above  explanation  was  helpful.  1  have  done  what  I 
consider  to  be  a  thorough  job  of  collecting  and  checking  all  the 
information  for  these  proceedings.  Errors  will  occur,  however,  and  while 
I  will  entertain  any  comments  and  criticisms  on  this  issue,  these 
proceedings  will  stand  as  published. 

Thank  you  for  your  participation,  and  your  patience! 


Capt  Vishnu  V.  Nevrekar  November  1987 

Earth  Sciences  Division 
USAF  Geophysics  Laboratory 
Hanscom  AFB,  MA  01731 


ABOUT  THE  GRAVITY  GRADIOMETRY  CONFERENCE 


The  First  Gravity  Gradiometry  Conference  was  held  at  the  Air  Force  Cambridge 
Research  Laboratory  (AFCRL,  now  AFGL)  in  1973.  Its  purpose  was  to  provide 
a  forum  to  evaluate  and  compare  the  efforts  of  three  vendors  (Charles  Stark 
Draper  Lab,  Hughes  Research  Lab  and  Bell  Aerospace  Textron)  in  still-emerging 
areas  of  gravity  gradiometry.  About  15  people  attended,  most  of  them  from  the 
companies  mentioned  above  or  the  Terrestrial  Sciences  Division  at  AFCRL.  In 
contrast,  the  1987  Conference  had  a  guest  list  of  70  plus  attendees,  with 
participation  from  academia  (foreign  and  domestic),  private  industry  and 
government.  The  papers  presented  were  not  restricted  to  gradiometry  alone. 
Indeed,  the  scope  of  this  annual  event  has  broadened  considerably  since  1973. 

With  the  exception  of  the  first  two  conferences,  all  the  others  have  been 
held  at  the  US  Air  Force  Academy  in  Colorado  Springs,  Colorado.  The  Geodesy 
and  Gravity  Branch  of  the  Earth  Sciences  Division  of  the  Air  Force  Geophysics 
Laboratory  (AFGL),  Hanscom  AFB,  Massachusetts,  has  always  organized  the  event, 
which  usually  takes  place  around  the  second  week  in  February.  This  trend  is 
expected  to  continue. 

If  you  are  not  already  on  our  mailing  list  and  would  like  to  attend  the 
1988  Conference,  or  if  you  have  any  questions,  please  write  to: 

Ms  Qaire  McCartney 
AFGL/LW 

Hanscom  AFB,  MA  01731 

Due  to  space  constraints,  we  restrict  the  size  of  our  Conferences  to  about  75 
people.  Attendance  will  generally  be  on  a  "first-come,  first-served  basis"  once 
the  completed  registration  forms  are  returned  to  us.  We  shall  mail  these  forms 
later  this  year. 

While  we  have  a  limited  number  of  copies  of  the  proceedings  for  non-attendees 
of  the  1987  Conference,  copies  of  proceedings  for  prior  years  are  not  available. 
Also,  we  appreciate  any  comments  or  suggestions  you  may  have  regarding  this 
document . 


ACKNOWLEDQIENTS 


We  couldn't  possibly  organize  a  conference  the  scope  and  size  of  our  forum 
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Conference . 
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acknowledge  her  assistance  in  handling  all  the  conference  arrangements, 
including  hotel  accomodation  for  the  attendees,  transportation  for  the 
Conference  and  NORAD/ CSOC  tours,  luncheons  during  the  conference  and  the 
"mixers"  later  in  the  evening,  all  of  which  were  set  up  with  great  skill 
and  professionalism. 
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at  the  complex  space  defense  environment. 

Next,  we  thank  all  the  speakers  for  taking  the  time  to  compile  and  present 
their  papers  for  the  benefit  of  the  Conference  attendees.  As  always,  the  broad 
mix  of  topics  went  a  long  way  towards  making  the  Conference  an  intellectually 
stimulating  event.  Indeed,  the  high  quality  of  the  research  material  presented 
"made"  this  Conference. 

Finally,  we  thank  Col  J.R.  Johnson,  then  Commander,  AFGL,  Dr  Donald  H.  Eckhardt , 
Director,  Earth  Sciences  Division  and  Dr  Thomas  P.  Rooney,  Giief,  Geodesy 
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Fifteenth  Gravity  Gradlometer  Conference 
United  States  Air  Force  Academy 
Colorado  Springs,  Colorado 

CONFERENCE  AGENDA 


Tuesday,  10  February  1987 

1900  -  2200  -  Pre- (bnference  Get-Together  at  Hilton  Inn 

Early  Registration 

Wednesday,  II  February  1987 

0700  -  Depart  Hilton  Inn  for  Fairchild  Hall 

0730  -  Registration  -  3rd  floor  Fairchild  Hall,  South  End 

0745  -  Welcome/Introduction  -  Copt  Terry  J.  Rindak 

0815  -  Presentation  by  Dr.  Georges  Balmlno  of  the  ONERA  (Office  National 
d' Etudes  et  de  Recherches  Aerospatlales) • 

"GRADIO  Project:  A  SGG  Mission  Based  on  Microaccelerometers” 

0845  -  Presentation  by  Dr.  G.  Ian  Moore  of  the  University  of  Queensland. 

"A  Mercury  Manometer  Gravity  Gradlometer” 

0900  -  Presentation  by  Mr.  Ernest  H.  Metzger  of  Bell  Aerospace  Textron. 

"Bell  Aerospace  Gravity  Gradlometer  Survey  System  (GGSS)  -  Program 
Review" 

0925  -  Presentation  by  Dr.  Frank  J-  van  Kann  of  the  University  of  Western 
Australia . 

"A  Prototype  Superconducting  Gravity  Gradlometer  for  Geophysical 
Exploration” 

0952  -  Presentation  by  Dr.  Warren  G.  Heller  of  The  Analytic  Sciences  Corp. 

"Gravity  Gradlometer  Survey  System  (GGSS)  Data  Processing  and  Data 
Use” 

1016  -  Break 

1035  -  Presentation  by  Mr.  A1  Jircitano  of  Bell  Aerospace  Textron. 
"Self-Gradient  Calibration  of  the  GGSS  in  a  C-130  Aircraft" 


1058  -  Presentation  by  Dr.  Sam  C.  Bose  of  Applied  Sciences  Analytics,  Inc • 

"Gravity  Gradiometer  Data  Procescini^  Using  the  Karhunen-Loeve 
Method" 

1120  -  Presentation  by  Mr.  David  M.  Gleason  of  the  Air  Force  Geophysics 
Laboratory. 

"Numerically  Deriving  the  Kernels  of  an  Integral  Predictor  Yielding 
Surface  Gravity  Disturbance  Components  from  Airborne  Gradient  Data" 

1130  -  Presentation  by  Mr.  A1  Jircitano  of  Bell  Aerospace  Textron. 

"Stage  II  Simulation  Results  Using  the  NSWC  Synthetic  Gravity 
Field" 

1150  -  Depart  Fairchild  Hall  for  USAFA  Noncommissioned  Officers' 

(NOD)  Qub 

1200  -  Lunch  -  USAFA  NOD  Qub 

1245  -  Depart  NCD  Qub  for  Fairchild  Hall 

1330  -  Presentation  by  Dr.  Richard  H.  Rapp  of  Ohio  State  University. 

"Gradient  Information  in  New  High  Degree  Spherical  Harmonic 
Expansions " 

1354  -  Presentation  by  Mr.  John  J.  Graham  of  the  Defense  Mapping  Agency 
Aerospace  Cfenter. 

"The  Effect  of  Topography  on  Airborne  Gravity  Gradiometer  Data" 

1357  -  Presentation  by  Mr.  Mike  Sideris  of  the  University  of  Oblgary. 

"Effect  of  Terrain  Representation,  Grid  Spacing,  and  Flight  Altitude 
on  Topographic  Cbrrectlons  for  Airborne  Gradlometry" 

1417  -  Presentation  by  Dr.  Rene  Forsberg  of  Geodetic  Institute  (Denmark) 
(currently  at  the  University  of  Cblgary,  Cbnada). 

"Topographic  Effects  in  Airborne  Gravity  Gradlometry" 

1434  -  Presentation  by  Dr.  Alan  H.  Zorn  of  Dynamics  Research  Corporation. 

"Observability  of  Laplace's  Equation  Using  a  Torsion-Type 
Gravity  Gradiometer” 

1510  -  Break 

1530  -  Presentation  by  Dr.  Citl  Bowin  of  Woods  Hole  Oceanographic  Institute. 

"Ratios  of  Gravity  Gradient,  Gravity,  and  Geoid  for  Determination  of 
Qrustal  Structure" 


1550  -  Presentation  by  Dr.  Rene  Forsberg  of  Geodetic  Institute  (Denmark). 

"Combining  Gravity  Gradiometry  with  Other  Exploration 
Methods  for  Geophysical  Prospecting" 

1600  -  Presentation  by  Dr.  Rene  Forsberg  of  Geodetic  Institute  (Denmark). 

"Computation  of  the  Gravity  Vector  from  Torsion  Balance  Data  in  S. 
Ohio" 

1615  -  Presentation  by  Dr.  Hans  Baussus  von  Luetzow  of  the  U.S.  Army 
Engineer  Topographic  Laboratories. 

"Estimation  of  Gravity  Vector  Components  from  Bell  Gravity  Gradiometer 
and  Auxiliary  Data  under  Gbnsideration  of  Topography  and  .Associated 
Analytical  Upward  Continuation  Aspects" 

1635  -  Depart  Fairchild  Hall  for  the  Hilton  Inn 

1700  -  Reception  -  Hilton  Inn 

Thursday,  12  February  1987 

0700  -  Depart  Hilton  Inn  for  Faircnild  Hall 

0755  -  Presentation  by  Dr.  M.  VoL  Moody  of  the  University  of  Maryland. 

"Development  of  A  Three-Axis  Superconducting  Gravity  Gradiometer 
and  a  Six-Axis  Superconducting  Accelerometer" 

0835  -  Presentation  by  Dr.  Bahram  Mashhoon  of  the  University  of  Missourl- 
Golumbia . 

"The  Gravitational  Magnetic  Field  of  the  Earth  and  the  Possibility 
of  Measuring  It  Using  an  Orbiting  Gravity  Gradiometer" 

0905  -  Presentation  by  Dr.  Ho  Jung  Palk  of  the  University  of  Maryland. 

"Tests  of  General  Relativity  in  Earth  Orbit  Using  a  Superconducting 
Gravity  Gradiometer" 

0928  -  Presentation  by  Dr.  Dave  Sonnabend  of  Jet  Propulsion  Laboratory. 

"Magnetic  Isolation  -  Qoslng  the  Loop" 

0941  -  Presentation  by  Dr.  Dan  Long  of  Eastern  Washington  University. 
"Laboratory  G(R)  Experiment  -  Progress  Report" 


1004  -  Break 


1030  -  Cheyenne  Mountain  Complex  Overview 
Briefing  by  Maj  Bill  Cbrver,  USAF 
(Chief,  NORAD  Presentations  Division). 

1110  -  Form  Groups  A  &  B 

1115  -  Depart  Fairchild  Hall  for  USAFA  NOD  Club 
1130  -  Lunch  -  USAFA  NGD  Qub 

1200  -  Depart  USAF  Academy  for  Falcon  Air  Force  Station 

1245  -  Arrive  Falcon  APS  for  briefing  on  2nd  Space  Wing 

Tour  of  the  Cbnsolidated  Space  Operations  Center  (CSOC) 

(Group  A) 

1415  -  Depart  CSOC  for  Cheyenne  Mountain  Complex  (CMC) 

1500  -  Arrive  CMC 

1505  -  Security  in-processing  and  process  through  metal  detector 
1525  -  Travel 

1530  -  Tour  NORAD  Command  Post 
Tour  Industrial  Area 

1620  -  Travel/question  and  answer  session 
1630  -  Depart  for  Hilton  Inn 

1715  -  Arrive  Hilton  Inn 

(Group  B) 

1415  -  Depart  USAF  Academy  for  Peterson  Air  Force  Base  (AFB) 

1445  -  Arrive  Peterson  AFB  museum 

1600  -  Depart  Peterson  AFB  for  Hilton  Inn 

1630  -  Arrive  Hilton  Inn 

Ftiday,  13  February  1987 


0800  -  Tour  of  JILA,  Boulder,  Qilorado 


Papers  Included  in  VOLUME  I  of  the  Oanference  Proceedings 


1.  *Dr.  Georges  Balmlno,  C-N.E.S. /Bureau  Gravlmetrlque  International,  France 
Dr.  Alain  Bernard,  ONERA  (Office  National  d'Etudes  et  de  Recherches 

Aerospatlales ,  Ftance) 

Dr.  Pierre  Touboul,  ONERA,  France 

"GRADIO  Project  :  A  SGG  Mission  Based  on 
Mlcroaccelerometers" 


2.  *Dr.  G.  Ian  Moore,  University  of  Queensland,  Australia 
Dr.  Frank  D.  Stacey,  University  of  Queensland,  Australia 
Dr.  Gary  J.  Tuck,  University  of  Queensland,  Australia 
Dr.  Barry  D.  Goodwin,  University  of  Queensland,  Australia 

"A  Mercury  Manometer  Gravity  Gradiometer" 


3.  Mr.  Louis  L.  Pfohl,  Bell  Aerospace  Textron 
*Mt.  Ernest  Metzger,  Bell  Aerospace  Textron 

"Bell  Aerospace  Gravity  Gradiometer  Survey 
System  (GGSS)  -  Program  Review" 


4.  *Dr.  Frank  J.  van  Kann,  et  al.  University  of  Western  Australia 

"A  Prototype  Superconducting  Gravity 
Gradiometer  for  Geophysical  Exploration" 


5.  *Dr.  Warren  G.  Heller,  The  Analytic  Sciences  Corporation 

"Gravity  Gradiometer  Survey  Systtu  (GGSS) 
Data  Processing  and  Data  Use” 


6.  Dr.  W.  John  Hutcheson,  Bell  Aerospace  Textron 

(Paper  presented  by  Mr.  Al  Jircitano  of  Bell  Aerospace  Textron) 

"Self-Gradient  Calibration  of  the  GGSS 
in  a  0-130  Aircraft" 


7.  *0r.  Sam  C.  Bose,  Applied  Science  Analytics,  Inc 

Mr.  Glenn  E.  Thobe,  Applied  Science  Analytics,  Inc 

"Gravity  Gradiometer  Data  Processing  Using 
the  Karhunen-Loeve  Method" 


*  Denotes  Speaker  at  Co'^iference 


8.  *Mr .  David  M.  Gleason,  Air  Force  Geophysics  Laboratory 


"Numerically  Deriving  the  Kernels  of  an  Integral 
Predictor  Yielding  Surface  Gravity  Disturbance 
Cbmponents  from  Airborne  Gradient  Data" 


9.  Dr.  W.  John  Hutcheson,  Bell  Aerospace  Textron 

(Paper  presented  by  Mr.  A1  Jircitano  of  Bell  Aerospace  Textron) 

"Stage  II  Simulation  Results  Using 
the  NSWC  Synthetic  Gravity  Field" 


10.  *Dr.  Richard  H.  Rapp,  Ohio  State  University 

"Gradient  Information  in  New  High 
Degree  Spherical  Harmonic  Expansions 


*  Denotes  Speaker  at  Conference 


Papers  included  In  VOLUME  II  of  the  Conference  Proceedings 


1.  *Mr.  John  J.  Graham,  Defense  Mapping  Agency  Aerospace  Center 

Mr.  Joseph  L.  Toohey,  Defense  Mapping  Agency  Aerospace  Center 

"The  Effect  of  Topography  on  Airborne 
Gravity  Gradlometer  Data" 

2.  Dr.  Klaus-Peter  Schwarz,  University  of  Chlgary,  Chnada 
*Mr.  M.G.  Slderls,  University  of  Calgary,  Canada 

Dr.  I.N.  Tziavos,  University  of  Chlgary,  Canada 
(Dr.  Tziavos  on  leave  from  the  University  of  Thessaloniki,  Greece) 

"Effect  of  Terrain  Representation,  Grid  Spacing,  and 
Flight  Altitude  on  Topographic  Corrections  for 
Airborne  Gradiometry" 


3.  *Dr.  Rene  Forsberg,  Geodaetisk  Institut,  Denmark 

"Topographic  Effects  in  Airborne  Gravity  Gradiometry" 


4.  *Dr.  Alan  H.  Zorn,  Dynamics  Research  Corporation 

"Observability  of  Laplace's  Equation  Using 
a  Torsion-Type  Gravity  Gradlometer” 


5.  *Dr.  (irl  Bowin,  Woods  Hole  Oceanographic  Institute 

"Ratios  of  Gravity  Gradient,  Gravity,  and  Geoid 
for  Determination  of  Qrustal  Structure" 


6.  Dr.  Anthony  A.  Vassillou,  University  of  Calgary,  Canada 

(Paper  presented  by  Dr.  Rene  Forsberg,  Geodaetisk  Institut,  Denmark) 

"Combining  Gravity  Gradiometry  with  other 
Exploration  Methods  for  Geophysical  Prospecting 


7.  Dr.  D.  Arabelos,  University  of  Thessaloniki,  Greece 

Mr.  Christian  Tschernlng,  Geodaetisk  Institut,  Denmark 

(Paper  presented  by  Dr.  Rene  Forsberg,  Geodaetisk  Institut,  Denmark) 

"Computation  of  the  Gravity  Vector  from  Torsion 
Balance  Data  In  Southern  Ohio" 


*  Denotes  Speaker  at  (bnference 


8.  *Dr.  Hans  Baussus  von  Luetzow,  US  Army  Engineer  Topographic  Laboratory 

"Estimation  of  Gravity  Vector  Components  from  Bell  Gravity 
Gradlometer  and  Auxiliary  Data  under  Consideration  of 
Topography  and  Assoclat-u  Analytical  Upward  Continuation 
Aspects" 

9.  Dr.  H.  A.  Gian,  University  of  Maryland 
Dr.  Q.  Kong,  University  of  Maryland 

*Dr.  M.  Vol  Moody,  University  of  Maryland 
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ABSTRACT 

The  reduction  and  conversion  of  airborne  gravity  gradlometer  data  to  ground 
level  estimates  of  the  gravity  disturbance  vector  Is  currently  of  considerable 
Interest  In  support  of  short  wavelength  gravity  modeling.  A  pressing  problem 
Is  the  need  for  an  accurate  procedure  for  the  downward  continuation  of  data 
acquired  at  altitude  by  the  airborne  Gravity  Gradlometer  Survey  System  (GGSS). 

As  part  of  ongoing  Investigations,  a  prism  method  has  been  used  to  calculate  the 
effect  of  topography  on  the  gravity  disturbance  vector  and  the  five  Independent 
second-order  gravity  gradients.  Calculations  of  the  contribution  of  topography 
to  the  magnitudes  of  these  gravimetric  parameters  were  made  at  both  surface  and 
elevated  points  In  the  Wichita  Mountains  of  Oklahoma.  Cbmputatlons  were  made 
utilizing  Digital  Terrain  Elevation  Data  (DTED)  with  an  assumed  constant  density 
of  2.67  grams/centlmeters^  for  the  topographic  masses.  Results  are  presented 
which  reflect  the  use  of  DTED  sets  of  different  horizontal  extent  and  grid 


Interval. 
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ABSTRACT 


The  reduction  and  conversion  of  airborne  gravity  gradiometer  data  to  ground 
level  estimates  of  the  gravity  disturbance  vector  is  currently  of  considerable 
interest  in  support  of  short  wavelength  gravity  modeling.  A  pressing  problem 
is  the  need  for  an  accurate  procedure  for  the  downward  continuation  of  data 
acquired  at  altitude  by  the  Airborne  Gravity  Gradiometer  Survey  System  (GGSS). 
As  part  of  ongoing  investigations,  a  prism  method  has  been  used  to  calculate 
the  effect  of  topography  on  the  gravity  disturbance  vector  and  the  five  inde¬ 
pendent  second-order  gravity  gradients.  Calculations  of  the  contribution  of 
topography  to  the  magnitudes  of  these  gravimetric  parameters  were  made  at  both 
surface  and  elevated  points  in  the  Wichita  Mountains  of  Oklahoma.  Computations 
were  made  utilizing  Digital  Terrain  Elevation  Data  (DTED)  with  an  assumed  con¬ 
stant  density  of  2.67  grams/centimeters'^  for  the  topographic  masses.  Results 
are  presented  which  reflect  the  use  of  DTED  sets  of  different  horizontal  extent 
and  grid  interval . 


THE  EFFECT  OF  TOPOGRAPHY  ON  AIRBORNE  GRAVITY  GRADIOMETER  DATA 


I.  INTRODUCTION 

The  field  testing  of  the  airborne  Gravity  Gradiometer  Survey  System  (GGSS), 
being  built  by  Bell  Aerospace/Textron  for  the  Defense  Mapping  Agency,  is 
scheduled  to  commence  before  midyear.  Therefore,  methods  for  validation  of 
the  system's  ability  to  map  the  local  gravity  field  to  high  detail  through 
gravity  gradient  measurements  and  their  subsequent  downward  continuation/ 
conversion  into  surface  gravity  disturbance  components  are  of  pressing 
interest.  Validation  can  in-part  be  accomplished  by  estimating  the  influ¬ 
ence  of  local  topography  on  the  radial  disturbance,  the  deflection  compo¬ 
nents,  and  the  second-order  gravity  gradients  at  both  surface  and  aloft 
stations.  This  report  summarizes  the  results  of  computing  topographic 
terrain  effects  from  the  topography  above  mean  sea  level.  The  terrain 
effects  are  calculated  on  the  basis  of  homogeneous  rectangular  prisms  which 
model  the  terrain  masses  with  an  assumed  constant  density. 

The  surface  computation  points  coincide  with  two  astro-geodetic  stations 
located  in  the  Wichita  Mountains  of  Oklahoma.  One  is  near  Sunset  Peak  and 
the  other  is  on  Mount  Scott  as  shown  in  Figure  1.  Terrain  effects  were 
computed  at  these  surface  stations  and  at  points  directly  overhead  at 
altitude  5000  feet  above  the  geoid.  Digital  Terrain  Elevation  Data  (DTED) 
supplied  the  topographic  model  needed  to  compute  the  terrain  effects.  The 
main  goal  of  the  study  was  to  determine  which  DTED  field  should  make  up  an 
inner  grid  zone  and  to  what  radial  extent  outward.  There  was  also  a  need 
to  establish  the  coarsest  DTED  representation  permissible  for  the  outer 
zone  and  its  span  of  coverage  for  adequate  modeling  of  the  local  terrain 
effects  on  selective  gravimetric  quantities. 


II.  DISCUSSION 
a.  Objective 

A  major  part  of  the  short  wavelength  variation  of  the  gravity  field  in  a 
local  area  is  due  to  topography.  Therefore,  especially  in  mountainous 
areas,  one  would  anticipate  the  need  for  the  best  available  topographic 
elevation  data  in  the  immediate  vicinity  of  a  computation  point.  At  some 
further  distance  beyond  this,  a  coarser  terrain  representation  could  be 
utilized  to  make  the  computational  process  more  efficient  with  a  minimal 
effect  on  results.  With  this  premise,  terrain  effects  were  computed  for 
the  stations  shown  in  Figure  1  to  establish  the  radius  of  an  inner  zone 
for  3"  point  DTED  which  is  our  finest  grain  DTED.  The  computations  also 
allow  the  determination  of  the  largest  mean  DTED  representation  that  is 
permissible  for  an  outer  zone  and  its  span  of  coverage  for  adequate 
modeling  of  the  local  terrain  effects  on  selective  gravimetric  quantities. 
These  quantities  include  the  radial  gravity  disturbance,  the  deflection 
of  the  vertical  components,  and  the  second-order  gravity  gradients. 

Inner  zone  modeling  with  mean  DTED  will  indicate  whether  there  is  a  need 
for  the  exclusive  use  of  point  DTED  for  this  area.  Insight  into  upward 
continuation  effects  is  afforded  with  the  inclusion  of  aloft  computation 
points. 


b.  Method  of  Analysis 


A  modified  version  of  the  Rene  Forsberg  prism  program  was  run  to  compute 
the  topographic  terrain  effects.  T!  ^  program  is  discussed  in  Reference  1. 
The  right-handed  coordinate  system  was  centered  at  each  computation  sta¬ 
tion  with  X  pointing  east,  Y  pointing  north,  and  Z  pointing  up.  Prisms 
were  formed  from  the  geoid  base  up  to  the  DIED  topography  of  an  assumed 
density  of  2.67  gm/cm^  as  illustrated  in  Figure  2.  The  integration  was 
performed  numerically  using  the  distances  from  station  coordinates  to  the 
eight  corners  of  each  prism.  The  surface  station  height  was  offset  1  cm 
to  avoid  the  central  prism  from  being  automatically  eliminated  from  the 
program  and  to  avoid  the  station  from  being  located  within  the  central 
prism  boundaries. 

c.  DTED  Fields  Used 

A  3"  point  DTEO  field  was  built  around  the  stations  from  the  same  data 
that  OMAAC  sent  to  Bell  Aerospace/Textron.  Also  created  were  9",  12", 

15",  30",  r,  and  3'  mean  fields  from  the  3"  point  data.  When  using 
both  an  inner  field  and  an  outer  field  around  a  station,  the  prisms  must 
properly  fit  together.  If  the  inner  field  is  the  3"  DTED,  then  the  outer 
mean  field  must  be  a  multiple  of  6"  for  the  prisms  to  fit  together. 

d.  Formulation 

Terrain  mass  modeling  is  accomplished  with  homogeneous  rectangular  prisms 
in  the  calculations.  Gravitational  formulas  for  such  prisms  are  known 
from  MacMillan's  work  on  potential  theory  in  Reference  2.  The  prism 
dimensions  for  this  study  were  controlled  by  the  description  of  the  topo¬ 
graphy.  Figure  3  shows  the  indefinite  integral  solutions  that  are  used 
to  compute  exact  values  for  ten  gravimetric  quantities  at  the  computation 
point  P  in  Figure  2  for  a  single  prism.  When  the  separation  between  a 
prism  and  the  computation  point  permits,  approximate  prism  formulations 
are  used  instead  of  the  exact  ones.  Figure  4  gives  details  on  how  a 
series  expansion  for  a  prism's  potential  may  be  derived  by  formulating 
the  reciprocal  distance  "r"  as  an  infinite  series  in  terms  of  the  Legendre 
polynomials.  Using  only  the  first  few  terms  of  the  potential  series,  all 
desired  gravity  quantities  may  be  found  by  simple  differentiation.  The 
exact  formulation  is  normally  used  for  the  central  prism  while  approximate 
formulation  is  used  for  all  remaining  prisms  to  obtain  the  desired  gravity 
quantities. 


III.  RESULTS 

a.  Tabular  Output 

Tables  1  and  2  show  the  effects  of  topgraphic  modeling  by  different  DTED 
fields  within  a  12'  radius  about  the  two  surface  stations.  The  various 
terrain  fields  used  were  3",  9",  12",  15",  30",  1',  and  3'  DTED.  The 
same  investigation  is  repeated  in  Tables  3  and  4  for  the  computation 
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points  at  elevation  5000  feet  above  the  geoid  and  directly  overhead  the 
ground  stations.  Terrain  effect  variations  due  to  unit  step  increases 
over  the  interval,  6'  to  13',  in  the  radius  of  a  central  zone  of  3"  point 
OTEO  are  shown  in  Tables  5  and  6  for  the  ground  level  stations.  In  the 
remaining  tables,  terrain  effects  are  accounted  for  by  employing  the  3“ 
point  DTEO  in  an  inner  zone  and  one  of  the  mean  DTED  fields  in  an  outer 
zone.  It  waw  discovered  that  the  prism  program  requires  the  outer  grid 
field  to  be  a  multiple  of  6"  for  the  prisms  to  properly  fit  together. 

This  is  the  same  as  saying  that  if  one  extends  the  3"  field,  then  the 
outer  mean  field  grid  points  must  coincide  with  a  3"  point.  Examples  of 
outer  fields  would  then  be  12",  30",  1',  and  3'. 

Tables  7  and  8  indicate  variation  in  terrain  effects  caused  by  mean  DTED 
representations  in  an  outer  zone  of  radius  30'  and  an  inner  zone  radius 
of  12*.  The  means  considered  for  the  outer  zone  terrain  modeling  were 
12",  30",  r  and  3'  DTEO.  Each  table  displays  results  for  one  of  the 
surface  points.  Sunset  Peak  or  Mount  Scott,  and  the  related  overhead 
point.  Tables  9  and  10  reflect  the  usage  of  1 '  mean  DTED  in  two  expanded 
outer  zones,  one  of  radius  2^  and  the  other  2.5°  with  an  inner  zone  radius 
of  12'.  These  are  again  composite  tables  for  a  pair  of  ground/aloft 

stations.  For  the  ground  level  stations.  Tables  11  and  12  indicate  ter¬ 

rain  variations  due  to  inner  zone  modeling  with  different  radial  extents 
over  the  interval  1'  to  10'  while  the  outer  1'  DTED  zone  was  extended  out 
to  a  radius  of  2.5^.  Tables  13  and  14  are  composite  tables  for  a  pair 
of  ground/aloft  stations  where  the  terrain  modeling  was  accomplished  by 
an  inner  region  of  3"  point  DTEO  and  an  outer  zone  of  3'  mean  DTED.  The 

radii  of  the  inner  and  outer  regions  were  5'  and  2.5°,  respectively,  in 

one  case  while  6'  and  3^  in  another. 

b.  Analysis 

From  Tables  1-4  the  pronounced  decline  in  the  ground  elevation  at  Mount 
Scott  with  the  increasingly  smoother  DTED  modeling  in  a  12'  central  zone 
indicates  steeper  terrain  here  than  near  Sunset  Peak.  A  200  meter  change 
in  elevation  occurred  when  the  1'  OTEO  field  was  used  instead  of  the  best 
3"  point  DTEO  representation.  Tables  1  and  2  clearly  show  the  need  for 
the  best  DTEO  modeling  in  the  immediate  vicinity  of  each  ground  station 
as  large  gradient  variations  resulted  even  with  the  9"  mean  DTED  field. 

The  same  requirement  applies  to  the  aloft  stations  even  though  the  vari¬ 
ations  are  not  as  striking  due  to  attenuation  effects  with  altitude. 

Tables  5-6  show  diminishing  contributions  from  successive  additions  of  1' 
bandwidths  of  3"  point  DTED  to  a  central  zone  with  an  initial  radius  of 
6'.  The  aggregate  contribution  of  the  seven  1'  bands  is  between  -23  E 
and  -25  E  on  Tzz  and  11-13  E  on  Txx  and  Tyy,  These  tables  also  show  that 
the  off-diagonal  gradients  experienced  magnitude  changes  of  less  than  1  E. 

Later  tables  will  indicate  an  improbable  need  for  continuance  of  best 
terrain  modeling  beyond  the  6'  radius.  From  Tables  7-8  the  inclusion  of 
an  outer  band,  extending  from  12*  to  30',  has  an  impact  of  approximately 
-12  E  on  Tzz,  5-7  E  on  Txx  and  Tyy,  and  between  -0.02- E  and  0.7  E  on  the 
off-diagonal  gradients  irrespective  of  terrain  modeling  by  12",  30",  1', 
or  3'  mean  OTEO.  This  beckons  the  use  of  either  the  1'  or  3'  mean  DTEO 
field  for  outer  zone  modeling.  Tables  9  and  10  show  that  outer  zone 


modeling  with  the  1'  mean  DIED  field  must  extend  to  2.5°  for  the  diagonal 
gradients  to  exhibit  satisfactory  convergence.  For  all  cases,  the  span  of 
1'  mean  DIED  from  30'  to  2.5°  contributes  between  -6  E  and  -7  E  to  Tzz 
with  approximately  -0.5  E  of  that  '.^_'ning  from  the  band  2°-2.5°.  This  same 
half-degree  width  band  accounts  for  less  than  0.4  E  to  the  Txx  and  Tyy 
values.  Even  less  contribution  to  the  results  would  be  expected  from  any 
band  beyond  2.5°  due  to  the  nature  of  the  kernel  functions.  With  entries 
in  Tables  9-10  as  benchmarks,  results  in  Tables  11-12  show  that  the  inner 
grid  radius  for  the  3“  point  DTEO  may  be  reduced  to  at  least  5'  or  6'  for 
a  2.5°  outer  radius  with  the  gradient  changes  being  0.24  E  or  less  in 
absolute  value.  Comparing  Tables  13-14  with  Tables  11-12  shows  that  outer 
zone  modeling,  from  5'  to  2.5°  with  3'  mean  DTED  instead  of  1'  mean  OTED, 
causes  gradient  changes  with  magnitudes  no  greater  than  0.41  E.  Selection 
of  the  larger  mean  field  would  reduce  computer  processing  time  in  that 
fewer  prisms  are  needed.  It  should  be  noted  that  each  table  contains  re¬ 
sults  for  the  gravity  disturbance  and  deflection  of  the  vertical  components 
that  the  reader  may  examine.  The  deflection  components  were  more  affected 
by  distant  topographic  masses  than  the  radial  gravity  disturbance. 

When  an  outer  grid  is  used,  the  inner  grid  values  are  extended  further 
than  they  were  in  Tables  1  through  4.  For  Instance  when  the  3'  field  was 
the  outer  grid  and  the  inner  radius  was  5'  as  in  Table  13,  the  portion  of 
fine  grid  used  in  computations  agrees  with  that  of  the  8'  radius  of  Table 
5.  Another  example  is  when  the  1'  field  was  the  outer  grid  and  the  inner 
radius  was  6'  as  in  Table  11.  This  agrees  with  the  portion  of  fine  grid 
used  in  computations  for  the  7'  radius  of  Table  5.  The  inner  field  is 
then  extended  by  one  unit  of  the  outer  grid  spacing  to  form  prisms  that 
fit  next  to  the  prisms  from  the  outer  field.  This  explains  why  the  outer 
grid  must  be  a  multiple  of  6"  for  the  prisms  to  fit  together.  This  also 
explains  the  increase  in  number  of  prisms  for  Tables  7  and  8.  The  3" 
field  was  extended  further  out  by  1'  and  3'  which  caused  an  increase  in 
number  of  prisms.  However,  as  one  increases  the  outer  radius  to  3°,  a 
decrease  in  number  of  prisms  would  be  seen  by  using  a  3'  outer  mean  field 
rather  than  an  outer  field  of  higher  density.  Computer  time  is  thus 
reduced  by  using  the  3'  outer  field. 


IV.  CONCLUSIONS  AND  FUTURE  CONSIDERATIONS 

The  inner  field  must  be  the  3"  point  DTED  with  an  inner  radius  of  at  least 
5'.  The  outer  field  may  be  3'  mean  OTEO  with  an  outer  radius  of  at  least 
?.5°.  If  the  field  is  big  enough  it  would  be  advisable  to  use  an  outer 
radius  of  3°  and  an  inner  radius  of  6'  in  general.  This  would  allow 
convergence  of  the  topographic  terrain  effects  within  acceptable  limits. 
The  3"  inner  grid  contributes  the  major  portions  of  the  gradient  values 
which  means  that  the  3"  elevations  near  the  station  coordinates  must  be 
as  accurate  as  possible.  An  inner  grid  that  is  a  less  dense  mean  field 
would  not  yield  the  proper  elevations  near  the  station.  The  outer  grid 
may  be  thought  of  as  fine-tuning  the  values  until  convergence  is  achieved 
within  acceptable  limits.  The  3"  OTED  tapes  are  written  in  3'  blocks 
which  makes  it  convenient  to  form  3'  mean  DTED  fields  for  the  outer  grid. 
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One  must  consider  the  computer  time  saved  by  using  3'  mean  DIED  fields 
for  the  outer  grid  as  large  numbers  of  stations  are  computed.  The  number 
of  prisms  is  cut  down  without  sacrificing  accuracy.  About  90  seconds  of 
CPU  time  is  added  per  station  with  3'  mean  DIED  out  to  a  radius  of  3°  and 
3"  point  DTEO  out  to  a  radius  of  6'.  The  radial  gravity  disturbance  and 
second-order  gradients  converged  within  desired  limits.  The  deflection 
components  will  not  converge  as  the  outer  grid  radial  extent  is  increased. 
Conversion  of  the  deflections  into  the  other  two  disturbance  components 
does  not  give  near  the  magnitude  of  the  radial  disturbance.  We  see  the 
high  frequency  nature  of  the  second-order  gradients  by  their  rapid  con¬ 
vergence.  The  much  lower  frequency  nature  of  the  deflection  components 
will  not  allow  convergence  in  the  computation  of  the  terrain  effects. 

In  the  future  it  would  be  desirable  to  augment  the  analysis  with  more 
stations  in  the  Wichita  mountains.  The  spectral  characteristics  and 
covariance  funcions  for  both  local  topography  and  terrain  effects  should 
be  defined.  Terrain  effect  computation  by  the  FFT  method  may  be  investi¬ 
gated  in  the  future  for  reduction  of  computer  processing  time.  Comparisons 
of  G6SS  gradiometry  data  against  gradients  derived  from  DTED  fields  will 
need  to  be  made.  It  would  be  desirable  to  incorporate  error  propagation 
into  the  topographic  terrain  effects  programs  as  errors  in  the  OTED  and 
density  assumption  become  better  defined.  A  major  goal  is  to  be  able  to 
add  long  and  intermediate  wavelength  information  to  terrain  effects  for 
prediction  of  more  accurate  gravimetric  quantities. 
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FIGURE  3.  Exact  Prlsn  FornuUcion  for  Different 


Cravimctric  Quantities 
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SUNSET  PEAK  EFFECTS  USING  ONLV  INNER  RADIUS  OF  12' 
AND  STATION  ELEVATION  LOCATED  ON  TERRAIN  SURFACE 


CO 

m 

o 

CD 

H* 

• 

o 

00 

CM 

rH 

VD 

CO 

n  in 

CO  cri 

n 

ID 

w 

o 

m  CO 

CM 

0>  CM  O  (M  cn  rH 

ocj 

(N 

C^l 

rH 

M 

rH 

iH  in  rH  o 

Qj 

CM 

>1 

• 

•  •  •  •  • 

H 

iH 

o  o  o  o  o 

■ 

1 

till 

(M 

CM 

in  o 

in 

o 

iH  (M 

CM 

CM 

H 

w 

• 

• 

•  • 

• 

• 

U 

U  rM 

rH  iH 

rH 

rH 

CM  cn  o  n 

(U 

1 

1 

1  1 

1 

I 

N 

n 

rH  o  rH  O 

in 

X 

• 

H 

(M 

1 

0 

0 

0 

0 

M 

u  cn 

wO  ^  "d* 

H*  C^l 

1 

1  i 

CO 

ij  ^ 

m  m 

• 

• 

•  • 

• 

• 

in 

rH 

r-H 

rH 

iH  fH 

o 

1 

1 

1  1 

1 

1 

^  iH 

00  ^  00  oi 

>1^  in 

00  iH  CD  c^ 

X  0 

• 

*  •  •  •  • 

Eh  (u  rH 

in  n  CM  vD  o 

in 

f— < 

CM  ro 

rH 

CM 

n 

rH  iH  rH 

r-< 

cn  o 

r>.  CD 

c 

o 

flj 

• 

# 

•  • 

• 

• 

•H 

a 

tnn 

m 

n 

CH 

n 

S  LD 

in 

in  in 

in 

in 

in 

cn  CM  in  CM  in 

w 

N  X)  ^ 

in  CD  n  rH 

N  ’H 

• 

EH  C  CM 

CM  iH  rv  in 

> 

• 

D 

t:J«  ID  n  CM 

uJ 

• 

u 

n 

in 

r>*  00 

in 

CO 

CO 

CO 

CO  00 

CD 

CD 

H* 

rH 

rH  in  1.0  ^  o 

in 

>1 

rH 

a>  00  cn  o 

a 

Sj 

{>1 

• 

aj 

H 

H*  cn  rH  m  CM 

jj 

•H 

1  rH  m  CM  rH 

> 

(U 

1 

till 

u 

£ 

CjJ 

n 

in 

CO 

in 

CO 

CO 

CD  00 

00 

CD 

CD  00  CD  00  in 

'Z 

•d* 

H* 

X 

CD  vD  QO  r*^  rH 

H 

X 

• 

•  •  •  •  • 

CO 

Eh 

rH  CM  n  in  n  ro 

rH  n  ^  rH  H  H 

1  I  1  t  1 

Q 

D 

uJ 

z 

• 

s  s 

m 

z 

U 

ca 

m 

cr» 

CM  in 

b 

o 

U 

n  cn  CM  in  o  o 

M 

rH  i-H 

n  KD 

M 

iH  rH  cn  iD 

Cm 

Cm 

- 

It] 


es 

< 

(- 


r<J 


D-i 

O 


in 


<-1 

o 


DJ 
r  1 

Ij^ 

a: 

LD 

CO 


Q  is 

^  I — I 

a:  < 
a: 
oi  o: 
DJ  CJ 
E-i 


iJ 


c 


Q 
CO 

o  u 


o 

a 


CO 

Id 


.o 

O 


CO 

H 

r  1 
CO 
Cu 
Cu  lO 
CO  CO 


H  is 
H  o 


cj  H 

E-<  CO 


D  D 


CO 

ro 

o 

p* 

W.0 

'X. 

o 

CD 

PI 

1—4  wD 

CO 

m 

in 

CO 

(Ti 

ro  wO 

w 

o 

ro 

CO 

CM 

CD 

CM  CO  rH 

in  WD 

cr: 

C'l 

P3 

r-4 

N  WD 

p-  ro  P^ 

rH  CO 

CL 

,CM 

>1 

«  •  • 

«  • 

Eh  wo 

p^  in 

WD  O 

• 

1 

CM  CM 

CM 

in 

CD 

in 

CO 

rH  ^ 

<c 

^  CD 

ro 

ro 

wD  o 

H 

in  • 

• 

• 

• 

•  • 

CO 

U  r-< 

rH 

1— 1 

i-H  (M 

CO 

P3  CM  WD 

in  ^ 

QJ 

N 

in  00  o 

O  rH 

in 

X 

•  •  • 

•  • 

Eh  pi 

WD  ro  ^ 

in  rH 

M 

u  ro 

D 

r-l 

wO 

<T>  cr> 

rH 

CM  iH  1 

CO 

Su  CD 

r-f 

pH 

P) 

^  rH 

^  1 

1  1 

x 

rtl  • 

• 

* 

« 

•  • 

in 

W  o 

o 

o 

o 

o  o 

O 

1 

! 

>  -^  in 'd*  in 

CO  in 

psi)  p-  wD  WD  P- 

Ch  CO 

X  o  • 

»  •  • 

•  • 

Eh  0  wO 

CM  CM  P^ 

in 

in  ^ 

rH 

o 

r-  CO 

o  in  in  in  in  i 

—i  cr* 

in 

ro 

cr» 

CD  o 

C  rH 

I  1  1 

1 

o 

fTj  • 

• 

• 

• 

«  • 

•r^  1 

Q 

o 

cr« 

pj  in 

•D 

wD 

WD  in 

in  WD 

ro  o  ^ 

CD  ^ 

bJ  X)  -«d* 

rH  WD  CM  WD  WD 

N  "iH  • 

•  •  • 

•  • 

Eh  d  (M 

cr»  CD  CM 

CM  CD 

> 

CM  ^  rH 

WD  in 

a 

P- 

WD  in  CM 

CO 

r-H 

o 

ro 

ro 

O  CO 

H 

Q 


o 

r>* 


00  ^.0  00  cr> 

<J0  '>0  lO 


:> 

ij 

CO 


in 

Vj 

0 

Jj 

(D 


IT)  cr>  CN 

r— I  rH  ^  ^  iH 


rj 

I 


<r>  o  in  ^  o 

^  rH  n 

n  m  CN  CM  t 
I  I  I  I 


CO 


o 


o  n  n  o  CT' 
00  <.0  ^  CD  cn 
lO  in  “d* 


wo 


CO  iH  rH  P^  rH 
CM  <7>  Ti’  CD  CM  in 


X 


Q 

U 

CO 

M 

Cu 


n  cr>  (M  in  o  o 
r-<  r-i  ro  wo 


Q 

U 

CO 


1-4  cr»  CO  wo  (D  00 

p*  WD  CM 

n  (M  1-4  r-(  I  \ 

till 


Lu 


n  o>  CM  in  o  o 

f-4  r-H  m  WO 


sa'fSET  PEAK  EFFECTS  USING  ONLY  INNER  RADIUS  OF  12' 
AND  STATION  ELEVATION  5000  FT  ABOVE  MEAN  SEA  LEVEL 


CO 

cn 

o 

CD 

s 

o 

CD 

CM 

rH 

CD 

CO 

ro 

in 

CO 

CTt 

ro 

CD 

in 

in 

w 

o 

'Cf* 

n 

00 

CM 

o 

<7l  rH 

Oi 

CM 

CM 

rH 

N 

cn 

in  in  rH 

o 

04 

CM 

• 

• 

# 

• 

• 

• 

. 

EH 

7 

7 

1 

1 

1 

1 

1 

1 

in 

CD 

CD 

CD 

CD 

o 

o 

O 

o 

O 

o 

rH 

• , 

Eh 

in 

• 

• 

• 

• 

• 

• 

cr» 

S 

u 

rH 

r-\ 

rH 

rH 

CO  in  o  CM  CM 

n) 

1 

\ 

1 

« 

1 

I 

N 

CD 

cr»  rH  CM  CD 

p 

in 

X 

• 

• 

• 

• 

cd’ 

• 

H 

in 

in 

CD 

CD 

M 

u 

•d* 

n 

CM 

CM 

CD 

1 

1 

1 

1 

1 

1 

CO 

CM 

CM 

(M 

CM 

rH 

rH 

Oj 

• 

• 

• 

• 

• 

• 

in 

rH 

rH 

rH 

rH 

0 

CD 

1 

1 

1 

1 

1 

1 

>  00 

in 

CM 

CO 

rH 

«d« 

CO 

0^ 

in 

X  0 

• 

• 

# 

• 

• 

• 

Eh 

fO 

ro 

in 

CD 

rH 

O 

CD 

CD 

CD  00 

<T» 

o 

d 

O 

fc! 

• 

• 

• 

• 

« 

« 

•H 

Q 

tno 

o 

o 

o 

O 

rH 

CD 

CM 

pi  in 

in 

in 

in 

in 

in 

in  CO 

rH 

N  4J  00 

cr» 

o 

O 

CD 

CD 

N  ‘rH 

• 

• 

• 

• 

• 

• 

EH  C  ^ 

in 

in 

in 

> 

CM 

CM 

CM 

CM 

(M 

a 

w 

g 


u 

Cd 


cn 


Q 

u 

Cd 

M 

Cu 


n  in 

CD  CO  00 
^ 


CD  Lll  CD 
CD  CD  CD 
-d* 


in 

kl 

tu 

JJ 

(U 

E 


>>1 

H 


^  ^ 

CM  (Nl  CM 

in  in  m 


^  ^  ^ 
(N  (N  CM 

in  in  in 


W 

EH 


cn  0^  CM  in  o  o 
rH  f-l  n  VD 


Q 

iJ 

Cd 

M 

Cm 


c^  in  r-j  o  in  o 
cr>  o  iH  r^  o  n 

m  ^  ^  ^  ^ 

rH  iH  rH  I— l  rH  rH 

I  I  I  i  I  I 


CD  IM  m  CD  CO  CM 

CD  <J'  cn  cr>  ro 

•  ••••• 
O  O  O  O  O  iH 

r-i  r-\  r-i 

I  I  I  I  I  I 


n  cr>  CM  in  o  o 
fH  iH  n  'X 


DJ 
U 
< 
Q  Cu 

a  O' 
CU  ID 

cn 
Cc  ^ 

ai  M 
CU  < 
s  ct: 
s  q: 
M  CU 


Ul 

rn 

CTt 

'ST' 

<->4 

kO 

CO 

rH 

0 

cn 

M 

in 

in 

CX^ 

in 

cr> 

Oi 

in 

f-M 

<j  ^ 

CO 

cr> 

0”^  in  ••  • 

CU  u  o  o  o 

n)  t  I  I 
in 


KD  CO  O  n  KD 
vjo  04  m  o 
cr»  o  rH  m  in 
n  rn  in  o  CO 
04  in  CD  (M  in 

r-<  rH  fH  CS  (N 


in  CD  O  CN 
CT>  O  O  O 


O  O  rH  rH  r~l 

>  I  I  I  i 


U'U'inintNcD^O'ce* 
O'!  KjOr-<cvjmro^'^in 

I  >  I  t  I  I  I  I 


H 

c 


H 
u 

CU 
Cj:. 

(I4  uJ 
CU  CU 

IS 
<  o 
CU  M 

E-<  H 
CU  CO 
CO 

^  Q 
D)  S 
-yj  < 


m 

tninriinkD'd*r-ip>.cN 

0 

r~lOOOr— f<  LO 

0 

(Tf . 

U  Q 

Q 

tpnJfimrnmmrnm 

2  CU 

gininintnLninLnin 

0  H 

< 

0  U 

is  0 

> 

M  1—3 

u 

CO 

CU 

rn  uo  D  in  'D  kO  kO 

D  S 

CDCDCDCDCDCDCDCD 

0 

CO  M 

Sh 

r.n 

Q 


> 

uJ 

CU 


ki 

(U 
iJ 
(U 
6 

w  rn 
CO 
is 

CO 


CO 

D 

M 

Q 

< 


in  o-  CD  in  kO  'JD  kO 

CO  (D  CD  CD  CD  CD  CD 
^  ^  ^  ^ 


kDr^cDa>Or-<cvim 


cDcDCDcno^cTtcricr) 

^  rH  pH  rH  f-H  i-H  rH 


'^'-4rHrHrHrHrHrH 
I  I  I  I  I  I  I  I 


r^r^r^r^o»rN.f>Hr^ 

N  rnrorororornroro 

H  r^(NJr4(N(NICM(NCN 
'  I  I  I  I  I  I  f 


in 


t>i 


o 

>r^r^a>dkD^rHD 

•ui^OrHCNro^inin 

. . 


ly  O  rH  rH  rH  rH  rH 

rn  ro  ro  ro  ro  ro 

C 

•H 


n  n 


inorHrHfnrnt^'d^in 

N4Jinr^r^cno^^,;*«t«^ 

N  ‘tH  •••••*»» 

H  Cincomo^kDr^f^o 

D  rH  f-H 


C^lOrHOlTlLnrHn 

>1  rHD^kDrHOrHp*) 


H  inrHcriO'kDin'd'm 
CNCNrHrHrHrHrHrH 
I  I  I  I  I  I  I  I 


CD(7>0^r^kD03r^D 

rnokDvor'JinkOin 


ominr^coorHCH 

I  r~4  r-l  F— 1 


kDo»cDcr>OrH(NJn 


UI 

e3 

< 

H 


CDiH^ini^r^CDCD 
N  LnvDvjDvD^DvDvjDvD 


H  UD(^<J3vjDvJ3vO*s0vD 
I  I  I  I  I  I  i  I 


oDinr^rHocrtcor^ 

N  in  in  in  in  Lf)  ^  ^  nj* 

^  ••••••»« 

Eh  r'3cv)(MCNrocNfs4c>J 

'^rHrHiHrHrHrHrHrH 
I  I  I  I  I  I  I  I 

W 

o 

jnx)  cr>otMminuDn-Q3 

X  . . 

Eh  (Uinv^kjDUDvDvJDi^kD 

oooooooo 

C  rHrHrHrHiHrHrHrH 
•HI  II  I  I  I  II 

WrHOnrHnJGDUDUS 
N4J  ncnnJi^r^rj^CD 

N  . . . 

D  «>Dinin«d^'d«^^^i« 
p^r>*r^r^ps.r^t^r^ 


^aio^^onr>-uD 

CDinrHTf^rHOrH^ 


iHa3vD^rf*nnJrHO 

cor^r^r^r^f^r^r>* 

mnnnnnron 

I  I  I  I  I  I  I  I 

r^rH^ooniLnaio 

rtnorHvonna^ 


nJcr»r^innc^rHO 

cDr^rvf^r^r^rHpN. 

mnnmnnmn 

I  I  I  1  I  I  I  I 


<^p^a3a>OrHr<irn 

rH  rH  rH  rH 


o 


ro 


CO 

ro  (N 

Cu 

s 

m  CD 

'  O 

CO 

dS  o 

(N 

M 

CM  o 

rH  CO 

X 

in  o 

D 

Co  I — 1 

Ob 

m  m 

O  Q 

H  rH 

C 

c 

^  CO  CD 

Q  iX 

Eh 

in  •  • 

uU 

Cd 

u  o  o 

Cd  c 

0  I  1 

M 

Cu  X 

0 

EH 

hb 

u 

X  M 

CO 

Su  o  o 

Cd  s: 

bo; 

fij  •  • 

V-/  CM  CM 

"z  o 

M  uJ 

1  1 

Cd 

■ — ' 

=  VH 

in  CD  m 

ro  Cu 

i-i  cr^  <T> 

O 

Cj  •  • 

<  a: 

Q 

m 

Cd 

e  in  in 

O  H 


rH  •d'  ro  ro  r-< 

(Tt  <T>  n  CD  cr>  cn 

r>v  CTt  o 

rv  «jD  OJ  o  U) 

O  in  ID  O  O  LD 

m  ro  CD  CD  n  m 

rH  CTl  n  CD  ^ 

CD  CD  OD  00  CO 


o  o  o  o  o  o 
I  I  I  I  I  I 


r;f*  ni  oo  r'J  cn 

o  O  CD  CD  CD 


(NJ  (N  iH  rH  I— <  rH 

I  I  t  I  I  ) 


rn  ^  r-<  rH  r-< 

cr>  CT*  CC  CD  iD  vD 

.  •  •  .  •  « 

(D  m  IN  c^  IN 

in  in  in  in  in  in 


s 

M  O 
CO 

D  Cd 
U 
CO  CP 
H  < 
U 

Cd  cr; 

Cxj  <C 
Cu  > 
Cd 


< 


C  Q 
Cd  S 


uJ 

Cd  m  n  CD 

CD  CD  CD 

^  ^  ^  "'o' 

in 

Q 

(U 

OJ 


n  CD  n  m  fo 
CD  CD  CD  CD  CD 
^  ^ 


>  0 
tJ  a 
Ci3 

Eh 

CO 


CDCOCDCOC^C^rON 

T^'-rf'-^^mLninin 

rH  r— <  rH  rH 


fvjcNr^r^oooo 

N  ......•• 

»>,  rHrHrHiHCOOOCDCO 
H  I  I  i  t  I  t  I  I 


inrHCDCDinininin 

CMnCNC'JCDCDCDcD 

N  . . 

X  rvjfsjtMromiDLnin 

H  ^  I  I  I  «  ’  ‘  • 

w 

o 

>cDcDCD^^^f^'~* 

x)r^r*^r^r^r^r^r^r*^ 


Jx* 

H  ro  n  n  ro 
C 
*H 

tn^c^oQi—iiHor^o 

oj^r^cDt^^^c^c^ 

N  ‘H  •••••••  • 

N  CID^C^lQ^0^C^C^C^D4 
^  D  I  I  I  I  rH  rH  rH  rH 


r^^rO(Hcricr>D>co 

fHrHrHrHOOOO 

>,  . . 

CDODCDCDCDCDCDCD 
^  I  I  I  I  I  I  I  I 


coroocMCNOCDro 
iniTiCDCDnmc^JC^ 
M  . 

iH  rH  <H  rH  I  I  I  1 


9-?  CO 

Ld  r5  Q 

CO  M  iJ 

s  q  Cd 

D  ^  M 

CO  Di 


cMOiHmcJOrHro 

H  ro  rH  n 


Q 

u 

DJ 

M 

Cl. 


CMOrHrOtNOrHrO 

rH  n  rH  n 


» 


00 

hJ 

«J 

CO 

< 

H 


o 

cn 

'  O 

(N 

rH  CO 

D 
(1*  M 


CO 

S 

CO 

M 

<X 

Oi 


rOCJrH^rOC^i-Hrd^ 

rncDcncr>rnoocr>cr> 

cr>or^r^o^or^r^ 

cMor^ijDr^or^vr) 

inooininooLn 

nnnnrnromn 

ininintjDinLnLnuD 


Q  oc: 

H 

. . .  . 

tj 

Id 

Ur^CMCMCMCMCMCMCM 

U  < 

(U 

M 

tn 

Cm  X 

EH 

M 

(jocriC0‘>fr^f^kDrM 

cd  M 

CO 

Sj'ct^nnncMCMCMCM 

CJ  S 

1^ 

X 

'-'rHrHrHrHrHrHiHrH 

X  Q 

1  1  1  1  1  1  1  1 

M  vj 

Cd 

=  M 

{nLDLnLOLnkDkDkDr^ 

CO  Cm 

rH'^kDkDkOnnmn 

<c  o:: 

o 

Q 

vrcMcMCNCNfnconn 

Cd 

gr^r^t^r^LOLninin 

O  Eh. 

S  D 

> 

M  O 

id 

CO 

Cd 

rHrHrHfHrHrHrHrH 

D  Cd 

OOOOOOOO 

uj 

X 

^r^r^r^rsr^rvt^t^ 

CO  CQ 

EH 

tn 

H  < 

Q 

Vj 

CJ  M 

0 

Cd  a 

4J 

Cm  < 

> 

d) 

Cm  > 

lJ 

B 

Cd 

Cd  • 

x-^rHrHrHrH^T^^Tii'^ 

< 

OOOOCMCMCMCM 

Eh 

t-^r^i^r>inLnincn 

H  Q 

EH 

rH  rH  rH  rH 

O  S 

CO 

tj  < 

CO 

CO 

- 

E-i  Id 

Q 

Im  m 

tJ 

^  9 

Cd 

cMOrHncMOrHrn 

O 

M 

rH  n  rH  n 

X  iX 

Cm 

ncn^nocno^cn 

vdudvdvdvoCdvovd 

N  ••»•.... 

VjDvOisOvfi^*ci«Tj<Tt' 
H  I  I  I  I  I  I  i  I 


N 


tn 


O^rHl^lDfOrHCMOJ 

f-4fnr^r^r-t.rHrMi-4 


c^r^oJCNo^o^o^cn 


till 


-.0 

>cN3«jDnncriCDrvr^ 

»•••*••• 

X  oj^^uDvot^nronn 
H  o  o  o  o  I  I  I  I 

C  rH  iH  ^  rH 
•Hill  I 


tninocorocTtcnaDvD 

•i->cr>cN30ocr»cncncn 

N  . . 

N  ccnoooc^c^c^c^J 
Eh  Dr^Jroronr^r^i^h^ 

fs* 


CDl^P^VDC^CDCOr^ 

c^r^r^r^oooo 

^  »•••••»• 

>.  LninLnin»jD«^vDvD 

Eh 

ro  n  n  n  I  i  i  i 
I  I  I  I 


t^^CNJf^rHrHrHCO 

rHrj«r)CNa>cr>cr»cD 

X  . 

X  ^^^•;J'kDvDkOVD 
H  '^kOkD'JDCNCMCMCS 
cn  fo  n  n  I  i  i  i 
I  I  I  I 


Q 

uU  sr'.«.sr%^ 

W  CM  O  rH  ro  CM  O  rH  cn 

HH  rH  n  rH  ro 

t* 


U^fSET  PEAK  EFFECTS  USING  A  3"  INNER  FIELD  WITH  A  12' 


CM  CM  CD  CD 

N  CM  CM  O  O 

>t  .  .  .  . 

CO 

CM  c\>  C'J  r<J 

H  rH  rH  CD  CD 

CO 

r^  CD  iH  CD 

1  I  I  1 

D 

CO 

Ti*  n  -"d*  rn 

M 

M 

cn  o  cri  o 

'  _  > 

q 

cx 

in  cn  LH  cr> 

CL 

ro  ro  n  fo 

ID  in  O  CD 

(X 

N  CM  CM  ID  in 

CD  O  kD 

X  .  .  .  . 

CJ 

< 

^  in  ^  ^ 

H  CM  CM  in  in 

u 

E-< 

tn  •  •  •  • 

1  1  1  1 

CQ 

DJ 

C 

1 

2 

1 

2 

C 

(U 

in 

M 

to 

.  o 

>  ^  CM  CM  O 

M 

u  o  cr>  in 

>iaj  in  in  in  in 

> 

CO 

sj  in  CM  ^ 

X  o  •  *  •  • 

trj  •  •  •  • 

^  QJ  rH  rH  ^  ^ 

c 

w  rn  n  ro  ro 

ro  cn 

till 

c 

X 

•H 

H 

M 

tn  uD  CM  rH  m 

m  CD  CTI  CD  ID 

s 

r— 1  O  rH 

N  OJ  rH  ID  OD  n 

O 

•  •  •  • 

N  ‘H  »  •  •  • 

Q 

Q 

in  ^ 

Eh  c  id  id  in  in 

U 

g  in  in  in  in 

D  rH  rH 

U 

1  t 

M 

> 

Du 

ul 

W 

n  m  rn 

^  in  IH  CM 

(X 

CD  CD  CD  CD 

t>,  O  CD  O  CD 

W 

^ 

>1  .... 

H 

tn 

H  in  in 

tD 

Q 

1  1  1  1 

O 

OJ 

OJ 

> 

(U 

U 

g 

CM  in  ^ 

U 

w  ro  n 

CM  in  CD  in 

C 

CD  CD  CM  CM 

M  .... 

Is 

in  in 

H  iH  rH  O  O 

Q 

H 

r— 1  rH 

CM  CM  1  1 

2 

CO 

CO 

CO 

CO 

:d 

:d 

^  ^  ^ 

JD  %.  v  ^  » 

M 

o  o  o  o 

M  O  O  O  O 

q 

Q 

CM  in  CM  in 

q  CM  in  CM  in 

1— 1  rH  (— <  1— 1 

^  rH  pH  rH  rH 

(X 

cc 

<x 

(ij 

>J 

ca 


C^3 

rH 

CO 

IN  CM 

kD  n  n 
N  vD 

>s,  •  •  •  • 

H  VD  IJD  ^  ^ 

CO 

X 

rH  CD  rH  CD 

till 

<  to 

CO 

fo  m 

M 

)H 

cn  o  cn  o 

XQ 

cc; 

in  cn  in  c^ 

H  ^ 

Dj 

n  ro  n  m 

cn  CM  <JD  CD 

M  a 

s 

o  o 

N  CM  CM  IH  H 

X  •  *  •  • 

M 

< 

CM  CM  i-m  r-l 

9-<  CM  CM  cn  cn 

D  uJ 

H 

tn  •  •  •  • 

rH  rH  iH  iH 

j  m 

CJ 

c 

4 

5 

4 

5 

I  1 

U  <C 

OJ 

in 

M  M 

in 

•0 

Cm  K 

M 

u  n  o  o  CD 

>  cr»  rH  n  in 
j>,x)  0^  o  0^  cr> 

a:  > 

CO 

sj  cn  vD 

W  0  •  •  •  • 

CJ 

(^  •  t  •  • 

Em  0  ID  n  n 

X  <3 

>-•  CM  CM  CM  CM 

O  O  1  1 

'z 

till 

C  iH  rH 

M  X 

H 

•H  I  1 

SM  inuDrHC^WD  tnOrJrHn 

cnS  iHincar^o  NJJoincri^ 


Q 
U 
O  U 


M 

cn 


M 

Cq 

Di 

U 


O 

Q 


U 


• 

tpn  rj  ur» 
e  in  in 


rH  rH  iH 

O  O  O  O 


N 

H 


>1 


C  ^  n  VD  ID 
P  CSJ  r3  VD  ID 
rv 


CO  O  (T>  fsJ 
CD  H  O 


CO  Eh 

H 

in 

EH 

CM  CM  n  n 

Eh  X 

P 

Su 

ID  ID  ^  ^ 

U  O 

(U 

n  n  1  1 

X 

OJ 

1  1 

Cm  ' 

> 

OJ 

Cm 

P 

6 

'  CM  CM  (N  t-l 

Cl] 

W 

>-  i-l  t-l  ^ 

« 

rH  ra  rv  ^ 

< 

O  O  CM  CM 

X 

•  •  •  • 

EH 

X 

in  in 

Eh 

iH  o  rn  ro 

EH  P 

Eh 

i-H  iH 

ID  VD  CM  CM 

SS 

CO 

n  m  1  1 

1  1 

CO 

CO 

CO 

CO 

Eh  P 

P 

%  %  ^  ^ 

P 

%  %  %  ^ 

S  M 

o  o  o  o 

M 

o  o  o  o 

P  P 

P 

CM  in  CM  in 

Q 

CM  in  CM  in 

O  < 

rH  1— 1  rH  rH 

r-l  rH  rH  rH 

X  P 

cc; 

a 

SOMSET  PEAK  EFFECTS  USING  A  1'  OUTER  FIELD  WITH  A  150' 
RADIUS  AMD  3"  INNER  FIELD  WITH  A  VARIABLE  RADIUS, 

THE  STATION  ELEVATION  IS  ON  THE  TERRAIN  SURFACE 


fo  n  fo  m  (N  (^4 


N 

CM  CM  CM 

CM  CM  CM  CM 

4 

4 

4 

4  4 

4 

4 

CO 

•ch 

in 

in 

o 

ID 

r) 

cr> 

E-< 

rH 

rH 

rH 

1 

1 

1 

1 

ro 

pH 

CD 

in 

cr» 

o 

1 

1 

1 

1  1 

1 

1 

to 

m 

CD 

CD 

n 

in 

cn 

M 

o 

CD 

CD 

m 

in 

vD 

c^ 

a; 

a> 

cn 

o 

CM 

ID 

c^ 

a. 

fH 

rH 

rH 

rH 

CM 

p-  in  in 

in  in  in  in 

N 

CM  CM  CM 

CM  CM  CM 

CM 

in 

in 

CD 

cn 

o 

4 

4 

4 

•  4 

4 

4 

< 

in 

H 

CM 

CM  CM 

2 

2 

2 

CM 

in 

• 

• 

• 

• 

4 

• 

• 

1 

1 

1 

1  1 

1 

1 

Cd 

u 

2 

2 

2 

2 

CM  CM 

CM 

0 

tn 

in 

0 

> 

0^  CD  CM 

CM  n 

rH 

M 

u 

cr> 

r-l  CM  m 

in  ID 

o 

>iX) 

m  CM  CM 

n  ^ 'd'  in 

CO 

Su 

in  VD  UD  iJD 

ID  VD 

p- 

x  o 

4 

4 

4 

4  4 

4 

4 

ni 

• 

• 

4 

4 

• 

• 

4 

fH  0 

fH 

fH 

rH 

rH  rH 

pH 

rH 

fO 

ro 

n  m 

m  m 

n 

m 

m 

ro 

ro  m 

D 

fO 

1 

1 

1 

1 

1 

1 

1 

c 

•H 

in 

ro 

CM 

CM  CM 

CM 

CM 

CM 

in 

cr> 

ID 

rH  in 

rH 

pH 

o  o  o  o 

o 

o 

o 

b4  jj 

CJ^ 

CD 

o 

fri 

$ 

• 

« 

• 

4 

4 

• 

M  ‘H 

4 

4 

4 

•  4 

4 

4 

Q 

b^LO 

in 

in 

in 

in 

in 

in 

H  C 

ID 

ID 

VO 

VO  VD 

VO 

vD 

6  ID 

in 

in 

in 

ITi 

in 

in 

pH 

rH 

rH 

rH  rH 

rH 

rH 

1 

1 

1 

1  1 

1 

1 

> 

u 

VD 

Cd 

ro 

m 

n 

m 

n 

m 

n 

in 

cr> 

O 

^  n 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

>1 

rH 

n 

n 

CM  O 

cr> 

CD 

S 

4 

4 

4 

4  4 

4 

4 

in 

EH 

in 

in  in 

5 

5 

4 

S 

1 

1 

1 

1  1 

1 

1 

0 

jj 

Z> 

iJ 

0 

B 

ID  VO  m  C7^ 

Cd 

m 

m 

n 

n 

m 

rn 

m 

CD 

m  rH 

a»  VD 

in 

CJ 

CD 

CD 

CD 

CD 

CD 

CD 

4 

4 

4 

4  4 

4 

4 

tJ* 

H 

rH  CM  CM 

rH  iH 

rH 

rH 

CM  CM  CM 

CM  CM  CM 

CM 

CO 

CO 

CO 

D 

D 

M 

M 

«4 

% 

Q 

1 

2 

3 

4 

in  vD  o 

Q 

rH  CM  CO  ^  in  ID  o 

< 

rH 

3  • 

rH 

cc: 

(X 

MOUNT  SCOTT  EFFECTS  USING  A  1'  OUTER  FIELD  WITH  A  150' 
RADIUS  AND  3"  INNER  FIELD  WITH  A  VARIABLE  RADIUS, 

THE  STATION  ELEVATION  IS  ON  THE  TERRAIN  SURFACE 


I 


csi  03  in  in  in  ID  ijD 


>1  . 


CO 

in 

in 

O 

VO 

n 

cn 

H 

VO 

VO 

VO 

VO 

VO 

VO 

VO 

s 

m 

rH 

00 

in 

cn 

o 

1 

1 

1 

1 

1 

1 

1 

CO 

m 

rH 

00 

00 

n 

in 

00 

M 

o 

03 

cn 

m 

in 

VO 

<T\ 

a 

cn 

cn 

o 

CM 

VO 

on 

CLi 

iH 

rH 

rH 

rH 

CM 

cn 

CM 

CM 

CM 

CM 

CM 

N 

rn 

iH 

CM 

CM 

CM 

CM 

CM 

00 

CM 

o 

rH 

iH 

O 

o 

X 

• 

4 

4 

4 

4 

4 

4 

< 

rH 

CM 

CM 

CM 

CM 

CM 

CM 

EH 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

Eh 

w 

• 

• 

• 

• 

• 

4 

4 

iH 

iH 

rH 

rH 

rH 

•H 

iH 

W 

u 

in 

in 

in 

in 

in 

in 

in 

l 

1 

1 

1 

1 

I 

1 

a 

w 

tn 

0 

> 

CM 

VO 

rH 

VO 

cn 

cn 

M 

u 

cn 

cn 

tD 

VO 

N 

o 

>1 

jj 

VO 

VO 

03 

00 

00 

cn 

CO 

u 

00 

CX) 

00 

00 

00 

00 

cn 

X 

0 

• 

4 

4 

4 

4 

4 

4 

fd 

• 

• 

• 

• 

• 

• 

• 

Eh 

(U 

in 

VO 

VO 

vO 

VO 

VO 

VO 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

o 

O 

o 

o 

O 

o 

o 

1 

1 

1 

1 

1 

1 

1 

c 

rH 

iH 

rH 

iH 

rH 

rH 

rH 

•H 

1 

1 

1 

1 

1 

1 

1 

WLninmnf^CMrH 

iHCD03  03COa3a3a3 

O  td . 

Q  tj'n  n  n  n  n  n  cn 


EH 

CO 


rH  rH  iH  rH  rH  rH  rH 

o  o  o  o  o  o  o 

r>* 


o  o  o  o  o  o  o 


win^mncs^n 

Nojnor^on^in 

••••••• 

B  CNCMcsncnncn 

DC>ICSCSCSJCNJCMCS 


cs  CO  vD  00  ro  o 

t>i  vDOO^vDi^r^ 

>1  . 

H  oc^c^r^r5c^cs 
VO  VD  U3  U3  VD  vjO  (D 
n  n  r)  n  n  n  n 

I  I  I  I  t  I  I 

n  00  ^  o  CM 

«  r^o^iDiniDi^oD 


H  rHcnooooo 
vo  in  vo  ^  ^  U3  VD 
n  n  m  n  n  n  n 
I  I  I  I  I  I  I. 


TABLE  L3 


Q 

•Sh 

< 

'  Q 
U  U 
W  iJ 
M  CCl 
Ct4  <! 

M  o 

2  M 

ro  CO 
cq 

< 

M  . 
O  M 
S;  Q 
M 

2?  a 

X 
CQ  H 
Eh  O 
O  03 

CO  X 

Eh 

M 

s 

X 

a.  CO 


CO 

CO 


Q 

uj 


CO 

X 

CO 

M 

tX 

Oi 


VJD  03  LD 
^  ^  ^  ^ 
LO  03  IJO  CD 
O  rH  O  rH 
03  (N  03  CH 
f-H  .— 4 


ID  CD  CM 

^  in  ^  ^ 

Eh  03  •  •  •  « 

CO  u  CM  m  CM  ro 
(U 
to 


M 

to 

iX 


u 

Su 

ni 


VD  rH  CM 

Lo  n  ir> 
•  •  •  • 
ro  CD  n  m 

)  >  t  I 


03 


g 


po  <H  in  ch 
rH  o  n  ^ 

trm  in  ^  ^ 
H  in  in  in  in 


> 

►J 

W 

(U 

JU 

(U 

s 


CO  n  n  fo 
CO  CD  CD  CD 
ri*  ^ 


Eh  CO 
CO 

03  0^ 
X  CO 
X  fH 
CO  X 

o 

u 

X  fO 
M 

Cu  i< 


X 

u3 

w 


N 


X 

Eh 

CO 


cn  n  ^  ^ 

X  CD  CM  CM 

^  in  in 

rH  rH 


N 


ro  CM  03 

CM  CM*C33  o 
•  •  •  « 

rH  rH  Cv  CD 
I  !  I  I 


cn  CD  o 
CM  CM  in  CO 
•  •  •  * 
CM  CM  in  in 
>111 


tn 

o 

.  ^  rH  rH  03  03 

C^  O  CM 

WO*  •  «  • 

Eh  (U  rH  rH  ^  ^ 

^  n  n 

C 


W  VO  in  VD  CD 

N  •*«>  CO  CM  VO 

N  , 

Eh  c  VO  in 

ZJ  pH  rH 
I  ( 


>1 

fcs 

Eh 


w 

Eh 


fO  o  o 
ro  ro 
•  •  •  • 

in  ri*  m  *qj« 

>  »  I  I 


o  ^  CM 
CM  rH  rHrH 

•  •  •  • 

CM  CM  O  O 

CM  CM 


M 

^  ^ 

o  o  o  o 

w 

% 

o 

M 

in  CD  in  CO 

M 

in 

CJ 

< 

rH  rH  rH  rH 

*  ^  V  V 

§ 

rH 

a 

^  ^  ^ 
in  VO  in  VO 

a 

in  VO  in  so 


TABLE  14 


O 

W 
W  iJ 
M  CQ 
Cl,  < 
M 

S  < 

S  > 

M  O 

=  M 
CO  CO 


M 

o 

S  Q 

M  Kti 
CO 

^  X 
CO  H 
H  O 
a  cQ 
CO 

Ct,  X 
Cl,  H 
CO  M 


H  a 

O  lJ 

u  w 

CO  M 
Cl, 

E-<  _ 
s:  iX 
X  CO 
O  H 
X  3 
O 


Cc,  c 


ld  cn  lD 

^  ^  ^ 

UD  CD  VD  W 
O  pH  O  rM 
(VJ  O 


in  cr>  in 

<(J  ^  ^  rH  O  O 
£-<  in  •  •  •  • 
CO  u  in  LO  in  vD 


M  U  r-i  (Tt  ^ 

CO  su  CO  cn  X 

jiO  fd  •  •  *  • 

w  fsj  CS  CM 
III! 


in  m  o  pH  cr> 

pH  CD  pH  pH  n 

O  td  •  *  •  • 

Q  lT»ro  ^  in  in 

fi  in  in 


pH  pH  pH  pH 
O  O  O  O 


CO  ^  pH  pH  ^ 

O  O  CM 

s  in  in 


^  ^  ^ 
o  o  o  o 
in  CD  in  CD 

pH  pH  pH  pH 
«  ^  ^  ^ 
%  »  ^ 
in  ID  ld 


in  ir>  ’o>  cy 

L0  LO  V0  VD 

•  •  •  * 

UD  vO  ^  ^ 
I  I  I  I 


CM  ^  O* 
N  CM  pH  pH  CM 
{xj  ,  ,  .  . 

^  CM  CM  CTl  0> 
^  ^  rH  pH  pH 
I  I 
in 
o 

>  VD  0^  pH  cn. 
>1X9  ID  CD  LO  ID 

i*j  o  •  *  •  • 

^  (U  CD  VO  n  n 

o  o  I  I 

C  pH  pH 
•H  I  I 


in  ^  CD  CD  O 

bj  jj  pH  cr»  o  cn 

N  ‘pI  •  •  *  • 

Eh  c  n  cm  CD  in 
3  CM  CM  CD  CD 


O  CD  CM  CD 
CD  in  pH  CD 

t  •  •  •. 

CM  CM  m  CM 
CD  CD  ^  ^ 

men  \  I 


^  O  !>•  CM 

n  ^  0^  o 

•  •  •  • 

o  o  CM  m 

CD  CO  CM  CM 

n  n  I  I 

t  I 


o  o  o  o 

in  CD  in  CD 

pH  pH  pH  pH 
».  »  k  ^ 

«  «  %  . 
in  co-invo 


ANALYSIS  OF  RESULTS 


uJ 

0“; 

CU 

cr. 


a 


Q 

Co 

o 

00 

r-< 

cc: 

U 

Cl] 

W 

Q 

OR 

o 

V 

Cl4 

U 

o 

Q 

CU 

Ln 

n 

ii‘? 

cu 

H 

Du 

CO 

CO 

Q 

Q 

CU 

E-i 

m 

ra 

Q 

Q 

<Li 

Q 

u 

uJ 

M 

x 

E) 

o 

S 

r> 

o 

Oi 

< 

x 

X 

a 

cr; 

CO 

m 

m 

s 

CO 

X 

CO 

r 

H 

< 

E~f 

lx; 

•*  'I 

Cl] 

CJ 

CO 

CO 

CQ 

lJ 

u 

I>« 

< 

M 

M 

(n 

Q 

Q 

H? 

:e 

< 

C 

^Li 

a: 

(X 

Q 

Q 

Q 

Q 

^-l 

CJ 

M 

CL 

IH 

iX 

Cti 

r*i 

Cu 

o 

o 

D.' 

a 

Q1. 

(X 

[J 

DJ 

CO 

CO 

H 

R 

D 

E) 

■— 1 

O 

M 

O 

:j 

DJ 

ClI 

CO 

E 

X 

X 

X 

-< 

H 

R 

-fc- 

(X 

CQ 
CO  CO 

o 

cn 

CO 

U  D 

M 

< 

Q 

X  i-J 

R 

CX 

CO  < 

X 

CO 

O 

o  > 

00 

o 

I 

cx 

»-< 

cx  Q 

Q 

CO  R 

(X 

CO  CO 

X  CQ 

>  U 

R 

O  R 

X  CO 

X 

X  <rl 

U  M 

O  Du 

o 

O  CO 

CO  X 

U  Du 

CJ 

CJ  ct: 

CQ  M 

CO 

CJ 

lO 

Ct; 

ct:  CQ 

R  X 

D 

O  X 

O  [H 

O  M 

X  Q 

Du  HI 

CJ 

X 

<  CO 

< 

<  CO 

CQ 

a: 

o  cx 

S!  Du 

U  M 

CO  M 

CO  ct: 

Du 

Uj 

CJ  CQ 

Q  CO 

CO  CO 

M  R 

X  CO 

CO  R 

X 

X 

Q 

CO 

R  X 

CO 

CQ 

X  CJ 

M 

q:  X 

M 

CQ  < 

R  ^ 

D  M 

CQ  X 

X  w 

R  X 

M  d, 

R  ct: 

CO 

CQ  R 

< 

X  CO 

X  ij 

M  M 

o  cx 

P3  R 

o  < 

QS 

CO  o 

M 

Du  M 

R  O 

ct:  u 

><  CO 

Q  CU 

R  w 

o  < 

R  O 

o 

X  X 

o  ct: 

M  Ct: 

=  R 

O  Oi 

>  l>  3 

ro 

o  < 

X  Q 

<  > 

R 

ct: 

O  M 

cx  X 

CO  U 

Q  o 

M  a 

o  o 

R  CO 

u  o 

R  o 

u 

<  CX 

CO  Qj 

o 

ud 

cx  cx 

M  o 

CO  X 

<  CQ 

X  o 

Du  R 

uJ  CO 

t!  H 

U  CJ 

Cu  R 

Q  X 

u 

=  CO 

CO  D 

CO 

<  CO 

n  X 

Q  O 

Ct:  M 

X 

R 

Q 

>  R 

CO 

CO  CO 

CO  < 

(X 

X  o 

X  X 

X 

CO  o 

R  R 

R  R 

R  O 

>  R 

■K 

■K 

•< 

■K 

FUTURE  CONSIDERATIONS 


I 


CD 

s 

o 

M 

I 

CD 

< 

S 

O 

M 

EH 

M 

P 


X 

Eh  CD 
M  S 

CD  Eh 


D 

o 

< 

Eh 

M 

X 

U 

H 

X 

W 

X 

EH 

X 

M 


CD 


w 

X 

EH 

EH 

2 


O 

P 

< 


u 


hH  kC 
X  M 


X 

o 

o 

EH 

M. 

EH  O 

p 

P 

<  X 

o 

X 

X 

EH 

M  H 

P  w 

m 

Eh 

W  P 

X 

-  CD 

o 

X 

X  X 

Cl, 

>• 

P 

CD 

EH 


CD 

X 

M 


CJ  CD 
W  CD 


a  EH 


Cm 


M  X 


HC  ■ 


PREDICT  MORE  ACCURATE  GRAVITY  QUANTITIES  BY 
ADDING  LONG  AND  INTERMEDIATE  WAVELENGTH 
INFORMATION  TO  TOPOGRAPHIC  TERRAIN  EFFECTS 


EFFECT  OF  TERRAIN  REPRESENTATION,  GRID  SPAQNG, 
AND  FLIGHT  AI.TITUDE  ON  TOPOGRAPHIC  ODRRECTIONS 


FDR  AIRBORNE  GRADIOMETRY 
by 

Klaus-PeCer  bchwarz 
Michael  G.  Sideris 
I.N.  Tzlavos 

The  University  of  Calgary 
2500  University  Drive  N.W. 

Calgary,  Alberta 
CAIiADA  T2N  1N4 

ABSTRACT 

The  fast  Fourier  transform  (FFT)  algorithm  is  used  to  efficiently  compute 
terrain  corrections  for  airborne  gravity  gradients.  The  formulation  of  the 
equations  is  given  in  detail,  deriving  the  spectra  of  the  gradient  components 
of  the  gravitational  tensor  directly  from  the  spectrum  of  the  gravitational 
potential.  The  terrain  is  represented  by  either  line  masses  or  by  prisms. 
Results  of  the  method  are  given  for  a  very  rough  digital  terrain  field  in  the 
Kananaskis  area.  Comparisons  of  the  terrain  corrections  in  common  points  for 
grid  spacings  of  0.1,  0.2,  0,3,  0.4,  0.5  and  1  km  are  made  using  the  line  mass 
and  prism  representation  of  the  terrain,  for  two  different  flight  altitudes. 
Results  indicate  that  for  a  flight  height  of  I  km  above  the  terrain,  a  0.5  km  x 
0.5  km  grid  of  elevations  is  adequate  for  an  RMS  accuracy  of  1  Eotvos,  while  to 
retain  this  accuracy  for  a  flight  height  of  0.6  km  a  grid  spacing  of  about  0.25 


km  is  required. 


EFFECT  OF  TERRAIN 
REPRESENTATION,  GRID  SPACING, 
AND  FLIGHT  ALTITUDE  ON 
TOPOGRAPHIC  CORRECTIONS  FOR 
AIRBORNE  GRADIOMETRY 


BY 

I.N.  TZIAVOS^  M.G.  SIDERIS  AND  K.P.  SCHWARZ 


THE  UNIVERSITY  OF  CALGARY 
DEPARTMENT  OF  SURVEYING  ENGINEERING 
CALGARY,  ALBERTA,  CANADA 


*  On  leave  from 
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Presented  at  "  The  I5th  Gravity  Gradiometer  Conference  ",  Colorado  Springs,  Feb.  11-13, 1987 


NEED  FOR  ACCURATF  TERRAIN  CORRECTIONS 


MAIN  SOURCES  OF  NOISE  ( instrument  independent) 

-  Topographical  noise 

Due  to  surface  topography 

Can  be  minimized  using  OEMs  and  DTMs 

-  Geological  noise 

Due  to  density  variations  in  the  upper  crust 
Additional  information  needed  for  its  estimation 

-  High  frequency  effects  "filtered  out"  by  flight  altitude 


REQUIRED  TERRAIN  CORRECTION  ACCURACY 
•  1  E  RMS  or  better  at  flight  altitude 
-  Depends  on 

Terrain  sampling  rate  (grid  spacing) 

Terrain  representation  (prsms  or  point  heights) 
Flight  altitude 

(Extent  and  topographic  features  of  the  area) 


FFT  EVALUATION  MOST  CONVENIENT 

•  Fastest  approach  for  large  data  files 

•  Availability  of  DEMs 

-  Homogeneous  coverage  in  results 

-  Specral  analysis  and  Interpretation 


OBJECTIVES 


DERIVATION  OF  CONVENIENT  COMPUTATIONAL  FORMULAS 

•  2D  convolution  integrals 

Fast  evaluation  using  2D  FFT  routines 

-  Relationship  between  spectra  of  prisms  and  point  heights 


ESTIMATE  THE  EFFECT  OF  GRID  SPACING 

-  0.1km  X  0.1km  grid  results  used  as  control  values 

-  0.2km,  0.3km,  ...,1km  grid  results  compared  to  those  of  0.1km 
grid  for  1.1km  flight  altitude  (Zf  r  1.1km) 


ESTIMATE  THE  EFFECT  OF  TERRAIN  REPRESENTATION 

-  Prism  representation  results  versus  point  height  results 

-  Comparison  for  a  grid  spacing  of  1km  and  Zf  =  1.1km 


ESTIMATE  THE  EFFECT  OF  FLIGHT  ALTITUDE 

•  Zf  =  1.1km  and  Zf  =0.6km  were  used 

•  Comparison  for  various  grid  spacings 


PROPOSE  THE  GRID  SPACING  AND  TERRAIN  REPRESENT¬ 
ATION  NECESSARY  FOR  A  IE  OR  BETTER  RMS  ACCURACY 
AT  FLIGHT  LEVEL 


COMPUTATIONAL  FORMULAS 


flight 

level 


sea 

level 


0  dx  dy  dz 

V(Xo,yD.Zo)  =  Gp  J I J sNx-Xp)2+(y-yp)2  (i) 

E  h  Is2  +  (z.Zo)2]1'2 


GRAVITY  TERRAIN  CORRECTION  AT  Pq 

«V{Xp,yp.Zo)  h  *  •  *0 

t(xp,yp,Zo)  - =  Gp  J I J - dx  dy  (2) 

flZo  E  0  [82  +  (z-Zo)2]3^ 


-  Series  expansion  of  (2)  around  Zav>  neglecting  orgers  higher 
than  2nd  ,  gives  in  convolution  form 


t(Xp,yp,Zo)  =  -  2rtGpZav 

+  Gp[tii(xp,yp)*ki(xp,yp)+h2(xp,yp)"k2(xp,yp)l  (3) 


( =  Bouguer  affect  +  effect  of  residual  topography ) 


(4) 


hi  (x,y)  =  h(x,y)  -  Xgv  ,  h2(x,y)  =  ( h(x,y)  -Zgv  ]  ^ 

l<i(x,y)  =  ■  (Zo  av) '  +y  *  +  (*0  •  *av)*  1^'^ 

k2(x,y)  =  1/2  [x2  +  y2  +  (2o  -Zav)®] 

-  3(Zo  -  Zav)2  /  [x2  +  y2  ♦  (z^  -Zav)^]  5/2 


-  FFT  evaluation  of  (3) :  t2(Xp>yp^o)  =  (5) 

TZ(u,v)  =  F{t2(Xp,yp,Zo}  =  -27cGpZav8(u,v) 

+  Gp[Hi(u,v)Ki(u,v)  +  H2(u,v)K2(u,v)]  (6) 

Hi  (u,v)  =  F{hi  (x.y)}  ,  H2{u,v)  =  F{h2(x,y)} 

K-j  (u,v)  =  F{ki  (x,y)}  =  -2ite*2x<l(Zo  -  Zav)  (7) 

K2(u,v)  =  F{k2(x,y)}  =  -2jt^qe‘2x‘l(Zo  -  *av) 


GRADIENT  TENSOR  COMPONENTS  AS  DERIVATIVES  OF  OR  t 


«zx(x.y)  =  F  {TZX(u.v))  ,  TZX(u.v)  =  2iiiuTZ(u,v) 


tzy(x,y)  =  F-1{TZY(u.v)}  , 

tzz()'.y)  =  F*1{TZZ(u,v))  . 

txx(x.y)  =  F-’{TXX(u,v))  , 

txy(x.y)  =  F  ■''{TaY(u,v)} 
tyy(x,y)  =  F-1{TYY(u,v)}  , 


TZY(u,v)  =  2xivTZ(u,v) 

TZZ(u,v)  r  27rqTZ(u,v)  (8) 

TXX(u,v)  =  -  2jtu2q-lTZ(u,v) 
TXY(u,v)  =  -  2nuvq*'lTZ(u,v) 
TYY(u,v)  =  -  2itv2q-lTZ(u,v) 


TXX(u,v)  +  TYY(u,v)  +  TZZ(u,v)  =  0  or  -u2.v2  +  q2  =  o  (9) 


•50 


STATISTICAL  ANALYSIS  OF  TOPOGRAPHY  AND  TERRAIN 
CORRECTIONS  IN  THE  36  km  x  56  km  TEST  AREA  OF  ROUGH 
TOPOGRAPHY  USING  A  0.1  km  x  0.1km  GRID  OF  POINT  HEIGHTS 


Stati¬ 

stics 

Topo 

Ag  corr. 

Corr.  for  tij  at  zo 

=  4.5  km(Zf=  1.1km) 

graphy 

z=hp 

z=zo 

txx 

txy 

txz 

tyy 

tyz 

tzz 

metefs 

mgal 

Eotvos 

min 

1378.0 

-313.1 

46.6 

-76.7 

■40.7 

-74.9 

^5.9 

-62.8 

-129.7 

max 

3413.0 

-141.0 

14.6 

72.0 

38.9 

82.4 

70.5 

69.2 

86.5 

mean 

2129.6 

-226.7 

-15.1 

0.0 

0.3 

4.6 

0.8 

-3.5 

-0.8 

s.  d. 

350.2 

37.3 

12.2 

35.7 

14.6 

37.2 

18.5 

23.9 

45.5 

RMS 

2158.2 

229.7 

19.4 

35,7 

14.6 

37.5 

J 

18.5 

24.1 

45.5 

EFFECT  OF  GRID  SPACING 

(  USING  POIKT  HEIGHTS  AND  Zf  =  1.1km  ) 


RMS 

differences 

[E] 


^2X 


•  z-components  more  sensitive  than  the  horizontal 

-  Effect  on  the  horizontal  components  dependent  on  the  terrain 
frequency  content  along  the  x  and  y-direction 

-  1km  grid  insufficient  for  <  IE  RMS  accuracy 

-  0.5km  grid  sufficient  for  <  IE  RMS  accuracy 

-  0.25km  -  0.35km  grid  spacing  necessary  for  <  0.5  E  RMS 
accuracy 

•  Bias  always  <  IE  for  1km  grid 

<  0.35  E  for  0.5  km  grid 


EFFECT  OF  TERRAIN  REPRESENTATION 

( FOR  Zf  =  1 .1  km ) 


DIFFERENCES  BETWEEN  THE  RESULTS  FROM  PRISMS  AND  FROM 
POINT  HEIGHTS  AT  FOR  A  1km  x  1km  GRID 


Terrain  correction  differences,  in  E ,  for  tjj  at  Zf 

=  1.1km 

Statistics 

^xx 

^xy 

txz 

tyy 

tyz 

tzz 

min 

-6.2 

-1.5 

-7.0 

-1.6 

-2.1 

-7.3 

max 

6.5 

1.1 

8.1 

1.5 

1.6 

6.5 

mean 

0.4 

0.5 

0.0 

s .  d . 

3.5 

0.4 

3.5 

0.6 

3.6 

RMS 

3.5 

0.4 

3.5 

0.6 

3.6 

-  If  TZ(u,v)  is  the  spectrum  of  t2(x,y)  computed  from  prisms 

and  TZ(u,v)  is  the  spectrum  computed  from  point  heights , 
then 

_  sin(uiiAx)  sin(V7cAy) 

TZ(u,v)  - TZ(u,v)  (1 0) 

utcAx  VTcAy 

-  Prism  representation  has  smoothing  effect  on  high  frequencies 

•  x-components  more  affected  than  the  y-components  because 
this  particular  terrain  has  more  high  frequency  content  in  the  x- 
direction 

•  Differences  of  up  to  a  few  E  are  to  be  expected  in  the  more 
densitive  components 


EFFECT  OF  FLIGHT  ALTITUDE 

(U^!NG  POINT  HEIGHTS) 


RMS 


grid  spacing 
Ikml 


high  values  :  Zf  s0.6km  (ZQs4km) ;  low  values  :  Zf  sl.lkm(Zo=4.5km) 


-  Attenuation  stronger  for  z-components 

-  x-components  attenuated  more  than  y-components  due  to 
unisotropic  terrain  frequency  distribution 

-  Attenuation  independent  of  of  grid  spacing 

•  Accuracy  of  results  dependent  on  grid  spacing  /  altitude  ratio 
0.25km  grid  to  be  used  with  zf  =  0.6km  above  very  rough 
terrain  for  <  0.5  E  RMS  accuracy 


CONCLUSIONS 


TERRAIN  CORRECTIONS  NEEDED  TO  CORRECT  FOR  THE 
’•  TOPOGRAPHIC  NOISE  ”  ON  THE  MEASURED  GRADIENTS 


REQUIRED  ACCURACY  AT  FLIGHT  LEVEL  :  0.5E-1E 

-  Depends  on 

Grid  spacing 
Terrain  representation 
Flight  altitude 

-  Effects  on  horizontal  components  depend  on  terrain 
frequency  content  along  the  x  and  y-direction 


GRID  SPACING 

-  Affects  more  the  z-components 

-  Accuracy  decreases  with  Increase  of  grid  spacing 

TERRAIN  REPRESENTATION 

•  Point  and  prism  results  can  differ  by  a  few  E  RMS 

•  Prism  representation  is  more  accurate 


FLIGHT  ALTITUDE 

•  Stronger  attenuation  for  z-components 
-  Attenuation  effect  rather  insensitive  to  grid  spacing 


RECOMMENDATIONS 


FFT  TECHNIQUES  RECOMMENDED  FOR  THE  COMPUTATIONS 

-  Convenient  formula  derivations  in  the  specrai  domain 
•  Speed,  easy  handling  of  iarge  data  sets 

-  Homogeneous  coverage 

-  Frequency  domain  analysis  and  interpretation  of  data  and 
results 


FOR  Zf  =  1.1km  AND  Zf  =  0.6km,  0.5  km  AND  0.25  km  GRID 

SPACINGS  RESPECTIVELY  ARE  NECESSARY  FOR  <  IE 
RMS  ACCURACY 


AX  /  Zf  =  1  /  3  appears  to  be  a  good  choice 


ACCURACY  EFFECT  OF  HIGHER  THAN  2nd  ORDER  TERMS 
IN  THE  EQUATIONS  TO  BE  INVESTIGATED 


TOPOGRAPHIC  EFFECTS  IN  AIRBORNE  GRAVITY  GRADIOMETRY 


by 

Rene  Forsberg 
Geodaeclsk  InstltuC 
Gamlehave  Alle  22 
2920  CharloCtenlund 
Denmark. 


ABSTRACT 

The  gravitational  signal  due  to  terrain  masses  will  play  a  dominant  role 
in  airborne  gravity  gradiometer  surveys  over  mountainous  areas.  The  variations 
in  the  gravity  gradient  tensor  elements  may  easily  attain  magnitudes  of  several 
hundred  Eotvos  units,  thus  being  much  larger  than  typical  signals  associated 
with  possible  geophysical  structures  of  interest. 

To  smooth  the  gradient  field  and  enhance  “geological"  gravity  gradient 
anomalies,  the  gravitational  "noise"  caused  by  the  topography  may  be  attenuated 
using  available  digital  terrain  models,  from  which  the  elements  of  the  terrain  - 
induced  gradient  tensor  may  be  computed  efficiently  at  aircraft  altitude  using 
either  space  domain  (Integration)  or  frequency  domain  (FFT)  methods. 

In  the  paper,  both  of  these  computation  methods  are  outlined  and  compared, 
and  typical  magnitudes  of  the  effects  are  illustrated  by  examples  from  the  Rocky 
Mountains.  In  general,  statistics  of  the  gradient  variations  may  be  inferred 
rather  easily  from  emperical  ACF  parameters  of  the  topographic  heights,  and  the 
paper  concludes  with  a  number  of  such  examples  for  areas  of  different  types, 
from  lowlands  to  high  mountains,  yielding  useful  "hand  rules"  for  the  gravity 
gradient  terrain  effects  in  future  gravity  gradiometer  surveys  in  particular 


areas . 


TOPOGRAPHIC  EFFECTS  IN  AIRBORNE 
GRAVITY  GRADIOMETRY 


Rene  Forsberg 
Danish  Geodetic  institute 
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TERRAIN  REDUCTION:  removal  of  gravitational  signal  due 

TO  TOPOGRAPHY  IN  A  CONSISTENT  MANNER 


z 


o 


0  Geophysical  Exploration 

-  ENHANCE  GEOLOGICAL  SIGNAL 

-  EXAMPLE:  SIGNAL  10-20  E^  TOPOGRAPHY  50-100E 

0  Geodetic  Gravity  Field  Modelling 

-  MORE  SMOOTH  FIELD 

-  MORE  STABLE  DOWNWARD  CONTINUATION 

-  REMOVE-RESTORE  TECHNIQUE 

0  Local  Terrain  Effects  on  Second-Order  Gradients: 

-  MORE  SMOOTH  FIELD 

-  MORE  STABLE  DOWNWARD  CONTINUATION 

-  REMOVE-RESTORE  TECHNIQUE 

0  Local  Terrain  Effects  on  Second-Order  Gradients: 

-  very  large  on  ground^  requires  very  detailed  heights 

-  airborne  gradiometry:  attenuation. 


0  Modifications  For  Optimal  Gravity  Field  Modelling 
(Collocation  etc.) 

-  IMPORTANT  TO  AVOID  BIAS  IN  COMPUTED  TERRAIN 
EFFECTS.  ESPECIALLY  FOR  COMBINATION  SOLUTIONS 
WITH  TERRESTRIAL  GRAVITY  AND  DEFLECTIONS. 

-  BIAS  OCCUR  IN  Tyy.  AND  T^^  WHEN  COMPUTATION 

AREA  FINITE  AND  MEAN  HEIGHT  >  0! 

-  REMEDY:  REFER  TO  MEAN  HEIGHT  OR  USE  ISOSTASy! 

-  BETTER  REMEDY:  RTM  -  RESIDUAL  TERRAIN  MODEL 


h 

*^ref 
0 

RTM  ADVANTAGES: 

-  HEIGHT  DATA  NEEDED  FOR  SMALLER  AREA 

-  CONSISTENCY  OF  NEIGHBOURING  SURVEYS 

-  CONSISTENCY  IN  PRISM  AND  FFT  METHODS 


OTPUTATION  METHODS; 


PRISM  INTEGRATION  AND  FFT 


1)  PRISM  METHOD 


NUMERICAL  METHOD 


Approximate  formulas  (McMillan): 


-  6K( 


rs 


T  8  «»h(7  i-jlTjf ) 


2Mr^L 


rH 


T„  =  ■. 


Implementation:  detailed/coarse  grid^  innerzone  spline. 

AUTOMATIC  FORMULA  SELECTION 


:  PARKER  EXPANSION 


4 


z 


o 


ref 


f  Kv)  - 


U 

F:  Fourier  transformation^  <«=  (u^+v^) 
Gradients;  T-v  * 

I«  *  -'**T 

»  «-‘f 

T„  =  -iuwT 


Implementation  -Hj  h^j  ...  real  transforms 

-COMPLEX  FILTER  INVERSE  TRANSFORM:  PAIRWISE 
RESULTS 


ncTCR  fcOTvos  eorvos 

1000  2000  3000  -300  0  300  600  -300  0  300 


Profile  test:  33.3“N  (passing  2680  m  peak) 


ORAOIENTS  AT  3000  H  LEVEL 


LCOCNOt 

- TXT 

— —  TTl 
- TTY 


OISTWICe  I  Kill 


R.M.S.  FIT  PRISM 


fft:  t^^:  1.9E,  T^y'  ('^  =  8) 


EXAMPLE:  CENTRAL  COLORADO 


Table  1.  Comparison  test  -  '4^1  points  (Eotvos) 


^xx 

Tyy 

^zz 

^XY 

Txz 

Tyz 

COMPUTED  MEAN 
GRADIENTS 

X 

-^.9 

-n.6 

1.4 

-.9 

-4.9 

1.0 

cr 

3G.0 

34.5 

57.0 

20.3 

43.2 

41.0 

RMS  DIFFERENCE 
FFT  -  PRISM 

Q.5 

0.5 

0.9 

0.5 

.8 

1.0 

(N=8) 

Table  2.  Dependence  on  no,  of  FFT  expansion  terms 


T  T 

R.M.S.  DIFFERENCE  ' ZZ  XY 


N 

=  1 

22.2 

8.3 

N 

=  3 

5.5 

1.6 

N 

=  8 

0,8 

1.1 

N 

=  15 

0.9 

0.6 

COVARIANCE  AflALV$IS 


0  Topography  ACF  Gravity  ACl  - (first  order  parker) 

Gradient  ACF  at  altitude 

0  Selected  topography  ACF:  SINGLE  DIFFERENCE  LOGARITHfilC 
(Forsberg,  1986) 

Consistent  gravity  ACF  -  closed  formulas  both  on  ground 
and  aloft.  E.G, : 

Result  at  elevation  H: 

<S'^W)  ~  (2.ir&f )  S'h  ((T)+2H>*  ~  CF^Ifn^P 

T  FIXED  D  function  OF  X;^  FOR  TOPOGRAPHY 


0  Typical  good  fit:  T  =  10  km^  d  =  0.5  km  (Jotunheimen^  Norway) 


NORMflLIZEO  covariance 


DISTANCE  fKH) 


0  Analysis  for  sample  areas  (FFT)^  1°rtm 

— ►  AT  600  M  CLEARANCE  ELEV, 

0  Relationship  of  variances: 


VAR 

=  VAR 

-H 

-< 

-< 

II 

8  VAR 

(T,,) 

VAR 

'Txz' 

=  VAR 

1 

(T,,)  = 

1 

-  VAR 
9 

4. 

(Tzz) 

VAR 

(IxY^ 

1 

=  Q  VAR  (T  ) 
^  ZZ 

Colorado 

R.M.S.  (E) 

^xx 

Tyy 

— 1 

M 

'*'xY 

Txz 

T« 

Actual 

:5b 

A3 

Ai 

Model 

38 

38 

62 

22 

AA 

AA 

Topographic  effects  effectively  handled  by 

RTM  -  REDUCTION  "REMOVE-RESTORE" 

Effects  large  in  Alpine  regions  (up  to  200-300E) 

Computations  by  FFT  or  prism  method 
(High-order  expansion  necessary  in  FFT) 

Easy  magnitude  estimation  from  topography  ACF, 

-  relatively  large  effects  even  for  lowlands  (due 
lower  effective  flight  level) 


T^ITLF:  _OF_R\P^R:  Topographic  "‘‘focts  in  Airborne  Gravity  G  rad  Lomet  ry 
SPEAKER :  Rene  Forsberg 

0JL;E  S_T I q VS  AVD  CO>PiEVTS: 

1 .  Quest  ion :  Carl  Bow in 

It  appears  that  your  conversion  of  topography  to  gravity  assumes  the  lack  of 
compensation  of  the  topogranhv ,  hence  the  predicted  gravity  may  be  larger 
than  observed. 

Response : 

Yes,  but  the  computed  gravitv  correlation  is  like  a  Bonguer  correction. 

(Carl  Bowin:  I'.'e  do  not  measure  a  Bouguer  anomalv,  onlv  the  equivalent  t’'’’ 
free-air  anomaly). 


OBSERVABILITY  OF  LAPLACE’S  EQUATION  USING 
A  TORSION-TYPE  GRAVITY  GRADIOMETER 


by 

Alan  H.  Zorn 

Dynamics  Research  Corporation 
60  Frontage  Road 
Andover,  MA  01810 


ABSTRACT: 

The  trace  of  the  gravitation  gradient  tensor,  theoretically  zero,  is  not 
directly  observable  by  a  torsion-type  gradiometer.  However,  changes  in  the  trace 
are  directly  observable  at  track  crossings  by  a  moving-base  torsion  gradiometer. 
The  practicality  of  measuring  these  trace  changes  is  illustrated  by  covariance 
analysis  of  CGSS  track  crossings.  The  sensitivity  of  these  results  to  flight 
conditions  and  noise  model  parameters  is  also  presented. 
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TENSION  GRADIOMETER  MEASURES  ALL  SIX 
COMPONENTS  OF  SYMMETRIC  GRADIENT  TENSOR 

TORSION  GRADIOMETER  MEASURES  ONLY  FIVE 


DEFINITIONS  AND  TERMINOLOGY 
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ABSTRACT 


OBSERVABILITY  OF  LAPLACE’S  EQUATION 
USING  A  TORSION-TYPE  GRAVITY 
GRADIOMETER 


ALAN  ZORN 

DYNAMICS  RESEARCH  CORPORATION 


The  trace  of  the  gravitation  gradient  tensor,  theo- 
reticaliy  zero,  is  not  directly  observable  by  a  torsion- 
type  gradiometer.  However,  changes  in  the  trace  are 
directly  observable  at  track  crossings  by  a  moving- 
base  torsion  gradiometer.  The  practicality  of  measuring 
these  trace  changes  is  illustrated  by  covariance 
anaiysis  of  GGSS  track  crossings.  The  sensitivity  of 
these  results  to  flight  conditions  and  noise  model 
parameters  is  also  presented. 
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TITLE  OF  PAPER:  Observability  of  Laplace's  Equation  Using  a  Torsion-Type 
Gravity  Gradiometer 


SPEAKER:  Alan  H.  Zorn 

QLT:  ST  ION'S  AND  COMMENTS: 


1.  Question:  Dave  Sonnabend 

Have  you  taken  into  account  the  errors  in  estimating  the  angular  rate 
corrections  to  the  gradiometer?  For  instance,  gyro  bias  causes  correlated 
errors  in  .ncasuring  A. 

Response : 

Rate  measurements  are  good  enough  for  instruments  used  in  aerial  surveys. 

2.  Comment :  Warren  Heller 

Dave's  point  was  that  your  technique  is  a  good  way  to  test  for  gradiometer 
errors  but  that  the  instrumentation  is  not  adequate  for  tests  of  inverse 
square  law. 

Response : 

I  am  not  advertising  ray  idea  as  a  test  of  the  inverse  square  law. 

3.  Comment :  Ho  Jung  Paik 

Your  signal  tor  inverse  square  law  test  Is  of  the  order  of  GM/R^,  which  is 
about  0.2  E/kra.  So  the  sensitivity  of  the  gradiometer  must  be  improved  by 
at  least  a  factor  of  10  before  you  can  have  a  test  of  the  inverse  square  law. 

Response : 

t 

My  simulation  of  the  Texas-Oklahoraa  area  shows  that  the  signal  is  about  5-10 
E/km. 

A .  Question :  Dan  Long 

1  suggest  there  is  more  likelihood  that  dX 

7z 

or  ^  would  show  interesting  results. 
dE 

Response : 


It  is  harder  to  measure  dX,  but  I  see  your 

dZ 


rather  than  ^ 
dS 


5.  Question: 


Alan  Rufty 


Noise,  as  measured  by  the  gradiometer,  does  not  form  a  conservative  field. 
How  do  you  separate  out  the  signal  from  the  noise  so  that  your  assumption 
of  a  conservative  field  can  be  used? 

Response : 

If  one  knows  something  about  the  nature  of  the  noise,  then  one  can  apply 
the  conservative  field  assumption  only  to  the  signal.  I  do  not  assume  that 
the  noise  is  conservative. 

6.  Question:  George  F.  Hinton 

Shouldn't  your  method  also  be  sensitive  to  the  Eotvos  effect? 

Response : 

No,  I  am  taking  this  into  account  with  my  measurement  of  the  velocity. 


RATIOS  OF  GRAVITY  GRADIENT.  GRAVITY.  AND  GEOID 
FOR  DETERMINATION  OF  CRUSTAL  STRUCTURE 


by 

Carl  Bowln 

Woods  Hole  Oceanographic  Institution 
Woods  Hole,  MA  02543 


ABSTRACT 

It  has  been  shown  (Bowin,  1983,  1985;  Bowin  et  al.,  1986)  that  the  simple 
ratio  of  gravity  ^o  geoid  anomalies  at  the  center  of  a  feature  provides 
information  comparable  to  that  obtained  by  the  Fourier  transforming  of  either 
field,  as  well  as  to  that  obtained  by  other  traditional  spatial  and  frequency 
analysis  methods.  If  perfectly  accurate  measurements  of  any  one  of  the  three 
fields  -  geoid,  gravity,  or  vertical  gravity  gradient  -  were  available  over  the 
entire  earth's  surface,  no  further  observations  would  be  necessary  or  useful. 
The  other  fields  could  be  derived  completely  by  either  integration  or  differ¬ 
entiation.  However,  both  caveats  fall  us  in  the  real  world:  the  existing  data 
are  neither  perfectly  accurate  nor  universal  in  coverage.  Thus,  for  the 
immediate  future,  knowledge  of  the  full  spectrum  of  the  Earth's  gravity  field 
will  come  from  combinations  of  the  various  measurement  types.  In  this  present¬ 
ation  we  summarize  our  ratio  results  to  date,  and  indicate  the  utility  of 
gravity  gradient  measurements  to  aid  the  determination  of  crustal  structure  and 
depth  of  causative  mass  anomaly  sources- 


Note:  The  following  references  include  most  of  the  illustrations  used  In 
Carl  Bowin's  talk  at  the  15th  Gravity  Gradiometry  Conference.  Also 
attached  are  copies  of  the  new  world  residual  geoid,  gravity  and 
vertical  gravity  gradient  maps  that  were  presented  at  the  Conference 
(Figure  1,  2  and  3  respectively ) - 

References : 

1)  Bowin,  Carl,  1983.  Depth  of  Principal  Mass  Anomalies  Contributing  to  the 
Earth's  Geoidal  Undulations  and  Gravity  Anomalies.  Marine  Geodesy,  V.7, 

p. 61-100. 

2)  Bowin,  Carl,  1985.  Global  Gravity  Maps  and  the  Structure  of  Earth. 

IN;  W.J.  Hinze ,  ed . ,  SEG  Volume;  Utility  of  Regional  Gravity-  and 
Magnetic-  Anomaly  Maps. 

3)  Bowin,  Carl,  1986.  Topography  at  the  Core  Mantle  Boundary.  Geophysical 
Research  Letters,  Vol.  13,  No.  13,  p. 1513-1516. 

■4)  Bowin,  C.,  Edward  Scheer,  Woollcott  Smith,  1986.  Depth  Estimates  from 
Ratios  of  Gravity  Geoid,  and  Gravity  Gradient  Anomalies.  Geophysics, 

Vol.  51,  No.l,  p . 123-136. 


COMBINING  GRAVITY  GRADIOMETRY  WITH  OTHER 
EXPLORATION  METHODS  FOR  GEOPHYSICAL  PROSPECTING 

by 

Dr.  Anthony  A.  Vassiliou 
Department  of  Surveying  Engineering 
The  University  of  Calgary 
2500  University  Drive  N.W. 

Calgary,  Alberta 
CANADA  T2N  1N4 

ABSTRACT 

The  Parker-Oldenburg  algorithm  for  fast  computation  of  potential  field 
anomalies  is  modified  to  allow  for  multilayer  inversion  of  gravity  data.  First 
the  algorithm  is  developed  for  the  inversion  of  gravity  anomaly  data  and  then  it 
is  further  extended  to  suit  airborne  gravity  gradiometer  data.  It  is  shown  that 
the  use  of  gradiometer  data  is  preferable  to  the  use  of  gravity  anomaly  data  for 
the  computation  of  the  anomalous  density  and  topography  of  subsurface  densities. 
The  solution  of  the  inverse  gravity  problem  is  constrained  by  density  and  if 
possible  layer  depth  Information.  Density  information  can  be  obtained  from 
borehole  surveys  or  from  computed  compresslonal  seismic  velocities  via  a  non¬ 
linear  formula.  Subsurface  layer  information  can  be  obtained  from  inverted 
seismic  reflection  data.  In  addition  to  gravity  gradiometer  data,  the  same 
inversion  algorithm  can  by  employed  to  determine  magnetic  susceptibility  for 
subsurface  magnetized  layers.  Using  aeromagnetlc  data  the  depth  of  these  layers 
is  determined  by  using  the  Spector-Grant  spectral  method. 


COMBINING  GRAVITY  GRADIOMETRY  WITH  OTHER 
EXPLORATION  METHODS  FOR  GEOPHYSICAL  PROSPECTING 


Anthony  A.  Vassiliou 

Department  of  Surveying  Engineering 
The  University  of  Calgary 
2500  University  Drive  N.V. 
Calgary,  Alberta 


Paper  presented  at  the  I5th  Gravity  Gradiometer  Conference 
Colorado  Springs,  U.S.A. 

February  11-12,  1987 


Objective : 


Develop  a  method  to  determine  the  anomalous  density  and 
depth  of  subsurface  layers  using  gravity  gradiometer  data 
and  incorporating  information  from  other  geophysical 
prospecting  methods. 


Content ; 

•  Forward  Parker's  algorithm. 

•  Single  layer  inversion  using  gravity  anomaly  and  gravity  gradiometer 
data. 

•  Multilayer  inversion  with  gravity  anomaly  and  g.radiometer  data. 

•  Constraints  on  the  multilayer  inversion  from  other  geophysical 
exploration  methods. 


1. 


Forvard  Parker's  algorithm 
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Fourier  transform  of  {  } 

Newton's  gravitational  constant 
density 

altitude  at  which  gravity  data  are  observed 

layer  topography 

anomalous  gravity  potential 

second-order  gradients  of  T 

spatial  frequencies  in  the  two  directions 

(u“  +  v“)  ' 


then : 
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-2‘nqz  oo 
2irGe  °  I 
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---qz. 
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z  z  '  ,  n : 
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-27qz  * 

F{T  (x,v)}  =  -2"Ge  °  Z  — ilS -  Ztruv  P{p  (x,y)h^(x ,y)  } 

XV  '  ,  n! 

n=l 


-27qz 


n-1 


F’T  rx.vD  =  -27Ge  °  Z  -  j  27u  F{p  (x.y)h"(x,v')  } 

xz  •  ,  n! 


n*l 


■2nqz  »  sn-1 


F{T  (x,y)}  =  -27Ge  °  Z  -  j  2'tv  F{p  (x  .y^h'^fx ,  v)  } 

vz  ,  n! 


n«l 


j  is  the  imaginary  unit. 


•  Assume  gravity  anomaly  data  are  available  and  have  been  corrected  for 
any  thin  sediments  or  water  overlying  the  layer.  Also  assume  that  the 
density  varies  in  some  smooth  fashion  reasonably  known.  Then  the  layer 
topography  h(x,y)  is  determined  from 

iTTqz^  »  n  I 

P{o  (x  ,y)h(x,y)  }  =  -  - E  -  F{d  (x  .y^h^Cx.y)  } 


•  Use  equation  (1)  in  an  iterative  inversion.  Start  with  some  estimate 
of  h(x,y)  (probably  given  from  borehole  surveys),  compute  R.H.S.  of  (O 
up  to  n=N,  then  compute  L.H.S.  and  take  its  inverse  2-D  FFT  to  get 
hfx,y).  Iterate  to  achieve  convergence.  The  degree  n=N  up  to  which  the 
sum  in  the  R.H.S.  of  (1)  is  computed,  is  determined  from  the  ratio 
(S  ,/S»)<l  .0x10  ^  (with  S  =  (2irq)''  Vn!  F{h^(x,y)p(x,y) }) .  Convergence 
in  the  iterative  algorithm  is  considered  when  the  r.m.s.  of  the 
differences  between  two  consecutive  estimates  of  gridded  h(x^,y^)  is 
smaller  than  0.3  m. 


•  Excessive  noise  amplification  in  downward  continuation  is  prevented 

by  low-pass  filtering.  For  gravity  anomaly  profiles,  it  has  been 

suggested  by  Oldenburg  to  use  cutoff  passband  frequency  u^ : 0 . OSSu^^O . 1 

(cycles/profile  spacing)  and  edge  of  the  stopband  critical  frequency 

O.lSSu  SO. 25. 
s 

•  For  gravity  gradiometer  data  higher  limits  for  cutoff  frequency  u^ 
and  stopband  critical  frequency  u^  have  to  be  used,  due  to  their  high 
frequency  signature.  Therefore  less  information  will  be  smoothed  out  in 
the  downward  continuation. 


•  The  iterative  inversion  algorithm  expressed  by  equation  (1)  converges 
onlv  for  densities  larger  than  a  minimum  density  P  ,  .  It  also 
converges  for  heights  z^>  minimum  z^.  In  short,  convergence  depends  on 
the  density  contrast,  depth  of  the  subsurface  layer  and  the  passband; 
ptopband  frequencies  of  the  low-pass  filter. 


Single  Laver  Gravity  Inversion  (continued^ 


•  Using  gravity  gradiometer  data,  we  have  six  data  sets  instead  of  one 
to  determine  the  subsurface  layer  topography.  Thus  the  determination  of 
h(x,y)  is  strengthened  due  to  the  additional  information,  and  thus  the 
solution  becomes  more  reliable.  The  equations  for  the  iterative 
inversion  of  h(x,y),  using  gradients  1..  are  derived  from  the  forward 
Parker’s  algorithm  in  the  same  way  as  for  the  gravity  anomaly  data. 


Gravity  Gradient  Single  Layer  Inverston 

•  For  simplicitv  assume  T  (x,v),  T  (x,v)  ,  T  (x,v)  data. 

X2  •  yz  '  zz  • 

•  From  the  forward  Parker's  algorithm,  after  reconstructing  the 
equat ions 

Imqz 

F{o  ^  fx  ,y)hj  (x,y)  }  =  -  (F{T^^(x,y)  }/2iTG]e 


r  (-Ttq) 


jZTTu  F{d  j  (x,y)h^(x,y)  } 


F{Oj(x.y)h^(x,y)}  =  -  [F{Ty^(x.y) }/2mG]e 


-  j2Trv  F{p  (x,y)h"(x,y)} 

n :  11 


iTtqz 

Ffo^(x.y')h^(x.y)}  =  -  [F{T^^  (x  .y) } /2rG  ]  e 

S  FCo,  (x,y)h"(x,>0  } 

n!  1  1 


•  Iterate  for  each  one  of  the  data  sets,  to  finally  compute 

three  sets  of  subsurface  layer  topography  h^(x,y).  Then  combine  the 

solutions  for  h.(x,y)  from  T  ,T  ,T  using  spectral  weighting 
L  X  Z  ^  Z  Z  Z 


F{h^(x,y);  = 


^  (u,v)F{hj  j(u,v)}  +  ^(u,v)P{hj  ^(u.v)}  +  2(u,v)F{h^  2(u,v)} 


S  (u,v)  +  S  (u,v)  +  S  (u,v) 

^.1  ‘'l.Z  ''l.S 

where  h,  h,  h,  ,  ...  solutions  for  h,(x,y)  from  T  ,  T  ,  T 

1,1  1,2  1,3  1  xz  yz  zz 

respectively 

S,  ;  S,  ;  are  the  corresponding  amplitude  spectra. 

^,1  N.2  ^.3 


This  iterative  weighting  procedure  can  be  applied  for  more  than  three 


Multilayer  Gravity  Inversion 


'  For  simplicity  ve  use  gravity  anomaly  data.  Assume  a  4  layer  model 
on  top  of  which  a  layer  (composed  of  sediments,  water)  with  density  c 

u 

is  located.  The  density  and  bedthickness  of  layers  1,2,3  are  known. 
Determine  the  topography  of  the  bottom  layer  h(x,y)  having  some 
knowledge  of  the  density  P. 

•  Use  Parker's  formula  for  more  than  one  layer 

4  -2ifqz.  *  .n-1 

P{Ag(x,y)^  -  2^G  L  e  E  -  F{a  (x,y)  h"(x,y)}  (2) 

k-1  n-1  n.  X  k 


(3) 


Multilayer  Gravity  Inversion  (continued) 


with  Ag(x,y)  determined  from 
others 

4  -2Trq2  »  .n-1 

P{Ag(x,y)}  =  -2itG  Z  e  Z  -  F{  o,  (x ,  y)h!^(x ,  y) }  (4) 

others  k=2  n=l 

•  Equation  (3)  is  used  in  an  iterative  manner  to  solve  for  the  unknown 
topography  h^(x,y).  Some  apriori  (or  model)  knowledge  about  h^(x,y)  is 
necessary.  This  ':nowledge  will  provide  the  depth  z^. 

•  The  depth  of  each  layer  is  then  determined  by 

2,  ,  =  2,  -  (bed  thickness),  ,  (5) 

k+ 1  k  k+1 


•  The  values  of  hj(x,y)  vary  at  each  iteration,  thus  the  depths  need 
to  be  updated.  Due  to  this  updating,  the  gravity  anomaly  Ag(x,y)  will 

others 

vary  as  well.  A  reasonable  starting  model  for  Ag(x,y)  would  prob  .bly  be 

others 


Ag<'x,y)  =  Ag(x,y) 
others  observed 


:  a.  (x,y) 
k=l 


(6) 


•  Combining  equations  (3),  (4),  (5),  (6)  an  iterative  procedure  for 

computing  the  topography  hj(x,y)  can  be  employed. 


Mulcilaver  Gravltv  Gradient  Inversion 


•  Assume  for  simplicity  that  onl  '  T  ,  T  ,  T  are  available. 

■  xz  yz  zz 

•  From  the  multi-layer  formulation  of  the  forward  Parker's  algorithm  ve 
get  after  reconstruction  of  the  formulas 

2r:qz 

F{a,(x,y)h  (x,y)}  =  -  [P{T  (x,y)  -  T  fx,y)}/2rG]e 

observed  others 


I  -  jP-ru  P{0j(x,y)hj(x,y)} 

n=2 


P{Oj (x,y)hj (x,y) ;  = 


-  [P(Ty^(x.y)  -  Ty^(x,y)}/2TtG]e 

observed  others 


j2TTV  F’C  j(x,y)h^(x,y) } 


P{o,(x,y)h  (x.y)}  =  -  [P{T  (x,y)  -  T  fx,y)}/2rG]e 

4  ^  Ac  C  C 

observed  others 


I  X±I£2_  if(0  (x,y)h”(x,y) } 
n=2  ^  ^ 


where  T  (x,v),  T  (x,v),T  (x.y)  are  the  gradients  T  T  ,T 

xz  '  yz  •  xz  ®  xz ,  yz  zz 


others  others  others 


computed  from  the  ocher  three  layers, 


•  The  same  Iteration  procedure  as  for  the  gravity  data  has  to  be 

applied  for  gravity  gradient  data.  To  combine  the  solutions 

h,  •,h,  ,,;h,  ,  for  the  subsurface  topography  from  T  ,T  ,T 
1,21,21,3  r  o  r  .  xz  yz  zr 

respectively  we  use  a  spectral  weighting  scheme  as  in  the  single  layer 


case . 


Multilayer  Gravity  Gradient  Inversion  (continued) 


•  After  stripping  the  overlying  layers,  single  layer  inversion  has 
essentially  to  be  performed.  Deeper  layer  inversion  can  then  be 
attempted  by  examining  a  lower  frequency  band. 

•  Again  the  advantage  of  using  gravity  gradiometer  data  is  that  five 
more  data  sets  are  available  with  higher  frequency  content.  Thus  the 
solution  for  h^(x,y)  is  more  reliable  and  more  detailed  for  shallow 
depths . 


Constraints  on  Multilayer  Gravity  Tr.version  fron  Other  Geophysical 


Methods 


•  The  Inverted  surface  need  not  be  the  bottom  one,  as  long  as 
Information  about  the  other  layers  is  available. 

•  Densities  from  the  other  layers  can  be  obtained  from  seismic  data,  or 
from  borehole  surveys.  Borehole  data  will  provide  density  information 
at  two,  or  three  points  (maybe  a  few  more)  in  the  area  of  interest. 
Then  the  densities  between  these  points  can  be  assumed  to  vary  linearly. 

•  Information  about  bedthickness  for  each  of  the  other  layers  can  be 
obtained  from  seismic  reflection  data  or  from  well  logs  or  from  sonic 
logs. 

•  If  detailed  information  about  the  shape  of  each  layer  is  not 
available  (as  is  often  the  case),  then  a  reasonable  approximation  is 
chat  the  topography  of  the  other  layers  varies  in  the  sane  way  as  the 
topography  of  the  inversion  subsurface. 

•  For  the  gravity  correction  due  to  the  sedimentary  layers,  the  density 
(for  Western  Canada),  is  taken  from  extensive  density  legging  as 

3 

p  =  2.266  +  0.146  z  t  0.210  (gr/cm  )  (z;  depth  of  sediments) 

•  .Assuming  known  conpressional  (P-wave)  velocities  with  depth,  from 
seismic  traveltime  inversion,  a  reasonable  density  approximation  can  be 
obtained  from  the  non-linear  formula 

0 .25 

p  “  0.23  V  '  (V  :  P-wave  seismic  velocity) 

•  T nfomatior.  about  densities  can  also  be  obtained  from  deconvolution 
of  reflection  seismograms  and  computed  reflection  coefficients.  In  this 
case  the  density  or  velocity  for  some  layer  should  be  known  (not  the 
upper  sedimentary  weathering  layer  though). 


Magnetic  Daca-Relation  to  Gravity  Data 


•  For  a  constant  ratio  of  density  P,  to  magnetic  susceptibility  k 

■T.  ° 

T  =  —  :  constant 

k 

the  magnetic  anomalous  potential  W  is  related  to  the  first  order 
gradients  of  the  anomalous  gravity  potential  T  by  the  Poisson's  equation 

with  :  magnetic  field  strength 

:  magnetic  permeability 


So 


7T 


T  u  G 
0 


•  A  similar  relation  can  be  developed  between  the  magnetic  field 
anomaly  3  (B  =  -  7W)  and  the  second-order  gradients  of  T.  Inversion  of 
magnetic  anomaly  data  resembles  very  much  the  inversion  of  gravity 
gradiometer  data. 


•  The  anomaly  in  the  direction  of  the  perturbing  field  'vb!  is  related 
to  the  magnetic  susceptibility  by  the  relation 


F{ ! abI 


ZmNCu ,v) q  ( H 


-Brqz 


,  n-I 


I  Ilia} -  F{k(x,v)h"(x,v)} 

,  n ! 


n=l 


where  the  dimensionless  vector  N(u,v)  is  a  function  of  magnetic 
declination  D  and  magnetic  inclination  I. 


•  Information  about  the  depths  of  the  magnetized  layers  can  be  obtained 
from  the  PSD  of  aeromagnetic  intensity  data,  by  using  the  method 
developed  by  Spector  and  Grant. 


Conclusions 


•  Forward  Parkers *s  algorithm  can  be  used  to  invert  gravity  anomaly  and 
gravity  gradiometer  data  in  an  iterative  mode. 

•  By  using  gravity  gradiometer  data,  a  solution  can  be  obtained  for  the 
shallower  layers  with  better  detail  than  for  deeper  layers  for  which 
gravity  anomaly  data  have  to  be  used. 

•  Knowledge  of  depth  and  density  of  the  overlying  layers  constrains  the 
solution  for  the  layer  topography. 

•  Depth  and  density  information  about  the  subsurface  layers  can  be 
obtained  from  either  borehole  surveys  or  inverted  seismic  reflection 
traveltime  data,  or  from  density  logs  or  from  sonic  logs. 

•  The  same  algorithm  as  for  gravity  anomaly  and  gravity  gradient  data 
can  be  used  for  magnetic  anomaly  data,  thus  yielding  magnetic 
susceptibility  values.  Information  about  the  magnetized  layer  depth  can 
be  obtained  from  analysis  of  aeromagnetic  data. 


COMPUTATION  OF  THE  GRAVITY  VECTOR 
FROM  TORSION  BALANCE  DATA  IN  SOUTHERN  OHIO 


by 


Dr-  D.  Araoelos 

Department  of  Geodesy  and  Surveying 
University  of  Thessaloniki 
Greece 

Mr.  Christian  Tscherning 
Geodaetisk  Institut 
Geodetic  Department  I 
Gamlehave  Alle  22 
2920  Charlottenlund 
Denmark 


ABSTRACT 

The  method  of  least  squares  collocation  was  used  for  gravity  vector 
estimation  from  torsion  balance  measurements  in  Southern  Ohio.  The  results 
were  dependent  on  the  covariance  function  used  and  on  the  selection  of  a  proper 
signal  to  noise  ratio.  Here  a  standard  deviation  of  the  noise  of  10  E.U. 

(10"'^  *  s”^)  gave  the  best  results.  Expressed  in  terras  of  standard  deviations 
of  observed  minus  predicted  differences,  the  best  results  were  0.4"  for  deflec¬ 
tions  of  the  vertical  and  2  mgal  for  gravity  anomalies.  This  compares  to  signal 
standard  deviations  of  4"  and  22  mgal  for  deflections  and  gravity  anomalies, 
respectively. 
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Abstract:  The  method  of  least  squares  collocation  was  used  for  gravity 
vector  estimation  from  torsion  balance  measurements  in  Southern  Ohio. 
The  results  were  dependent  on  the  covariance  function  used  and  on  the 
selection  of  a  proper  signal  to  noise  ratio.  Here  a  standard 
deviation  of  the  no«ise  of  10  E.  U.  <10  *s  “"I  gave  the  best  results. 
Expressed  in  terms  of  standard  deviations  of  observed  minus  predicted, 
differences,  the  best  results  were  0.4'*  for  deflections  of  the 
vertical  and  £  mgal  for  gravity  arK:malies.  This  compares  to  signal 
standard  deviations  of  4"  and  22  mgal  for  deflections  and  gravity 
anomalies,  respectively. 


1.  Introduction. 

The  method  of  collocation  and  similar  optimal  estimation  techniques 
have  been  proposed  as  possible  methods  for  gravity  vector  estimation 
from  gravity  gradiometer  data,  see  e. g.  Jekel i ( 1S35) .  Since  gravity 
grad  IOmeter  data  are  not  yet  available,  the  estimation  procedures  have 
only  been  tested  using  simulated  data.  However,  the  torsion  balance 
deliver's  measurements  of  the  same  kind  as  the  gradiometer',  namely 
linear'  combinations  of  second  or'der  derivatives  of  the  gr'avity 
potf.'ntial  (W). 

Using  torsion  balance  data  we  may  then  demonstrate  how  least— squares 
collocation  (USC)  can  bo  used  for  gravity  vector ■ est imat ion,  and  also 
show  which  type  of  difficulties  one  may  encounter  while  using  the 


riiet  hod  . 


I  ri  or'dei'  to  dt'aw  conclusions  from  the  use  of  torsion  balance  data  to^^ 
the  use  of  gr£«vity  gr'ad  i  oniet  et'  data,  the  data  d  i  st  t' i  but  i  on  should 
cor'r'espond  to  the  orie  planned  for  the  gradiometer.  Futher'more  gr-avity 
vector  data  must  also  be  available  in  the  same  a»'ea  as  the  tor'sion 
balance  data.  Such  an  optimal  situation  is  found  in  a  0. £5  deg.  ♦  0.75 
deg.  area  in  Southern  Ohio,  where  gravity  data,  deflect  io*ns  of  the 
vertical  and  torsion  balance  data  were  observed  and  collected  in  the 


period  p»'ior  to  1370.  The  data  is  somewhat  sparser  than  the  one 
planned  for  gradiometry,  and  the  area  covered  is  smaller  than  the 
typical  area  forseen  to  be  covered  by  so>me  hours  of  aerial  gradiometer 
measurements-  Furthermore  the  data  is  not  observed  at  altitude,  but  at 
ground  level.  Also  Southern  Ohio  is  not  typical,  since  the  area  has  no 
large  height  variations.  However,  in  an  operational  situation, 
topographic  effects  will  in  general  have  to  be  removed  (and  restored). 
In  this  way  all  areas  from  the  standpoint  of  the  gravity  field 
variation  will  look  like  Ohio  (Forsberg,  1984,  table  fc).  Also  the 
gravity  variation  at  an  altitude  of  0.5  -  £. 0  km  will  not  be  much 
different  from  the  variation  at  ground  level,  when  topogi'aphic  effects 
have  been  removed.  Furthermore,  as  we  shall  see,  an  increase  in  the  ^ 
data  density  would  probably  not  have  improved  the  quality  of  the 
prediction  vector. 


Deflection  of  the  vertical  estimation  have  been  successfully  tested  in 
Southern  Ohio  by  Badekas  and  Mueller  (19&3>,  using  a  simple  numerical 
integration  method.  Mith  LSC  we  should  be  (and  have  been)  able  to 
obtain  similar  results,  the  advantage  being  that  the  method  not  only 
enables  the  estimation  of  the  gravity  vector  at  points  where  the 
measurements  have  been  made,  but  also  in  all  points  within  a  certain 
distance  from  the  measurements.  LSC  permits  the  estimation  of  the 
prediction  ei't'or  as  well. 


In  section  £  we  will  introduce  the  available  data,  and  in  section  3  we 
will  dcscr'ibe  the  use  of  LSC,  and  the  associated  process  of  estimating 
the  so-called  empirical  covariance  function.  The  results  obtained 
using  various  combinations  of  the  data  as  well  as  different  analytic 
models  for  the  empirical  covariance  function  are  described  in  section 
4. 


The  Southern  Ohio  test  data  set. 


Let  W  be  the  gr'avity  potential  and  U  a  cor-respond  i  ng  t'efet'erice 
( nor'i.ia  1 )  pot  ent  i  a  1 .  Then  the  anomalous  gr'avity  potential  is  the 
diffet'ence  T  =  W  -  U. 

The  long-wavelength  part  of  W  (  or  T)  is  described  down  to  a  1  deg. 
r'esolution  by  one  of  the  available  spherical  harmonic  expansions 
complete  to  degr'ee  and  ot'der  N=1S0  such  as  the  OSUSl  and  the  GPM£ 
coefficient  sets  (see  Rapp(lSSl),  Wenzel ( 1925) ) . 


Free-air  gravity  anomalies  (2lg)  have  been  published  by  Heiskanen  and 
Dot  i  la  (1956),  and  deflections  of  the  vertical  (i^,  by  Mueller  and 
Preuss  (1965).  They  are  in  spherical  approx i mat  ion  related  to  T 
through  the  following  equations  (Heiskanen  and  Moritz,  1967). 
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where  «  is  the  latitude,  X  the  longitude,  V  normal  gravity  and  r  the 
distance  from  the  origin,  which  in  spherical  approximation  is  obtained 
as  the  sum  of  the  Earth  mean  radius  and  the  height.  Torsion  balance 
data  are  published  by  Badekas ( 1967) ,  and  their  relations  to  T  are 
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where  (K,y,  z)  are  the  coordiriates  of  a  local  — level 
with  H  east,  y  north  and  z  up.  The  distribution  of 
data  IS  shown  in  Fig.  1. 


coot'dinate  systedi 
the  torsion  balaril 


Digital  topographic  heights  was  not  available.  However,  topographic 
effects  have  been  published  for  the  toi-sion  balance  in  the  same 
publication  as  the  observations  and  for  the  deflections  in  Badekas  and 
Mueller  (1963).  The  topographic  effects  on  the  gravity  data  may  be 
calculated  using  the  simple  Bouguer  plate  reduction,  since  the  terrain 
corrections  af'e  rather  small  due  to  the  smooth  topography. 


The  use  of  the  spherical  harmonic  coefficients  within  the  framework  of 
LSC  is  equivalent  to  subtracting  the  contribution  of  the  coefficients 
from  the  data.  Since  the  data  used  for  the  determination  of  the 
coefficients  primarily  are  free-air  gravity  anomalies,  the 
coefficients  also  reprensent  the  attraction  of  the  topographic  masses 
down  to  1  deg.  resolution.  Consequently  only  the  effect  of  a  Bouguer- 
plate  referring  to  the  mean  altitude  of  the  topography  (£65  m)  must  be 
subtracted  fro*m  the  gravity  anomalies.  Howevei',  this  creates  an 
inconsistency  between  the  topographic  corrections  for  the  gravity  data 
and  for'  the  other  data  types,  which  as  we  shall  see  later  will  result 
in  biases  in  the  predicted  results  using  topographic  reduced  data. 

This  problem  is  easily  avoided  in  practice,  if  a  detailled  topographic 
model  is  available. 


Table  1  shows  the  statistical  characteristics  of  (1)  the  "raw"  data, 

(6)  the  data  from  which  the  contribution  from  a  set  of  spherical 
harmonic  coefficients  (GPM£  to  maximal  degree  ISO)  have  been 
subtracted  and  (3)  the  data  <2)  minus  the  topographic  effects.  Also 
statistics  for  a  subset  of  the  gravity  data  located  in  an  area  close 
to  the  torsion  balance  data  is  given  in  the  table.  These  1£  gravity 
values  have  also  been  used  in  the  following  when  comparing  observed 
and  predicted  quantities. 

Note  the  striking  difference  in  the  mean  values.  In  fact,  the  area  is 
located  on  the  top  and  east  of  a  large  anomalous  mass.  Also  nc*te,  that 
the  deflections  of  the  vertical  have  been  transformed  from  NAD1927  to> 
the  preliminary  NAD33  using  datum  shifts  parameters  provided  by  the 


U.  S.  Natiorial  Geodetic  Survey.  The  gravity  values  have  been 
t ransforriied  ft'orn  the  P'otsdarn  system  to  GRSISSO  using  a  bias  correction 
of  13. 7  mgal.  Hence  a  small  bias  with  respect  to  IGSN  1971  may  still 
be  left. 

3.  Least-squares  collocation  and  covariance  functions. 

The  theory  and  practice  of  LSC  is  described  in  numerous  publications. 

P  recent  survey  is  given  in  Tscherning ( 19S5b) . 

Suppose  we  have  given  n  observat ions  related  to  the  anomalous  gravity 
field  through  n  linear  functionals  (eq.  (1)  —  (7))  , 

L  (T)  =  y  +  e,  i  =  l,...-n. 

Here  e  is  the  noise  or  error.  Then  let  C.  .  denote  the  covariance 

between  L.  (T)  and  L  .  <T)  and  C,  the  covariance  between  the  value  of  a 
1  J  Li 

linear  functional  L  applied  on  the  anomalous  potential,  L(T),  and 
L^(T).  We  then  obtain  an  estimate  L(f)  of  L(T)  as 

L(f)={C,  .}^{C..-*-D  .rMy.}  (S) 

Lx  i.J  ij  ‘  ' 

where  D. .  is  the  covariance  of  the  noise  associated  with  the  i’th  and 
xj  _ 

j’th  observations.  Let  C  denote  the  sum  of  the  C.  .  and  D.  .  matrices. 

IJ  XJ 

Then  an  estimate  of  the  mean  square  error  for  a  linear  functional  is 
obtained  by 

£  -  ^  j  --1 

<r  <L(T)  -  L(T))  =  C  O) 

where  is  the  vector  of  covariances  between  the  observations  and  the 
quantity  L(T)  and  is  the  variance  of  this  quantity. 

Similar  equations  are  obtained  if,  instead  of  the  observations  y,  we 
use  the  values  from  which  we  have  subtracted  the  effect  of  the 
topography  or  the  cont  t' i  but  i  on  from  a  spherical  harmonic  expansion.  On 
the  right  hand  side  (  of  eq.  <S) )  the  covariance  function  will  be  the 
one  associated  with  these  t'esidual  observations.  The  vector  of 
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observations  will  consist  of  these  residuals  arid  on  the  left  hand  side 
we  qet  the  rtjsidual  po-tential,  to  which  in  the  final  step  we  should 
add  back  the  potential  irnplied  by  the  topography  and  the  sphei'ical 
harfiionic  enpansion.  In  the  case  we  only  have  subtr'acted  the 
contr'i  but  ion  from  GPMc:,  we  will  associate  a  superscript  "s"  with  the 
relevant  quantity,  and  if  also  the  topographic  effects  have  been 
subtracted  we  will  use  the  superscript  "t",  i.e.  /Ig^. 


The  starting  point  for  the  use  of  LSC  in  a  local  area  is  the 
computation  and  analytic  modelling  of  the  empir'ical  covariance 
function.  Its  estimation  is  discussed  in  detail  in  Goad  et  al.  (19S4). 
Numer'ically  the  estimation  is  simply  done  as  the  computation  of  the 
mean  value  of  products  of  quantities  lying  within  the  same  interval  of 
spherical  distance.  The  size  of  this  interval  is  called  the  sampling 
interval.  For  vectot'  quantities,  like  cotnporients  of  deflections  of  the 
vertical,  the  components  of  the  horizontal  gravity  (disturbance) 
gradient  (  eq.  <4),  (5))  and  the  miKed  second  order  derivatives  (eq. 

(6),  (7)),  the  products  are  formed  between  products  of  "along-track" 

and  "cross-track"  components.  For  details,  see  Krarup  and  Tscherning  .i 
(1934).  " 


The  quality  of  the  estimation  depends  on  the  size  of  the  sampling 
error,  the  regularity  of  the  data  distribution  and  the  data  error.  The 
sire  of  the  area  should  also  be  so  large  that  the  mean  value  of  the 
quantities  (from  which  the  contribution  from  the  spherical  harmonic 
expansion,  here  GPM2,  have  been  subtracted),  ,  is  zero.  For  a 

spherical  harmonic  expansion  to  degree  ISO,  this  generally  will  be 
achieved  for  a  1.5-2  deg.  equal  angular  area,  since  mean  values  of  1 
deg.  blocks  have  been  used  for  the  computation  of  the  coefficients. 

For  the  approx  imately  1  deg.  block  sur'rounding  the  torsion  balance 
data  the  mean  value  of  the  available  227  Jg  -values  is  -3.3  mgal  (cf. 
Table  1).  Since  the  standard  deviation  is  19.3  mgal,  this  mean  value 
is  not  significantly  different  from  zero.  The  gravity  values  in  this 
area  were  therefore  used  for  the  estimation  of  the  covariance 
function,  both  using  and  Jg^  values  (cf.  Fig.  2  and  3). 
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When  estiiiiating  CT'OSs-covar  i  ance  functions  with  the  torsion  balance 
data,  we  have  the  p*'obleni  that  the  data  covers  a  much  smaller  ar'ea. 
Futhermore  (cf.  Table  1),  the  gravity  variation  in  the  0.75  deg  *  '.'.85 
deg.  area  is  significantly  different  fr-i-m  the  var-iation  in  the  ar'ca 
immediatly  surt'ound  i  ng  the  tor'sion  balance  data.  Here  the  variation  is 
much  lai'ger,  and  the  mean  values  are  non-zero.  We  decided  anyway  to 
estimate  empirical  auto-  and  cross— covariance  functions  between  the 
(GPMc:  and  GF'M2  +  terrain  reduced)  gravity  and  torsion  balance  data. 
Figur'es  4  to  1£  show  the  covar'iance  functions  for  the  GPMIE'  reduced 
data.  It  should  be  noted  that  the  mean  values  have  not  been  subti'actcd 
fr.;.fii  the  data. 

The  modelling  of  the  empirical  covariance  function  is  based  on  the 
fitting  of  an  analytic  expression  to  the  empii'ical  values.  We  used  a 
model  for'  the  gravity  anomaly  covariance  function,  developed  in 
Tscherning  and  Rapp  (1374), 


C<P,0)  = 


I 


A(i-1) 


.  ( I -2 ) (1 +4 )  ^  r  r 


^8  .i+2 
rr  J 


P. (cos^pq), 


( 10) 


with  the  two  free  parameter's  A  and  R  ,  In  this  expr'ession  t'*  is  the 
distance  of  the  point  Q  from  the  or'igin,  P^  the  Legendr'e— polynomial  of 
degree  i,  '♦'p.Q  the  sphor'ical  distance  between  P  and  Q  and  is  the 
r'adius  of  the  so-called  Bjer'harnmar  spher'e.  The  summation  star'ts  fr'om 
ISl,  since  we  have  subtracted  the  effect  of  GPM2  to  degr'ee  ISO.  This 
is  equivalent  to  pr'esupposing  that  the  coefficients  are  err'or  free,  a 
hypothesis  which  is  certainly  not  true.  However,  the  consequence  of 
adopting  the  hypothesis  will  mainly  be  seen  in  pr'edicted  geoid 
heights,  see  Arabe  1  os  ( 19S0)  .  Also  the  start  of  the  summation  fr'om 
i  =  lS0  made  it  r'easonable  to  intr'oduce  a  modification,  namely  the 
division  by  i-*-4  instead  of  a  division  by  i-*-24  in  eq.  (10)  as 
recommended  in  Ibid.  (1374).  This  modification  gives  some  computational 
s  a  V  i  ri  g  s . 

The  factor'  A  in  eq.  (3)  is  linear'ly  r'elated  to  the  value  at  spher'ical 
distance  zero,  the  var'iance.  The  value  of  R  is  mor'e  difficult  to 
estimate  due  to  the  uncertainty  in  the  estimation  of  the  covar'iances 
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and  the  weak  t'e  1  at  i  onsh  x  p  between  and  the  covar  i  ances ,  see 
Tscher-nintj  (19S5b).We  decided  to  try  to  find  the  best  one  (within  the 
class  of  functions  given  by  eq. (10)),  simply  by  carrying  out 
prediction  exper inients  with  diffei'ent  values  of  R_,  and  then  selecting 
the  one,  which  gave  the  best  agreement  between  observed  and  pr'edicted 
deflection  values.  The  results  are  sunuiiai'i zed  in  Fig.  13  .  The 
analytic  covariance  function  found  using  this  procedure 
had  a  depth  to  the  BJerhammar-sphere  of  3.0  km.  Its  graph  is  shown 
with  the  empit'ical  values  in  Fig.  3  -  13. 


It  should  be  noted  that  this  procedut'e  of  selecting  a  covariance 
function  rnc-del  ver'y  well  could  have  been  used  if  only  torsion  balance 
or  gradiometer  data  had  been  available,  ft  subset  of  the  observations 
could  have  been  used  as  test -quant  it ies  and  the  remaining  observations 
as  the  given  data,  The  covariance  model,  which  gave  the  best  results 
in  a  least  squares  sense  must  then  be  the  one  in  between  the  used  set 
of  models,  which  is  most  close  to  the  true  empirical  covariance 
f unct ion. 

4.  Prediction  tests  and  results. 


For  the  evaluation  of  the  predictions  (eq.  (S) )  as  well  as  the 
calculation  of  error  estimated  (eq.  (9))  the  FORTRAN  program  GEOCOL, 
(Tscherning,  19S5a)  was  used. 

In  order  to  check  the  consistency  of  the  gravity  vector  data,  we  first 
predicted  the  deflections  at  the  S  sites  shown  in  Fig.  1  from  the  337 
gravity  values  and  subsequently  the  13  gravity  values  in  the  vicinity 
of  the  torsion  balance  data  from  15  pairs  of  deflection  points  in  the 
same  area  as  covered  by  the  337  gravity  values.  The  results  are  given 
in  Table  3.  Note  how  surprisingly  well  the  gravity  anomalies  are 
determined  ft'om  the  deflections,  ftlso  note,  that  the  deflections  are 
pr'odicted  just  as  well  from  the  gravity  data  as  from  the  to-i'sion 
balance  data,  as  we  shall  see  later. 

We  tnen  wanted  to  predict  also  the  torsion  balance  data  from  gravity, 
torsion  balance  from  torsion  balance  data  and  gravity  and  deflections 
from  torsion  balance  data.  We  assigned  a  noise  standard  deviation  of 


3 


1.0  E. U.  to  the  torsion  balance  data,  equal  to  the  value  given  in 
Badekas  and  Mueller  < 19&fl,  p.  £371).  However  all  our  results  were 
"bad",  in  the  sense  that  the  ei'ror  estimates  computed  using  eq.  0> 
wer-e  much  smaller'  than  the  var'iance  of  obser'ved  minus  computed  values. 
We  then  spent  14  days  intensivly  c:'>ocking  GEOCOL,  but  found  no 
significant  error.  This  made  us  go  back  to  the  publication 
Badekas ( 1967 ) ,  which  contains  detailed  maps  showing  the  individual 
data  in  vector'-like  for'm.  This  revealed  that  neighbouring  values  could 
be  r'ather'  differ'ent,  signifying  a  much  lar'ger'  et'r'or'  level  than  we  had 
expected.  Olso  the  empir'ically  estimated  auto— covar'i ance  functions 
indicated  that  a  r'ather'  lar'ge  noise  was  pr'esent.  This  can  be  seen  fr'om 
Figur'es  3  to  11,  wher'e  the  var'iance  is  much  lar'ger'  (ar'ound  100  E.  U.  ^  ) 
than  the  value  found  by  smoothly  extrapolating  fr'om  the  fir'st  £  —  3 
empir'ical  values  to  the  value  at  spherical  distance  zer'o  (the 
var'  i  ance )  . 


Then  changing  the  noise  var'iance  to  100  E.  U.  we  suddenly  obtained 
consistent  r'esults.  This  value  is  probably  too  pessimistic,  but  not 
impossible  ( I .  I .  Mue  1  1  er,  per'scnal  cornmunicat  ion,  1936).  We  then  tr'ied 
using  changing  noise  levels.  The  r'esults  ar'e  illustrated  in  Fig.  14  - 
16.  In  fact  values  of  the  noi se^etween  50  and  150  £.0.“^  all  give 
rather  similar  prediction  r'esults,  but  100  E.  U.  gave  the  best 
r'esults.  We  then  used  this  value  in  the  following  computations. 


Having  arr'ived  so  far',  we  decided  to  make  a  series  of  computational 
experiments,  which  should  illustr'ate  the  dependence  of  the  pr'ediction 
result  on 


(1)  the  model  selected  for'  the  covar'i  ance  function 
(£)  the  use  or'  non-use  of  the  topography 

(3)  the  data  type  and  che  combination  of  var'ious  data  types 

(4)  the  data  density. 


The  dependence  on  the  covar'iance  function  as  mentioned  alr'eady  has 
been  used  to  select  a  "best"  analytic  model.  The  dependence  on  the 
model  (the  depth  t.:<  the  B jer'hammar'— spher'e)  is  illustr'ated  in  Fig.  i;-;. 

The  use  of  t  opogr'aph  i  c  data  r'esults  in  a  slight  smoothing  of  the 
gr'avity  field  data  (cf.  Table  1),  and  we  should  expect  a  similar' 
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i fiiprovewcprit  of  the  pr'ediction  t'eiiult 
not  ver'y  much  is  gained,  but  this  is 
density.  Ttie  topographic  i nf orwat  ion 
the  data. 


s.  However,  as  seen  from  Table  3 
probably  due  to  the  high  data- 
is  a-priori  well  represented  by 


In  order  to  see  the  influence  of  the  data  density  we  used  either  all 
data  points  or  one  value  from  each  £’ ♦£’  cell.  Using  all  data  gave 
only  a  small  improvement  as  seen  from  Table  5.  More  interesting  is  the 
clear  difference  in  contribution  from  the  various  torsion  balance 
components  to  the  components  of  the  gravity  vector.  We  maybe  should 
expect  the  significant  contribution  o-f  the  horizontal  gradient  of  the 
gravity  disturbance  to  the  deflection  components  from  geometrical 
reasons  and  from  the  results  presented  in  Table  £,  where  gravity  data 
wet'e  used  to  predict  deflections.  The  mixed  second  order  derivatives 
contributes  little,  but  anyway  something  of  significance. 

It  is  clea?'  that  we  can  not  estimate  the  absolute  values  of  the 
gravity  vector  co<mponents  from  torsion  balance  data  alone.  This  is 
also  confirmed,  since  we  get  relatively  large  biases  (cf.  Table  3) .  ^ 

However,  gravity  vector  differences  are  well  determined  as  seen  from 
the  small  standard  deviations  in  the  table.  On  the  other  hand,  we  can 
easily  add  observed  gravity  vector  infoi'mation  to  the  input  data  in 
LSC.  We  selected  one  of  the  deflection  points,  which  coincide  with  one 
of  the  points  of  the  torsion  balance  network,  but  we  had  to  select  a 
gravity  station  at  some  distance  (  1  km  )  from  the  network.  The  result 
is  given  in  Table  4.  We  see  that  most  of  the  biases  have  decreased, 
but  the  bias  in  the  topographically  reduced  /ig  has  increased.  The 
reason  for  this  may  as  mentioned  be  due  to  the  inconsistency  of  the 
topographic  reductions. 

Contour  plots  showing  the  predicted  gravity  anomaly  and  deflection 
components  with  the  estimated  errors  of  prediction  (eq.  O) )  are  shown 
in  Fig.  17  £J^ 


.  Cone  1  US  ion. 
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We  have  seen  how  LSC  is  suitable  for  gravity  vector  estimation  in  a 
local  area  from  qrad i omet er- 1 i ke  data.  If  the  method  should  be  used 
economically  for  data  distributed  over  large  areas,  then  the  area  must 
be  split  in  sub-areas  of  reasonable  size. 

The  results  obtained  in  the  best  cases  (  O.A"  for  deflections  and  1.6 
mgal  for  gt'avity  anomalies)  is  very  encouraging  since  the  qradiometer 
is  supposed  t  .-  have  a  noise  level  below  1  E.  U.  However,  as  we  have 
seen  it  is  important  to  use  a  reasonable  signal  to  noise  ratio. 

The  proper  removal  of  biases,  using  observed  components  or  the  gravity 
vector  requires  some  futher'  i  nvest  i  gat  ions.  Here  a  data  set  covering  a 
larger  area  would  be  suitable,  permitting  the  proper  estimation  of  all 
cross-covar i ance  funct ions. 
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Figure  captions. 


Fig. 


Fig. 


Fig. 


Fig. 


Fig. 


Fig. 


Fig. 


Fig. 


Fig. 


Fig. 


Fig. 


1.  The  distribution  of  the  torsion  balance  (-•->,  gravity  data  (•> 
and  deflections  of  the  vertical  (o) ,  with  the  free-air 
gravity  anoriialy  field  contoured  at  2.  S  mgal  interval. 

2.  Eriipirical  (solid  line)  and  synthetic  (broken  line)  covariance 
functions  for  the  gravity  anomalies  in  Southern  Ohio.  GPM2 
t'eriioved  from  the  original  data.  Units  mgal‘". 


2.  Empirical  and  synthetic  covariance  functions  for  the  gravity 
anomalies  in  Southern  Ohio.  GPI<12  and  topographic  (simple 

Bouguer  plate)  effect  removed  from  the  original  data.  Units 

2 

mgal  . 


4.  Empirical  and  synthetic  crosa-covar iance  functions  between 
dg  and  T  .  Units  Aigal*E.  U. 

M  S 

5.  Empirical  and  synthetic  cross-covariance  functions  between 
dg  and  T^.  Units  mgal^E.U. 


&.  Empirical  and  synthetic  cross-covariance  functions  between 

2 

2*T  and  T  Units  E.  U.  . 
xy  A 


7. 


Empirical  and  synthetic  cross-covariance  functions  between 

Aq  and  EV'T  .  Units  mgal-o-E.  U. 
xy 


S.  Empirical  and  synthetic  auto-covariance  functions  for  the 

2 

T^_  component.  Units  E. U.  . 


9. 


Empirical  and 
T^_  component. 


synthetic  auto-covariance 
Units  E. U. 


f unct ions 


for  the 


10.  Empirical  and  synthetic  auto-covariance  functions  for  the 

2*T  component.  Units  E.  U.  ■"  . 
xy 

11.  Empirical  and  synthetic  auto-covariance  functions  for  the 

component.  Units  E.  U. 

A 
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Fig.  li.  Ewpi  t'iccx  I  arid  synthetic  cross-covariance  functions  between 

T  and  T  Units  E.  U. 

Ky  A 

Fig.  13.  Dependence  of  the  prediction  results  (expressed  in  terms 
of  standard  deviations  of  differences  observed  minus 
predicted  quantities)  on  the  covariance  function  (  the  depth 
to  the  Eijerhammar  sphere),  h  =  ,  A  -  t).  Topographic 

reduced  data  shown. 

Fig.  14.  The  effect  of  changing  of  the  noise  level  on  the  prediction 
results  of  g»'avity  anomalies  from  torsion  balance  data, 
(expressed  in  terms  of  standard  deviation  of  differences 
observed  minus  predicted).  A  =  unreduced,  x  =  topographic 
reduced  data. 

Fig.  15.  The  same  effect  as  in  figure  14,  for  the  tq  component 
of  the  deflection  of  the  vertical. 

Fig.  16.  The  same  effect  as  in  figure  14,  for  the  g  component 
of  the  deflection  of  the  vertical. 

Fig.  17.  Prediction  of  gravity  anomalies  from  torsion  balance  data. 
(Contour  interval  is  5  mgal). 

Fig.  IS.  Error  of  prediction  of  gravity  anomalies  from  torsion  balance 
data.  (Contour  line  interval  is  O.ds  mgal). 

Fig.  13.  Prediction  of  ig  from  torsion  balance  data.  (Contour  interval 
is  O. 5  arcsec) . 

Fig.  SO.  Error  of  prediction  of  from  torsion  balance  data.  (Contour 

f 

interval  is  0. iS  arcsec). 

Fig.  El.  Prediction  of  g  from  torsion  balance  data  using  eq. (3). 
(Contour  interval  is  0.5  arcsec). 

ES.  Error  of  prediction  of  g  from  torsion  balance  data.  (Contour 

t 

interval  is  arcsec). 


Fig. 
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Table  1.  Statistical  character isc ics  of  the  raw  and  reduced  data 
(  227  gravity  points,  1£  gravity  points,  3  pairs  of  deflections 
of  the  vertical  and  233  points  of  torsion  balance  observations) 


Raw  data 

Dat  a-GPM2 

Data-GPM2- 

topography 

St  andard 

mean  deviation 

st  andard 

mean  deviation 

standard 

mean  dev  i  at  i  or^ 

JJg 

(mgal)  * 

-1. 24 

19.  31 

16.  74 

Jg 

( rnga  1 )  ** 

22.  64 

£7.  66 

16.  31 

25.  93 

16.  24 

22.  41 

(arcsec) 

0.  92 

1. 44 

0.  93 

1. 49 

-0.  12 

1 . 70 

n 

( arcsec) 

4.  41 

4.  72 

4.  34 

4.  74 

5.  39 

4  7  ’ 

T 

X  z 

(E.  U.  ) 

~0.  56 

16.  £0 

~0.  52 

16.  24 

-1. 07 

14i9 

T 

(E.  U.  ) 

-3.  69 

20.  59 

-3.  73 

20.  63 

-7.  39 

19.  40 

(E.  U.  ) 

22.  62 

-0.  37 

22.  65 

-1.  12 

20.  10 

2*T 

xy 

(E.  U.  ) 

2.  33 

19.  70 

•  <1^3 

19.  74 

2.  43 

1 7.  53 

*  £27  points  *♦  12  points 
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Table  li.  P'redictiori  results 


of  the  two  cc<riiporierit 5 

tion  at  a  points  from 

po i nt  s 

of  the  deflec-  of 

£37  gt'avity  of 

1£  gravity  values  from 

deflection  values 

in  15  points 

n 

jUg 

obs. 

-  pred- 

obs.  -  pred. 

l  nput 

st  andard 

st  andard 

st  andard 

u  a  t  a 

mean  deviation 

mean  dev i at i on 

mean  deviation 

(arcsec) 

<  arcsec) 

(mgal ) 

'.irireduced  0.03  0.71 

topogr. 

-0.  1 7  0.  43 

5. £0  4. 3£ 

r'ed  uced 

-l.£7  0.  A7 

1.5£  0.47 

IS 

Table  3.  Prediction  of  gravity  anc-Malies  and  deflections  of  the  vertical 
frora  torsion  balance  observations  using  least  squares  collocation 


iJg 

(ftigal  ) 

(arcsec) 

0 

Input 

« 

standard 

standard 

standard 

components 

mean 

deviat ion 

mean  deviation 

mean  deviation 

remarU 

T 

10.  A3 

22.  AS 

-0.  27 

0.  73 

2.  70 

3.  69 

a 

T 

A.  97 

A.  37 

-0.  1.1* 

0.  76 

0.  AO 

0.  72 

- 

13.  7A 

20.  03 

-0.  95 

1. 70 

2.  51 

3.  17 

— 

2*T 

17.  23 

26.53 

0.  lA 

1.21 

2.  SA 

A.  00 

- 

xy 

T  .  T 

0.  79 

3-  AO 

-0.31 

0.  90 

0.  A  1 

0.  70 

- 

HZ’  yz 

T  2*T 

16.  35 

1A.9A 

-O.  A  5 

O 

Id 

• 

o 

1.71 

2.  92 

— 

A  xy 

ALL 

5.03 

3.  A1 

-0.  OA 

0.  7A 

0.  IS 

0.  56 

— 

T 

H  "Z 

10.  07 

22.66 

-0.  20 

0.  91 

2.  62 

3.  3S 

b 

T 

A.  A3 

5.  AA 

-0.  12 

0.  50 

0.  55 

0,  A7 

- 

ys 

12.  63 

19.26 

-1. 07 

1.  57 

2.  57 

2.91 

— 

2*T 

17.  59 

2A.  38 

0.  23 

1.21 

2.  72 

3.  93 

- 

xy 

T  .  T 

0.  OA 

3.  27 

-0.  20 

1.03 

0.  63 

u.  5c: 

— 

X 

tl 

u 

T  ..  2*T 

15.55 

3.23 

-0.  2A 

1.07 

1. 37 

2.71 

— 

1 

X 

ALL 

-0.27 

2.  35 

0.  AS 

0.  70 

0.  7A 

0.  31 

— 

T 

H  Z 

10.  39 

19.  32 

-1.  SA 

0.  AA 

3.  7A 

3.  &A 

c 

T 

5.  lA 

5.03 

-1. 32 

1.07 

1.  S3 

0.87 

- 

y= 

lA.  13 

16.  13 

-1.79 

1. 50 

3.57 

3.  2A 

— 

2*T 

xy 

16.79 

2A.03 

-1.05 

1.22 

3.  96 

3.  99 

- 

T  .  T 

1.51 

2.  7A 

-1. 6A 

0.61 

i.  37 

0.  77 

— 

N 

M 

X 

T  ..  2*T 

16.  AA 

13.21 

-1.  AA 

1.  AS 

3.  03 

3.09 

— 

Hy 

ALL 

3.  36 

2.  A7 

-1. 50 

0.  53 

1. 75 

0.  71 

T 

xy 

10.  6A 

19.  37 

-1. 35 

0.  63 

3.  70 

3.  SO 

d 

T 

y= 

A.  15 

5.  37 

-1. 2A 

0.  79 

1. 3A 

0.  60 

13.  9A 

15.  29 

-1. 39 

1 . 39 

3.  65 

3.  1 1 

— 

2*T 

xy 

17.  00 

21. 7A 

-0.  39 

1.27 

3.  83 

3.  9A 

- 

T  ,  T 

M2  y2 

0.  A2 

2.  A1 

-1.  A1 

0.30 

1. 95 

0.  56 

2*T 

A  xy 

15.91 

7.  21 

-1.  IS 

1.  lA 

2.  39 

2.96 

ALL 

-0.  93 

1.71 

-1.  17 

0.  AA 

1. 85 

0.  A3 

a  ;  without  topographic  reduction,  one  point  per  2’*  2’  cell  used 
without  topographic  reduction,  all  points  used 
with  topographic  reduction,  one  point  per  2’*  3’  cell  used 
with  topogr'aphic  reduction,  all  points  used 
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Table  4.  Prediction  results  using  all  torsion  balance  data  and 
in  addition  one  gravity  point  and  one  deflection  point 

p - 


^g 

(mg a  1 ) 


(arcsec) 


n 


i  nput 

standard 

standard 

st  andard 

data 

mean 

dev i at  ion 

mean  dev i at  ion 

mean  deviation 

un reduced 

0.  05 

£.  15 

0.  £6  0.  6S 

0.  0£  0.30 

topogr. 

reduced 

« 

O  I 

<  1 

OJ  1 

-0. 06  0. 49 

1 

1 

1 

u)  1 

•  1 

O  1 

1 

1 

I 

1 

UT  1 

UT  1 

•  1 

O  1 

1 

1 

1 

1 

1 

1 

-83 


£0 


'SI  '010.1 


Fig.  8 


Fig.  9 


'SI  '0U5.I 


Fig.  10 
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ESTIMATION  OF  GRAVITY  VECTOR  COMPONENTS 


FROM  BELL  GRADIOMETER  AND  AUXILIARY  DATA 
UNDER  GONSIDEPATION  OF  TOPOGRAPHY  AND 
ASSOaATED  ANALYTICAL  UPWARD  CONTINUATION  ASPECTS 


by  . 

Dr.  Hans  Baussus  von  Luetzov 
Department  of  the  Army 
U.S.  Army  Engineer  Topographic  Laboratories 
Fort  Bel  voir,  VA  22060-5546 


ABSTRACT 

Following  an  introduction,  the  paper  briefly  discusses  gravity  gradiometer 
applications.  It  outlines  the  estimation  of  first  order  derivatives  of  the 
anomalous  gravity  potential  from  Bell  gravity  gradiometer  and  auxiliary  data  in 
the  context  of  a  Welner-Kolmogorov  optimization  scheme  under  consideration  of 
computable  "topographic  noise",  accomplished  on  the  basis  of  the  Pellinen-Moritz 
solution  of  the  boundary  value  problem  of  physical  geodesy.  The  paper  also 
addresses  four  different  methods  of  analytical  upward  continuation  of  first 
order  derivatives  of  the  anomalous  gravity  potential  under  identification  of  a 
finite  difference  method  using  Laplace's  equation  as  the  most  economical  and 
efficient  one.  Relevant  conclusions  are  then  presented. 


1 .  INTRODUCTION 


The  estimation  of  gravity  vector  components  from  Bell  gravity  gradiometer 
and  auxiliary  data  has  been  discussed  by  a  number  of  authors,  including  Jordan 
(1982),  White  and  Golfjstein  (1984),  Center,  Jordan  and  Peacock  (1985),  and 
Baussus  von  Luetzow  (1985).  Apart  from  gravity  gradiometer  self  noise  power 
spectra  or  corresponding  covariance  functions,  the  collocation  theory  of 
physical  geodesy  was  employed  to  derive  spectra  and  cross-spectra  or  equivalent 
covariance  functions  involving  first  and  second  order  derivatives,  jof  the  anomalous 
gravity  potential.  The  basic  covariance  functions,  pertaining  either  to  a 
randomized  anomalo'us” 'gravity'  potent'ial  or  to  a  randomized  gravity  anomaly, 
accomplished  by  the  subtraction  of  low  frequency  harmonic  representations  of 
anomalous  potential  or  gravity  anomaly  data,  presuppose  the  existence  of  homogeneity 
and  isotropy.  The  application  of  collocation  theory  does  not  face  difficulties 
in  the  case  of  quasi-flat  terrain  and  is  a  prerequisite  to 
arrive  at  small  estimation  errors  for  gravity  vector  components 
in  the  context  of  an  airborne  survey  over  a  square  area  of  300  x  300  km*,  covered 
with  parallel  traverses  and  cross  traverses  5  km  apart,  conducted  at  an  elevation 
of  600  m  at  a  speed  of  about  360  km  hr  ,  and  employing  only  a  limited  amount  of 
measured  gravity  vector  component  data.  It  is,  however,  well  known  that  the 
assumption  of  homogeneity  and  isotropy  does  not  satisfactorily  hold  in  the  presence 
of  moderate  to  strong  mountainous  terrain.  Some  efforts  have  been  made  to  apply 
heuristic  topographic  corrections  to  gravity  gradiometer  measurements.  Still, 
the  problem  remains  to  provide  for  a  Wiener-Kolmogorov-type  estimation  process 
free  of  computable  topographic  noise  so  that  the  assumption  of  homogeneity  and 
isotropy  is  reasonably  fulfilled  with  respect  to  modified,  signal-type  measurements 
disturbed  by  gradiometer  self  noise  only.  The  solution  to  this  problem  is 
possible  under  application  of  approaches  to  the  geodetic  boundary  value  problem 
and  the  correlated  interpolation  of  gravity  anomalies  and  deflections  of  the 
vertical  in  mountainous  terrain  by  Pellinen  (1969),  Moritz  (1969),  and  Baussus 
von  Luetzow  (1971,  1981).  Section  3  as  the  main  section  of  this  paper  addresses 
the  solution  under  consideration  of  the  estimation  process  and  the  computation 
of  "topographic  noise".  The  preceding  section  addresses  in  a  compact  way 
gravity  gradiometer  applications  to  place  the  determination  of  gravity  vector 
components  over  relatively  large  areas  and  the  analytical  upward  continuation 
of  vertical  and  horizontal  derivatives  of  the  anomalous  gravity  potential  in 
the  proper  perspective.  The  relatively  short  sections  4-7  discuss  different 
methods  of  analytical  upward  continuation,  followed  by  several  conclusions. 

2.  GRAVITY  GRADIOMETER  APPLICATIONS 

■  "i  ~  •  - 

Geodetic  research  and  the  development  of  associated  technology  has  been 
considerably  stimulated  by  military  requirements  and  actual  or  potential 
military  applications  with  critical  funding  provided  by  the  U.S.  Government. 

The  result  thereof  was  an  increase  of  scientific  knowledge,  new  instrumentation 
and  techniques,  and  technology  transfer  to  the  civilian  sector,  particularly  in 
navigation  and  surveying. 

The  main  applications  of  gravity  gradiometers  follow: 

®  Establishment  of  the  spatial  gravity  field  in  ICBM^  launch  areas  (most  important 
application).  This  includes  gravity  vector  components  on  the  ground  for 
calibration,  local  to  global  coordinate  transformations,  geodetic  azimuth 
determination,  and  gravity  programmed  inertial  positioning  (beneficial  to 
Small  ICBM  or  Midgetman) .  An  accurate  spatial  gravity  field,  extending  to 
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Intercontinental  Ballistic  Missile 


2 


a  height  of  200  Icm,  provides  also  for  comparisons/calibration  for  the 
Geopotential  Research  Mission  (GRM)  and  a  Superconducting  Gravity  Gradiometer 
Mission  (SGGM).  Instrument:  Bell  Gravity  Gradiometer. 

2 

®  Accuracy  enhancement  of  SLBM's  .  Instrument:  Bell  Gravity  Gradiometer. 

®  Improvement  of  global  gravity  field  for  geophysical  applications  (GRM). 
Instrument:  Superconducting  Gravity  Gradiometer. 

Byproducts:  ICBM  target  error  reduction  (not  critical),  cruise  missile 
navigation,  improved  satellite  epheffl«.i:is. 

“  High  accuracy  space  navigation.  Instrument:  Superconducting  Gravity 
Gradiometer  (six  axis  version). 


*’  Gravity  vector  densif ication  over  large  areas  for  gravity  programmed 

inertial  navigation  and  for  geodetic  network  adjustments.  Instrument:  Bell 
Gravity  Gradiometer. 

”  Subterraneous  mass  detection  (cavaties,  oil,  water).  Instrument:  Super¬ 
conducting  Gravity  Gradiometer  (vertical  axis  version). 


”  Test  of  Newton's  square  law  and  tests  of  Einstein's  theory  of  relativity. 
Instrument:  Superconducting  Gravity  Gradiometer. 

3.  DETERMINATION  OF  DERIVATIVES  OF  THE  ANOMALOUS  GRAVITY  POTENTIAL  OVER 
MOUNTAINOUS  TERRAIN  USING  BELL  GRAVITY  GRADIOMETER  MEASUREMENTS  AND 
AUXILIARY  DATA 


A.  Notations  and  Relations  Used 


/ 

A.  2 

7^21 

* 


/. 


anomalous  gravity  potential 

derivatives  taken  along  the  local  horizon  in  a  northern  and 
eastern  direction 


n  s  Bell  gravity  gradiometer  instrument  noise  with  white  noise 

component  /J,  and  red  noise  component 


iT, 

T, 


3T 


signal  variable  indicator 

topographic  noise  component  of 

topographic  noise  component  of 

lower  order  spherical  harmonics 
subtraction 

lower  order  spherical  harmonics 
subtraction 

lower  order  spherical  harmonics 
for  subtraction 


representation  of 
representation  of 
representation  of 


used  for 
Tjfy  used  for 
used 


Sea  Launched  Ballistic  Missile 
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lower  order  spherical  harmonics  representation  of 
used  for  subtraction 


t  -  T,-ST, 

i  =  i-'t  = 

i  -t-t  =  4 - 


weight  factor  with  measurement  index  i  and  computation 
index  /(  ,  applicable  to  7^  in  Uiener-Kolmogorov  estimation 


weight  factor  with  measurement  index  J  and  Ci  mputation 
index  K  ,  applicable  to  in  Wiener-Kolmogorov  estimation 

earth's  mean  radius 


6- 

r 


'p 


Ah 


unit  sphere  (full  solid  angle) 

length  of  radius  vector  from  earth's  center  to  a  moving  point 
on  earth's  surface 

length  of  radius  vector  from  earth's  center  to  a  fixed  point 
on  earth's  surface 

length  of  radius  vector  from  earth's  center  to  a  fixed  point 
on  a  level  surface  above  earth's  surface 

angle  between  radius  vectors  r  and  K  oz 

f  A 

( -  ^.T-r  CASi')^ 

A  ^ 

azimuth  angle  counted  clockwise  from  north 

S (y  V*' TV  generalized  Stokes  Function 

'A'  / 

elevation  of  terrain  referring  to  a  movable  point  on  earth's 
surface 

elevation  of  terrain  referring  to  a  fixed  point  on  earth's  surface 
elevation  of  level  surface  above  terrain 

4 


4 


/  =  2 
“  JL 


planar  approximation  to 


e 


^3 

c 

^dp  ‘^3^^ 
b 

/ 


gravity  anomaly  on  earth's  surface 
topographic  correction  to  gravity  anomaly 
Faye  anomaly 

Bouguer  gradient 
gravitational  constant 
standard  density  ^ 


ST  =  jj/nA-L^  Js’ 

potential  of  the  topographic  masses 


with  V  as  an  arbitrary  variable 


<P 

o  -  f  -  <P 


geographic  latitude 
geographic  longitude 


B.  Estimation  Process 


A  Wiener-Kolmogorov  estimation  can,  under  consideration  of  ^ 
only,  be  characterized  by 


and 
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The  derivatives  ,  also  called  truth  data,  are,  of  course^ also  associated 
with  measurement  e^rrors.  With  respect  to  deflections  of  the  vertical, 
proportional  to  and  Ty  ,  efforts  are  being  made  to  determine  them 

astrogeodetically  with  an  error  of  0.1  arcsec  rms .  Gravity  anomalies 
can  be  measured  with  a  relatively  greater  accuracy.  In  order 
to  compute  the  vertical  derivative  ^  ^  .  high  degree  spherical  harmonic 

expansions  presented  by  Rapp  (1987)  or  fhe  GRIM3-L1  Model  described  byReigber, 
Balmino,  Muller,  Basch  and  Moynot  (1985)  may  be  employed.  These  models  are 
also  useful. for  the  subtraction  of  lower  harmonics  components  from  measured 
first  and  second  order  derivatives  of  7”  . 

Regarding  eq.  (1),  it  is  assumed  that  all  measured  quantities  refer  to  a 
level  surface.  In  the  case  of  quasi-flat  terrain,  //(ff  -estimates  have  been 
made  from  7^^  on  the  earth's  surface  and  from  airborne  Tj^k  -measurements, 
employing  spatial  covariance  functions.  This  is  not  possible  in  the  presence 
of  pronounced  topography.  In  this  case,  it  is 


where  ■  is  a  representative  gradient  obtained  from  gravity  gradiometer 

measurements.  There  is  no  significant  degradation  in  accuracy  if  eq.  (2)  is 
used . 

For  simplicity,  eq.  (1)  presupposes  the  prior  application  of  quasi-systematic 
corrections  obtained  from  lower  harmonic  expansions.  In  the  presence  of 
pronounced  mountainous  terrain,  eq.  (1)  is  reformulated  as 


C.  Computation  of  "Topographic  Noise" 

The  essential  task  is  the  computation  of  topographic  noise  components 
of  the  derivatives  ,  7^  ,  ^  •  In^this  respect,  the  solution  nomenclature 
employed  by  Moritz  (1969)  is  applied.  In  order  to  obtain  solutions  at 
elevation  instead  of  solutions  at  elevation  h  hp  ,  Stokes  generalized 

function  replaces  Stokes  function  the  solutions  where  applicable.^ 

First,  the  solution  for  T  at  a  point  is  written  as 


/ 


^  Moritz  ( 1969) ,  p.  31 . 

is  explicitly  formulated  by  Moritz  (1966).  See  p.  A9 . 


6 


The  solution  components  ace 


4  = 


A 


(6) 


6* 


/, 


4  *h 


4 


C7) 


Tf  represents  the  "signal"  term  of  /  since  its  structure  is  compatible 
with  the  existence  of  a  reproducible  kernel  in  the  sense  of  Krarup,  a 
prerequisite  for  the  application  of  the  collocation  method  of  physical  geodesy. 


/j  and  are  clearly  "noise"  terms  which  do  not  permit  the  establishment 
of  a  statistical  estimation  structure.  7^  appears  to  be  a  more  slowly  varying 
function.  Moritx  (1969)  states  that  the  te_rm  ST  may  be  neglected  in  the  context 
of  a  planar  approximation.^  Accordingly.  appears  to  be  mainly  responsible 
for  fluctuations  of  a  non-stationary  nature. 


An  improvement  of  the  above  analysis  is  possible  by  the  replacement  of 
by  the  isostatic  anomaly  ,  provided  that  isostatic  conditions  prevail. 

Neglecting  a  second  order  improvasott  involving  consideration  of  ^  -terms,  the 
additional  "topographic  noise"  can  be  written  as 


Without  consideration  of  the  term  7^  ,  the  "topographic  noise"  terms 
involving  vertical  and  horizontal  derivatives  of  /j^  and  ^  can  be  formulated  as 
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D.  Averaged  Second  Order  Derivatives 


Because  of  the  application  of  a  moving  average  of  gravity  gradiometer 
measurements  over  a  time  interval  4/  of  about  10  seconds,  it  is  not  necessary 
to  compute  second  order  derivatives  vith  respect  ,  Z  ot  and  7^  .  As 

an  example,  the  averaged  second  order  derivative  of  ^  vith  respect  to  y  is 


In  eq,  (13)  it  is  ^  ^  ^  ^/ )  ^ 

the  horizontal  aircraft  speed. 


Oi) 


"2^  being 


4 .  STANDARD  ANALYTICAL  UPWARD  COHTINUATION 


Assuming  that  the  three  first  order  derivatives  of  the  anomalous  gravity 
potential  have  been  estimated  from  gradiometric  surveys  at  rectangular  grid 
points  of  a  level  surface  at  elevation  -const.,  Poisson's  integral  formula 


C/¥) 


s 


can  be  used  for  the  analytical  upward  continuation  of  7* or  of  its  derivati 


ives 


'M.  .  According  to  Heiskanen  and  Moritz  (1967),  for  elevations 

smaller ^han  250  km,  the  planar  form  of  eq.  (14)  is  sufficiently  accurate.  For 
elevations  up  to  150  km,  which  apply  in  ICBM  launch  areas,  the  horizontal 
integration  distance  has  to  be  ten  times  the  elevation  to  assure  good  accuracy, 
i.e.,  a  very  large  information  base  is  required.  Additionally,  rectrangular 
grid  data  has  to  be  interpolated  to  conform  with  circular  zone  data  to  be 
employed  in  the  computation. 

5.  IMMEDIATE  SPATIAL  WitMER  KOtMOGOROV  ESTIMATION 


The  three  derivatives  of  the  anomalous  gravity  potential  may  be  estimated 
directly  for  grid  points  at  higher  elevations  in  the  context  of  a  Wiener-Ko Imogorov 
scheme,  using  spatial  covariance  functions  or  equivalent  power  spectra.  For 
higher  elevations,  data  derived  from  a  greater  number  of  adjacent  surveys  would 
be  required,  and  matrix  inversions  would  become  very  complicated.  Further 
complications  would  arise  in  the  presence  of  pronounced  topography  because 
of  the  non-applicability  of  effective  collocation  methods. 

6.  SPATIAL  COLLOCATION  UPWARD  CONTINUATION  OF  ,  Tx 

If  ^  are  first  estimated  at  grid  points  of  a  level  surface  with 

elevation  ,  a  collocation  upward  continuation  would  be  a  replacement  of 
the  Poisson  integral* formula  approach  by  statistical  methods.  Again,  a  large 
information  base  would  be  required,  inversions  of  matrices  associated  with 
higher  elevations  would  become  complicated,  and  complications  would  exist  in 
the  presence  of  pronounced  topography.  It  would  be  possi^e  though  to  apply 
the  collocation  extrapolation  only  to  signal  quantities  ^  ^  and 

to  compute  the  "topographic  noise"  components  separately  for  higher  elevations 
by  means  of  eqs .  (9)  -  (12). 


7;  NUMERICAL  UPWARD  CONTINUATION  USING  LAPLACE’S  DIFFERENTIAL  EQUATION 


The  upward  continuation  methods  discussed  in  sections  4-6  are  characterized 
by  the  use  of  a  specific  derivative  of  the  anomalous  gravity  potential,  given 
at  a  level  surface,  by  the  requirement  of  a  large  information  base  for 
extrapolations  to  higher  elevations,  by  associated  large  inversion  matrices 
in  collocation  approaches,  and  by  the  inadequacy  of  these  in  the  presence  of 
mountainous  terrain  in  the  absence  of  a  separation  into  signal  estimation 
and  "topographic  noise"  computation.  Fortunately,  the  availability  of  three 
vector  components  at  a  level  surface  makes  it  possible  to  use  information  at 
two  adjacent  level  surfaces  as  a  prerequisite  for  applying  Laplace's  equation. 
Hereby,  the  information  base  can  be  reduced  significantly,  and  inversions 
of  matrices  are  not  needed. 


In' Cartesian  coordinates,  the  availability  ot  ^ 

^  s at  the  levelZ^  makes  it  possible  to  compute 
the  level  -  2^  t  AZ  since 

1  2/  .  I  U 

iz  Jx  22  55 


(/fc) 
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The  derivative  can,  of  course,  be  determined  separately  through  use  of 

gradiometer  information,  and  a  representative  value  can  be  employed  for 
analytical  upward  continuation. 

If  a  variable  -f  satisfying  Laplace's  equation  is  given  at  two  adjacent 
levels  separated  by  •  Laplace's  equation 


l/^) 


permits  application  of  the  numerical  algorithm 

/.  '  -/»  - 

In  more  appropriate  spherical  coordinates  T  ,  0  ,  \  ,  eq.  (17)  has  to 
be  replaced  by 


where  the  derivatives  in  the  bracketed  term  are  to  be  evaluated  by  finite 
differences.  Higher  order  finite  difference  algorithms  may  be  applied  if 
considered  advantageous. 

It  is  possible  to  use  eq,  (18)  only  for  "signal"  components  and  to 
compute  the  "topographic  noise"  effects  separately,  using  grid  lengths  of  1  km 
instead  of  5  km. 

8.  CONCLUSIONS 

Wiener-Kolmogorov-type  estimation  of  gravity  vector  components  or  of 
related  first  order  derivatives  of  the  anomalous  gravity  potential  from  Bell 
gravity  gradicmet^r  and  auxiliary  data  in  the  context  of  an  airborne  area 
survey  can  also  be  accomplished  in  the  presence  of  pronounced  mountainous 
terrain.  A  specific  solution  of  the  boundary  value  problem  of  physical 
geodesy  permits  the  determination  of  "signal"  and  "topographic  noise" 
solution  terms.  "Topographic  noise"  terms  can  be  computed  for  each  first 
derivative  of  the  anomalous  gravity  potential,  and  their  consideration  permits 
the  application  of  algorithms,  including  numerical  weight  factors,  valid  in  the 
case  of  quasi-flat  terrain.  A  slight  degradation  in  accuracy  may  be  expected. 
The  standard  analytical  upward  continuation  method  and  two  spatial  collocation 
methods  for  the  estimation  of  first  order  derivatives  of  the  anomalous  gravity 
potential  can  be  replaced  by  numerical  upward  continuation  using  Laplace's 
equation  because  of  the  availability  of  three  gravity  vector  components  at 
the  information  base  level  surface.  In  this  case,  the  standard  information 
base  can  be  reduced  considerably,  and  no  matrix  inversions 

are  required.  Analytical  upward  continuation  of  "topographic  noise"  components 
may  also  be  accomplished  separately,  using  a  higher  resolution  grid.  The 
computation  of  "topographic  noise"  effects  is  mathematically  laborious,  can^ 
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however^ be  achieved  by  means  of  high  speed  computers.  The  alternative  in  the 
presence  of  mountainous  terrain  is  to  reduce  flight  traverses  from  300  km  to 
about  100  km,  to  use  highly  accurate  "truth"  data  at  both  ends  of  the 
traverse,  and  to  repeat  single  traverse  surveys. 
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ABSTRACT 

A  three-axis  superconducting  gravity  gradiometer  which  measures  the  three 
in-line  components  of  the  gravity  gradient  tensor  is  under  development  at  the 
University  of  Maryland.  The  instrument  is  being  developed  under  a  NASA  contract 
for  the  purpose  of  precision  gravity  experiments  and  gravity  field  mapping  from 
an  orbiting  platform.  The  design  of  the  gradiometer  employs  a  number  of  recently 
devised  techniques  which  rely  on  certain  properties  of  superconductors  to  obtain 
a  fundamental  noise  level  of  10“^  Eotvos  as  well  as  the  high  degree  of 

stability  necessary  for  such  extreme  sensitivity.  In  order  to  compensate  for 
errors  associated  with  angular  motions,  the  gradiometer  will  be  integrated  with 
a  six-axis  superconducting  accelerometer  which  is  being  developed  at  the 
University  of  Maryland  under  an  AFGL  contract.  The  accelerometer  will  sense  the 
rotational  and  translational  motions  of  the  gradiometer  platform.  The  principle 
and  design  of  the  instruments  will  be  discussed.  In  addition,  the  results  of 
tests  on  the  first  axis  of  the  gradiometer  will  be  presented. 
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TITLE  OF  PAPER:  Development  of  A  Three-Axis  Superconducting  Gravity  Gradiometer 
and  A  Six-Axis  Superconducting  Accelerometer 

SPEAKER:  Martin  Vol  Moody 

QUESTIOMS  AND  COMMEMTS: 

1.  Question :  Charles  F.  Martin 

The  non-linear  effects  of  levitation  at  11.5  Hz  caused  by  "d"  and  "D" 
dimensions  can  be  Ignored  in  space.  Can  you  adjust  adequately  for 
surface  or  near  earth  use? 

Response : 

Yes,  that  distance  can  be  adjusted  and  should  present  no  problem  in  near 
earth  operation. 

2.  Question :  Anthony  R.  Barringer 

Is  the  six  axis  accelerometer  applicable  in  a  1  "g”  environment? 

Response ; 

Yes.  We  plan  to  test  it  in  a  1  "g"  environment.  It  will  levitate  at 

1  "g". 

3.  Question :  James  E.  Fix 

What  Q  (quality  factor)  did  you  achieve? 

Response  : 


Q  is  pressure  limited.  Achieved  Q  =  10^  in  tests  at  Stanford. 


THE  GRAVITATIONAL  MAGNETIC  FIELD  OF  THE  EARTH 


’  AND  THE  POSSIBILITY  OF  MEASURING  IT 
USING  AN  ORBITING  GRAVITY  GRADIOMETER 


by 

Dr.  Bahram  Mashhoon 
Department  of  Physics  and  Astronomy 
University  of  Missouri- Columbia 
Columbia,  MO  65211 


ABSTRACT 

Einstein's  theory  of  gravitation  predicts  that  a  rotari ng  body  such  as  the 
Earth  carries  the  local  inertial  frames  around  it  ("dragging  of  the  inertial 
frames").  This  results  in  an  effective  gravitational  "magnetic"  field  due  to 
mass  current  in  addition  to  the  usual  -  i.e.,  essentially  Newtonian  -  gravita¬ 
tional  "electric"  field,  in  close  analogy  with  the  electromagnetic  field  around 
a  rotating  charged  body.  This  fundamental  prediction  of  general  relativity  has 
not  yet  been  tested.  A  gravity  gradiometer  in  Earth  orbit  measures  mostly  the 
components  of  the  tidal  field  (i.e.,  curvature  tensor)  of  the  Earth.  The 
gravitational  "electric"  and  "magnetic"  fields  of  the  Earth  are  reflected  in 
the  tidal  tensor.  The  tidal  accelerations  due  to  a  rotating  system  have  been 
studied  and  the  relativistic  corrections  to  the  Newtonian  results  have  been 
determined  (Mashhoon  and  Theiss,  Phys.  Rev.  Lett.  49,  1542-L545,  1982).  The 
possibility  of  measuring  relativistic  effects  using  Paik's  superconducting 
gravity  gradiometer  currently  under  development  is  investigated. 
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TITLE  OF  PAPER:  The  Gravitational  Magnetic  Field  of  the  Earth  and  the  Possibility 
of  Measuring  it  Using  an  Orbiting  Gravity  Gradiometer 

SPEAKER:  Bahrain  Mashhoon 

QUESTIONS  AND  COMMENTS: 

1 .  Question:  Dan  Long 

Have  the  gravitational  magnetic  effects  been  derived  classically  as  can  be 
done  in  Electricity  and  Magnetism? 

Response : 

One  would  get  about  the  same  results  but  that  theory  has  defects  and  I  have 
not  used  it- 

2.  Question:  Alan  Zorn 

Your  satellite  examples  show  that  the  relativistic  correction  terras  are 
both  conservative  and  harmonic  (symmetric  and  traceless).  Is  this  true 
in  general? 

Response ; 

Yes,  general  relativity  implies  this  in  general. 

3.  Question:  Charles  Finley 

Does  the  orientation  of  the  spacecraft,  i.e.,  the  orientation  of  the  SGG 
within  the  Spacecraft  (inertial  vs.  earth  pointing)  make  a  difference 
for  this  experiment?  If  so,  which  would  be  preferred? 

Response : 

In  theory,  it  would  not  matter  which  orientation  were  used  since  the 
deviation  from  true  gyroscopic  inertial  orientation  would  have  to  be 
measured  and/or  determined.  However,  in  practice  this  is  an  important 
question  in  that  they  would  want  the  orientation  to  be  such  as  to  minimize 
the  difficulty  in  measuring  the  deviation  from  the  true  Inertial. 

4.  Question:  Warren  Heller 

What  is  the  size  of  the  quadrilateral/magnetic  effect  in  ordinary  terms? 
Response: 


3  X  10~1-1-  of  an  orbit  time. 


TESTS  OF  GENERAL  RELATIVITY  IN  EARTH  ORBIT 
USING  A  SUPERCONDUCTING  GRAVITY  GRADIOMETER 


by 

Dr.  Ho  Jung  Paik 

Department  of  Physlca  and  Astronomy 
University  of  Maryland 
College  Park,  MD  20742 


ABSTRACT 

A  gravity  gradioraeter  measures  a  component  of  the  Riemann  curvature  tensor 
which  is  the  fundamental  gravitational  field  in  General  Relativity.  Both  the 
mass  M  and  the  angular  momentum  J  of  a  spinning  object  contribute  to  the  Riemann 
tensor.  For  the  Earth,  the  relativistic  corrections  due  to  M  and  J  amount  to 
7  X  10~^®  and  8  x  10“^^  of  the  Newtonian  gravity  gradient,  respectively.  These 
effects  could  be  resolved  with  a  signal-to-noi se  ratio  of  10  in  one  year  by  a 
three-axis  superconducting  gravity  gradlometer  of  10"^  E  sensitivity  in  a 

low-altitude  polar  orbit.  We  discuss  the  experimental  strategy  and  the  require¬ 
ments  that  these  experiments  put  on  the  orbit  and  the  altitude  control  of  the 
spacecraft . 
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TITLE  OF  PAyER;  Tests  of  General  Relativity  in  Earth  Orbit  Using  A.  Superconducting 
Gravity  Gradlometer 

SPEAKER:  Ho  Jung  Paik 

QUESTIONS  AND  COMMENTS: 

1.  Question:  Richard  Hansen 

What  about  interference  to  dragging  of  inertial  frames  from  higher  moments 
of  the  earth? 

Response : 


These  appear  at  higher  multiples  of  the  orbital  frequency.  The  only  effects 
which  can  interfere  are  those  with  components  at  orbital  frequency. 

2 .  Question:  Alan  Zorn 

Can  you  comment  on  the  relationship  between  the  current  state  of  scale 
factor  mismatch  that  you  expect  on  your  current  instrument  and  how  the 
requirements  for  null  test  (10“^*^)  and  magnetic  field  (10~^)  experiments  can  be 
met? 

Response ; 


We  expect  lO"^  on  the  current  3-axls  instrument,  so  the  magnetic  field 
experiment  can  be  done  pretty  easily.  We  recognize  that  the  null  test 
experiment  will  be  much  harder.  However,  post  processing  of  the  data  will 
help  us  reach  this  goal,  we  believe. 


MAGNETIC  ISOLATION-aoSING  THE  LOOP 


by 

Dr.  Dave  Sonnabend 
Mr.  A.  Miguel  •‘‘an  Martin 

Jet  Propulsion  Laboratory 
California  Institute  of  Technology 


ABSTRACT 

Progress  on  the  design,  construction,  and  testing  of  the  JPL  single  axis 
magnetic  Isolation  test  facility  has  been  previously  reported.  Our  more  recent 
work  has  been  to  close  the  loop  with  a  digital  computer,  in  order  to  try  out 
various  vibration  isolation  and  semi  drag  free  control  laws.  The  status  of 
this  work  will  be  reported.  For  vibration  isolation,  we  will  discuss  both  the 
theoretical  limitations  of  the  technique,  and  the  test  limitations  imposed  by 
the  facility.  Finally,  for  semi  drag  free  operation,  we  will  look  at  shaping 
of  the  current  impulses,  to  mini-nice  the  excitation  of  vibration  modes  of  the 


floated  Instrument. 


Magnetic  Isolation — Closing  the  Loop 

Dave  Sonnabend  and  A.  Miguel  San  Martin 

Jet  Propulsion  Laboratory 
California  Institute  of  Technology 

S  February  1987 


Abstract 

Progress  on  the  design,  construction,  and  testing  of  the  JPL  single 
axis  magnetic  isolation  test  facility  has  been  previously  reported.  Our 
more  recent  work  has  been  to  close  the  loop  with  a  digital  computer, 
in  order  to  try  out  various  vibration  and  semi  drag  free  control  laws. 
The  status  of  this  work  will  be  reported.  For  vibration  isolation,  we 
will  discuss  both  the  theoretical  limitations  of  the  technique,  and  the 
test  limitations  imposed  by  the  facility.  Finally,  for  semi  drag  free 
operation,  we  will  look  at  shaping  of  the  current  impulses,  to  minimize 
the  e.xcitation  of  vibration  modes  of  the  floated  instrument. 


1  Introduction 

Previous  work  on  magnetic  eddy  current  isolation  of  gradiometers  was  re¬ 
ported  in  Refs.  1  and  2.  A  single  axis  torsion  pendulum  facility  was  de¬ 
scribed,  on  which  vibration  isolation  and  semi  drag  free  operation  ideas 
could  be  tested.  The  facility  characteristics,  and  initial  testing  of  the  mag¬ 
netic  forcers  was  reported  in  Ref.  2.  The  abstract  above  tells  what  we  would 
like  to  have  said  today;  but.  as  we  have  not  made  that  much  progress,  I'll 
discuss  only  what  we’ve  actually  accomplished.  Mostly,  this  will  cover  our 
efforts  to  close  a  digital  control  loop  around  the  magnetic  test  facility.  At 
the  end  I’ll  also  mention  some  parallel  analytical  work,  at  Arizona  State 
University,  on  the  physics  of  eddy  currents. 
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2  Hardware  Description 

The  plant  (the  original  facility)  comprises  the  torsion  pendulum,  four  eddy 
current  actuators,  and  two  sensors.  A  photo,  from  Ref.  1,  is  shown  in  Fig.  1. 
The  actuators  are  mounted  externally,  on  brackets,  positioned  at  opposite 
faces  of  the  box  (where  the  gradiometer  would  reside),  two  on  each  side.  The 
coil  axes  are  normal  to  the  box  faces,  along  the  direction  of  the  box  motion. 
The  design  and  implementation  of  the  torsion  pendulum  and  actuators  are 
described  in  Ref.  2.  A  condensed  theoretical  treatment  of  eddy  current 
forcing  is  given  in  Ref.  1;  and  a  substantial  generalization  in  Ref.  3.  A 
synopsis  of  the  latter  work  is  presented  here,  in  the  final  section. 

The  two  sensors  measure  the  displacement  between  the  box  and  the  brack¬ 
ets.  They  are  Kaman  Instrument  Co.  eddy  current  sensors,  Model  KD-2400. 
Each  sensor  consists  of  two  subassemblies — the  sensor  head  and  a  signal  con¬ 
ditioning  module.  The  proximity  of  the  sensor  to  the  box  controls  a  variable 
gain  oscillator  within  the  conditioning  module.  The  oscillator  amplitude  is 
detected  to  provide  an  analog  signal  proportional  to  displacement.  Among 
the  sensor  features  are  low  cost,  no  contact,  10  kHz  frequency  response, 
adjustable  gauge  factor,  and  0.25  mm  resolution. 

The  function  of  the  control  loop  is  to  execute  a  control  law,  designed  to 
achieve  certain  performance  objectives.  The  controller  samples  analog  data 
from  the  sensors,  converts  these  to  digital  form,  computes  actuator  com¬ 
mands  according  to  the  trial  control  law.  converts  back  to  analog,  and  fi¬ 
nally  doles  out  power  to  the  actuators.  All  these  operations  are  carried 
out  periodically,  in  real  time.  The  two  most  important  requirements  of  this 
hardware  are  at  least  12  bit  resolution,  and  a  sampling  frequency  of  at  least 
10  Hz. 

Fig.  2  shows  the  block  diagram  of  the  control  loop.  The  computer  is  a  PC’s 
Limited  IBM  AT  compatible;  and  its  characteristics  are  listed  in  Table  1. 
The  function  of  the  computer  is  to  carry  out  the  calculations  dictated  by 
the  control  law,  control  the  overall  timing  of  the  loop,  and  generate  test 
evaluation  data. 

The  input/output  interface  board  is  located  in  one  of  the  expansion  slots  of 
the  computer,  and  is  connected  to  the  computer  output  bus.  Its  function  is 
to  perform  A/D  and  D/A  input  and  output  conversions.  Table  2  lists  the 
functions  of  the  interface  board.  The  board  consists  of  a  multiplexer,  one 
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A/D  and  two  D/A  converters,  a  clock,  <utd  a  microprocessor  that  controls 
the  operation  of  these  elements  and  communicates  with  the  computer  mi¬ 
croprocessor  in  a  high  level  command  language.  The  analog  outputs  are  of 
the  sample  and  hold  type.  Purely  digital  input  and  output  lines  are  also 
provided. 

The  signal  generator  provides  an  analog  sinusoidal  constant  amplitude  and 
frequency  signal.  We  are  presently  operating  at  40  kHz.  The  amplitude 
modulator  regulates  the  amplitude  of  the  sinusoid  according  to  the  control 
signal  from  the  interface  board.  At  the  time  of  the  talk,  this  was  a  multiply¬ 
ing  D/A  converter  (see  Fig.  3);  requiring  a  digital  output  from  the  interface 
board.  However,  due  to  the  unavailability  of  the  DAC1220  micro-chip,  this 
approach  was  changed  to  using  an  analog  modulator,  as  shown  in  Fig.  4. 

Finally,  the  function  of  the  power  amplifier  is  to  amplify  the  signal  from  the 
amplitude  modulator  to  provide  the  desired  current  to  the  forcing  coils.  It 
is  a  standard  stereo  audio  amplifier,  a  Phase  Linear  Model  400,  Series  Two. 
It  can  put  out  205  watts  rms  into  8  ft  loads,  and  has  a  bandwidth  of  50 
kHz.  Since  the  talk,  these  loops  have  all  been  closed;  and  simple  control 
laws  have  been  demonstrated. 


3  Eddy  Current  Analysis 


In  a  parallel  effort  to  achieve  a  better  understanding  of  the  physics  of  eddy 
current  interactions,  a  contract  was  let  to  Arizona  State  University.  Ear¬ 
lier  work,  summarized  in  Ref.  1,  examined  a  conductive  spherical  shell  on 
the  axis  of  a  circular  coil,  carrying  a  sinusoidal  current.  The  force  and 
dissipation  were  calculated.  The  newer  work,  reported  in  Ref.  3,  extends 
this  to  off-axis  cases,  and  computes  the  torque  on  the  sphere.  The  analyt¬ 
ical  results  are  much  too  complicated  to  be  presented  here;  however,  the 
qualitative  results  are  summarized  in  Fig.  5.  As  might  be  expected  from 
symmetry  considerations,  the  axial  repelling  force  does  not  depend  on  dis¬ 
placement  from  the  axis,  to  first  order  in  this  displacement;  and  there  is  no 
torque  on-axis.  However,  off-axis,  to  first  order  in  the  displacement,  there 
is  a  force  tending  to  restore  the  sphere  to  the  axis,  and  a  torque  tending  to 
roll  the  sphere  toward  the  coil,  both  proportional  to  the  displacement.  The 
dissipation  in  the  sphere  was  shown  not  to  vary,  to  first  order  in  displace¬ 
ment.  Future  work  is  planned  in  this  area  to  develop  an  electrical  model  of 


3 


the  coil-sphere  combination,  that  will  serve  as  a  basis  for  design  of  practical 
isolation  systems. 
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OPERATING  SYSTEM;  DOS 
PROGRAMMING  LANGUAGE:  STSC  APL<PLUS 


JPL  INPUT/OUTPUT  INTERFACE  BOARD 


Fig.  1.  Magnetic  Isolation  Test  Facility 


JPL  SYSTEM  BLOCK  DIAGRAM 


Figure 


PERFORMS  THE  FOLLOWING  PRODUCT; 


J 

& 

•  “I 


WHERE  AN  =  1  IFAN  DIGITAL  INPUT  IS  HIGH 
AN  =  0  IF  AN  DIGITAL  INPUT  IS  LOW 


Figure  4. 
Analog  Modulator 


EDDY  CURRENT  ANALYSIS 

A 

\J 

COIL  SPHERE 

Force:  Zero  Order  -  Away  from  Coil 
First  Order  -  Toward  Axis 

Torque:  Zero  Order  -  None 

First  Order  -  Roll  Toward  Coil 

Dissipation:  No  First  Order 
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Qienoy,  WA  99004 


ABSTRACT 

I  report  on  progress  in  the  three  areas  of  vibrations,  the  optical  lever, 
and  the  isothermal  environment.  All  vibrating  machinery  has  been  decoupled  from 
the  laboratory  floor  and  supported  from  the  second  floor  above.  The  apparatus 
table  has  been  rigidified  against  low  frequency  vibrations.  Preliminary 
measurements  indicate  a  reduction  of  about  a  factor  of  4  of  unwanted  vibrations 
coupling  into  the  torsion  pendulum.  The  new  optical  lever  is  now  fully  opera¬ 
tional.  It  appears  to  have  a  5  part  per  ten  thousand  per  hour  drift  which 
cannot  be  traced  to  a  particular  cause  such  as  the  light  source  power  supply. 

I  suspect  long  term  deformation  of  the  light  bulb  filament  Is  taking  place. 

This  drift  should  present  no  measurement  problem,  but  is  disappointing.  Cbn- 
structlon  of  a  new  (fireproof)  isottiermal  insulation  jacket  for  the  apparatus 
is  now  in  progress  along  with  the  ci)nstructlon  of  the  thermal  regulation  devlci-. 
I  hope  to  report  on  a  fully  operational  signal  to  noise  ratio  at  the  conference. 
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I  report  on  progress  In  overhauling  and  Improving  the 

Washington  University  Cavendish  Torsion  Balance.  Recent 
Improvements,  mainly  funded  by  AFGL,  are  Illustrated  In  the 
following  slides  presented  at  th^conference.  These  substantial 
modifications  addressed  reducing  the  effects  of  tilt  on  the 
output  signal,  reducing  vibrations,  and  improving  thermal 
stability.  The  last  illustration  shows  the  gravitational  force 
signal  now  put  out  by  the  apparatus.  It  Indicates  about  a  2  part 
per  ten  thousand  torque  signal  stability.  This  is  about  a  factor 
of  10  Improvement  over  previous  work. 


View  (1)  The  entire  apparatus  with  ring  masses  in  view.  Note 
the  thin  aluminum  legs  which  are  reinforced  in  the  next  picture. 


view  (2)  3  X  3  X  1/2  alunlnum  angle  has  been  used 

the  frame.  The  1/2”  aluminum  plate  below  the  main 


to  reinforce 
table  was 


later  loaded  with  200  lbs  of  lead 


View  (3) 


A  view  of  the  frame  showing  the  cross  reinforcing. 


View  (5)  The  new,  tilt  proof,  optical  lever  of  quasi  auto 
collimator  design  can  be  seen  in  the  right  portion  of  the  photo 


View  ( 6 ) 
external 


The  new  optical  lever  features  a  newly  constructed 
light  source  and  fiber  optics.  This  reduces  temperature 


drllL  In  the  vacuum  chamber. 


view  (7)  The  fan  for  cooling  the  optical  lever  light  source  is 
mounted  on  the  I  beams  of  the  second  floor  above  to  provide 
vibrational  isolation.  The  4  x  4  in  the  picture  can  exert  a  few 
tons  against  the  second  floor  to  tilt  the  floor  and  apparatus  at 
will,  without  introducing  further  gravitational  masses  near  the 


apparatus . 


view  (8)  The  newly  constructed,  thermally,  and  electrically 
shielded  output  box  for  the  optical  lever. 


view  (10)  After  leaving  the 
fore-line  hose  goes  to  a  400 


500  lbs  Inertial  block  mass,  the 
lbs  mass  mounted  on  lead  brlcks< 


view  (11)  The  newly  constructed  Isothermal  water  bath  needed  to 
supply  the  diffusion  pump  with  constant  temperature  water  has 
several  vibrating  components.  These  are  mounted  on  the  angle 
Iron  seen  hanging  down. 


View  (14)  Heater  and  control  thermlster  with  fail-safe 
thermostat  mounted  on  the  North-East  leg. 


View  (15)  Heater  and  control  thermlster  with  fail  safe 
thermostat  mounted  on  the  West  leg. 


view  (16) 


The  five  O.Ol  C  Yellow  Springs  Instrument  Company 


thermal  controllers  used  to  control  power  to  the  heaters. 


fitser  glass  felled  insulation. 


View  (18)  Insulation  of  the  sphere  top.  An  additional 
hood  covers  what  Is  seen.  Preliminary  results  suggest  0 
stability  over  a  few  hours. 


6  inch 
-00l°C 


view  (19)  A  thermlster  shown  on  Che  side  of  the  attracted  ball 
chamber.  This  therralster  will  monitor  any  changes  in  temperature 
produced  by  introducing  Che  attracting  masses  near  the  chamber. 

18  chermlsters  have  been  mounted  on  the  apparatus. 


View  (20)  Strip  chart  trace  of  the  optical  lever  output  of  the 
gravitational  signal.  10  intervals  is  1%  of  the  ring  torque 
signal.  The  centroid  of  the  sine  wave  seems  stable  to  2  or  3 


parts  per  ten  thousand  of  the  ring  torque. 


